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Theoretical  and  experimental  studies  of  many-body  effects  in  semiconductor  quan- 
tum well  (QW)  lasers  are  presented.  These  many-body  effects  arise  bom  the  Coulomb 
force  interactions  among  the  charged  particles  trapped  in  the  quantum  well  and  their  inter- 
action with  phonons.  Expressions  are  derived  for  the  dependence  of  optical  gain  and  spon- 
taneous emission  on  photon  energy  using  an  electron-hole  plasma  (EHP)  theory  with 
many-body  effects  taken  into  account.  Optical  gain  and  spontaneous  emission  spectra  are 
calculated  for  the  well-developed,  infrared  emitting,  InGaAs/GaAs  strained  QW  structure 
and  the  relatively-new  blue-green  CdZnSe/ZnSSe  strained  QW  structure.  The  results  show 
that  Coulomb  attraction  effects  (CE)  have  a  significandy  stronger  effect  on  the  spectral 
characteristics  of  the  wide-bandgap  CdZnSe  QW  structure  than  on  the  narrow-bandgap  In- 
GaAs  QW  structure.  The  room-temperature  dependences  of  lasing  energy  and  threshold 
current  on  cavity  length  for  InGaAs  and  CdZnSe  QW  diode  lasers  are  numerically  simulat- 
ed and  compared  to  experiment.  In  all  cases,  the  comparisons  are  found  to  be  in  better 
agreement  when  CE  is  included  in  the  calculations.  The  importance  of  CE  is  also  found 


when  comparing  the  experimental  and  calculated  lasing  energies  for  a  CdZnSe  QW  laser 
over  the  temperature  range  from  160  K  to  300  K.  For  photon  energies  somewhat  higher 
than  the  lasing  energy,  the  calculated  optical  gain  spectra  with  CE  included  show  reso- 
nance-like absorption  features  usually  attributed  to  excitons  in  the  QW.  A  comparison  of 
these  absorption  spectra  with  published  experimental  spectra  shows  good  agreement. 
Since  excitons  do  not  exist  in  the  CE  model  used,  this  success  refutes  published  work  which 
says  that  free  excitons  play  an  important  role  in  determining  the  magnitude  of  optical  gain 
in  room  temperature  CdZnSe  QW  lasers. 
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CHAPTER  1 
INTRODUCTION 


Semiconductor  diode  lasers  have  become  the  key  components  in  many  commercial 
products  such  as  compact  disk  players,  laser  printers  and  optical  fiber  telecommunication 
systems.  High-power  diode  lasers  are  also  used  in  applications  such  as  medical  surgery, 
free-space  communication  and  diode  laser  pumped  solid-state  lasers.  Diode  laser  design 
with  a  double  heterojunction  (DH)  active  layer  of  thickness  about  100  to  200  nm  was  per- 
fected in  the  late  1970s  and  early  1980s.  When  reducing  the  thickness  of  the  active  layer 
down  to  the  10  nm  regime,  electrons  and  holes  are  confined  in  the  active  layer  and  the  as- 
sociated states  become  quantized  in  the  potential  wells  along  the  direction  normal  to  the  ac- 
tive layer.  Quantum  well  (QW)  diode  lasers,  devices  with  this  very  thin  active  layer  design, 
were  realized  in  the  1980s  and  shown  to  have  higher  performance  than  their  DH  counter- 
parts. By  introducing  built-in  strain  in  the  QW,  these  strained  QW  lasers  were  found  to  be 
substantially  superior  to  unstrained  QW  lasers.  Although  these  improvements  in  diode  la- 
sers rely  heavily  on  advances  in  high-quality  crystal-growth  techniques  such  as  molecular 
beam  epitaxy  (MBE),  the  ideas  behind  them  came  from  thoroughly  understanding  the  basic 
mechanism  of  optical  gain  and  spontaneous  emission  in  diode  lasers.  Stimulated  electron- 
hole  recombination  in  the  active  layer  was  believed  to  be  responsible  for  the  laser  action  in 
diode  lasers  and  an  electron-hole  plasma  (EHP)  theory  was  used  to  model  the  laser  devices 
such  as  the  well-developed  infrared  (IR)  GaAs-based  QW  lasers  [Asa93,  Col95,  Coi93]. 
Soon  after  the  demonstration  of  the  world's  first  blue-green  diode  lasers  fabricated  using  a 
CdZnSe  QW  structure  by  3M  Company  personnel  in  June  1991  [Haa91],  confusion  about 
the  gain  mechanism  in  these  new  CdZnSe  QW  lasers  began.  J.  Ding  et  al.  reported  the  first 


observation  of  so-called  excitonic  gain  and  laser  emission  in  ZnSe-based  QWs  and  pro- 
posed a  phenomenological  model  of  optical  gain  based  on  partial  phase-space  filling  (PSF) 
of  an  inhomogeneously  broadened  exciton  resonance  [Din92].  This  phenomenological 
gain  model,  which  includes  the  effects  of  exciton-phonon  interaction,  was  used  to  explain 
their  pump-probe  experiment  at  cryogenic  temperature  up  to  220  K  Pin93].  In  1994,  J. 
Ding  et  al.  reported  the  measured  gain/absorption  spectra  of  CdZnSe  QW  diode  lasers  at  T 
=  300  K  in  which  two  resonance-like  absorption  peaks  appearing  at  low  excitation  level 
were  attributed  to  absorption  resonances  from  n  =  \  HH  and  LH  excitons  [Din94].  The  con- 
clusion drawn  was  that  the  origin  of  optical  gain  in  blue-green  CdZnSe  QW  lasers  is  exci- 
tonic in  nature  even  at  room  temperature.  Since  the  conventional  EHP  theory  did  not 
include  exciton  effects,  these  claims  attracted  considerable  attention  in  the  diode  laser  com- 
munity. 

Setting  aside  the  truth  about  the  lasing  mechanism  in  CdZnSe  QW  lasers  for  a  mo- 
ment, physicists,  working  for  many  years  in  the  research  field  of  fundamental  semiconduc- 
tor physics,  have  been  studying  the  mutual  interactions  among  charged  particles  (i.e. 
electrons  and  holes)  in  semiconductors.  These  interactions  in  an  EHP  system  include  car- 
rier scattering  (CS),  bandgap  renormalization  (BGR)  and  Coulomb  enhancement  (CE)  and 
are  usually  referred  to  as  many-body  effects  in  the  literatures.  Haug  et  al.  showed  that  gain 
spectra  of  a  semiconductor  changed  when  many-body  effects  were  taken  into  account 
[Hau90].  Although  CS  and  BGR  were  usually  included  in  the  conventional  EHP  theory 
[Asa93,  Col95,  Cor93],  there  seemed  to  be  a  tendency  to  ignore  CE  in  the  diode  laser  com- 
munity since  CE  calculations  are  more  complicated  than  those  for  CS  and  BGR.  A  theory 
with  such  an  ignorance  of  CE  was  shown  to  lead  to  predictions  that  were  in  worse  agree- 
ment with  experiment  than  predictions  of  a  theory  with  only  CS  included  [Che93,  Chi88]. 
The  errors  result  from  the  fact  that  it  is  inconsistent  to  account  for  the  Coulomb  repulsion, 
which  gives  rise  to  BGR,  and  at  the  same  time  ignore  the  Coulomb  attraction,  which  gives 
rise  to  CE.   Since  Coulomb  interactions  are  stronger  in  materials  with  smaller  dielectric 


constants,  CE  is  expected  to  be  more  important  in  the  wide-bandgap  CdZnSe  QW  lasers 
[Cho95b,  Ree95a]  than  in  the  narrow-bandgap  GaAs-based  QW  lasers  [Cho94,  Hau90]. 

Now  going  back  to  the  discussion  about  the  nature  of  the  gain  mechanism  in  CdZnSe 
QW  lasers,  two  questions  were  thus  raised:  "Can  the  operation  of  CdZnSe  QW  lasers  be 
simulated  using  EHP  theory  with  inclusion  of  CE?"  and  "Is  optical  gain  in  CdZnSe  QW 
lasers  of  excitonic  nature?"  The  approach  to  answering  these  questions  is  twofold:  theoret- 
ical and  experimental.  In  the  theoretical  part,  theories  of  many-body  effects  available  in 
the  literatures  were  studied  extensively  and  understood  thoroughly.  These  existing  theoret- 
ical works  were  all  done  under  the  assumption  of  a  perfect  two-dimensional  (2D)  QW  with 
zero  QW  thickness  and  infinite  barrier  height.  In  addition,  the  formula  derived  for  many- 
body  effects  are  limited  to  the  case  that  only  the  lowest  conduction  (CB)  and  valence  band 
(VB)  subbands  exist.  The  simple  perfect  2D  QW  model  has  no  physical  meaning  except 
for  easy  demonstration  of  many-body  effects  in  semiconductors. 

In  order  to  use  these  many-body  formulas  and  implement  them  in  the  device  simula- 
tion for  actual  QW  diode  lasers,  several  modifications  to  the  formulas  had  to  be  made. 
First,  the  finite  QW  thickness  and  barrier  height  in  real  QW  lasers  were  taken  into  account 
In  such  a  quasi- 2D  QW,  a  number  of  subbands  in  both  CB  and  VB  become  available  from 
solving  the  one-dimensional  (ID)  Schrodinger's  equation.  Second,  effective  mass  anisot- 
ropy  in  the  VB  due  to  the  large  built-in  strain  the  QW  was  considered,  where  the  heavy  hole 
(HH)  and  light  hole  (LH)  effective  masses  in  the  in-plane  direction  become  different  from 
those  in  the  growth  direction.  Third,  a  model  with  CE  included  for  the  spontaneous  emis- 
sion function  was  developed  in  order  to  estimate  the  threshold  current  required  for  lasing, 
since  most  of  the  existing  literatures  on  many-body  effects  really  dealt  only  with  the  gain 
function  [Cho94,  Hau90].  Four,  the  modified  many-body  formulas  were  integrated  with 
the  state-of-the-art  EHP  tiieory  where  CE  has  not  been  taken  into  account  [Asa93,  Col95, 
Cor93].  Efforts  were  also  placed  on  converting  formulas  from  wavevector  space  to  energy 
space,  from  summations  to  integrals  and  from  the  C.G.S.  units  used  by  many-body  physi- 


cists  to  the  M.K.S.  units  which  are  more  familiar  to  diode  laser  engineers.  As  a  result  of 
these  works,  the  complicated  formulas  derived  by  many-body  physicists  have  been  trans- 
lated to  a  form  more  user-friendly  so  that  electrical  engineers  can  implement  them  into 
computer  codes  for  device  simulation  of  QW  diode  lasers. 

In  the  experimental  part  of  this  work,  a  number  of  well-developed  IR  InGaAs 
strained  QW  lasers  and  relatively  new  blue-green  CdZnSe  strained  QW  lasers  were  fabri- 
cated. Since  the  laser  threshold  current  and  peak  wavelength  change  with  laser  cavity 
length,  different  models  for  the  physical  mechanism  of  optical  gain  and  spontaneous  emis- 
sion were  tested  by  investigating  these  length  dependences.  Since  the  emphasis  in  this 
work  has  been  placed  on  the  CdZnSe  QW  laser,  further  measurements  on  the  temperature 
dependence  of  lasing  wavelength  were  performed  on  the  CdZnSe  QW  lasers. 

This  dissertation  is  organized  such  that  the  theoretical  models  are  described  and  de- 
veloped first,  experiments  and  experimental  results  are  then  presented,  predictions  by  using 
the  theoretical  models  and  their  comparisons  with  experiments  are  made  and  finally  the 
conclusions  are  drawn.  Several  comparisons  and  comments  on  the  differences  between  the 
InGaAs  and  CdZnSe  QW  lasers  are  made  throughout  each  chapter. 

In  Chapter  2,  the  gain  and  spontaneous  emission  functions  based  on  a  single-particle 
model  are  discussed  in  detail.  It  begins  by  introducing  the  epitaxial  structures  for  the  two 
QW  lasers  of  interest:  III-V  InGaAs  and  II-VI  CdZnSe  QW  lasers.  Effect  of  strain  on  the 
bandgap  of  the  QW  material  are  discussed  and  the  active  region  model  used  in  the  theoret- 
ical simulation  for  the  QW  laser  structures  described.  Finally,  the  optical  gain  and  sponta- 
neous emission  spectra  in  the  single-particle  model  are  calculated. 

In  Chapter  3,  many-body  effects,  such  as  CS,  BGR  and  CE,  on  the  gain  and  sponta- 
neous emission  functions  are  discussed  sequentially.  Semi-derivations  of  BGR  and  CE  are 
presented.  Calculation  results  of  the  optical  gain  and  spontaneous  emission  spectra  with 
many-body  effects  included  are  shown  and  comments  on  how  CS,  BGR  and  CE  modify  the 
single-particle  spectra  are  made.  > 


In  Chapter  4,  experimental  measurements  of  laser  threshold  characteristics  for  In- 
GaAs/GaAs  QW  and  CdZnSe/ZnSSe  QW  lasers  are  presented  and  compared  to  predic- 
tions. The  first  question  mentioned  above  "Can  the  operations  of  CdZnSe  QW  lasers  be 
simulated  by  using  the  EHP  theory  with  inclusion  of  CE?"  is  answered. 

In  Chapter  5,  the  other  question  "Is  optical  gain  in  CdZnSe  QW  lasers  of  excitonic 
nature?"  is  answered  by  comparing  several  spectral  features  of  the  calculated  gain/absorp- 
tion and  spontaneous  emission  spectra  with  the  inclusion  of  CE  to  those  of  the  experimental 
observations. 

In  Chapter  6,  the  carrier-induced  antiguiding  factor  is  calculated  for  the  CdZnSe  QW 
lasers  and  its  effect  on  the  quality  of  laser  beam  is  described.  Near-field  patterns  are  mea- 
siu'ed  and  compared  qualitatively  to  the  predictions. 

Finally,  In  Chapter  7,  conclusions  are  drawn  and  future  works  are  suggested. 


CHAPTER  2 
SINGLE-PARTICLE  MODEL 


2.1  Introduction 


The  chapter  begins  by  introducing  the  epitaxial  structures  for  the  two  QW  lasers  of 
interest:  ni-V  InGaAs  and  II-VI  CdZnSe  QW  lasers.  Due  to  the  mismatch  of  the  in-plane 
lattice  constant  between  a  thin  QW  layer  and  thick  barrier  layers,  the  thin  QW  layer  is  sub- 
ject to  lattice  deformation  and  under  biaxial  strain.  Then,  effect  of  strain  on  the  bandgap  of 
the  QW  material  will  be  discussed.  A  common  approach  is  used  to  determine  the  bandgap 
lineup  at  the  interface  between  the  QW  and  barrier  layers  by  assuming  a  conduction  band 
offset  ratio.  Strain  calculation  will  be  performed  on  both  InGaAs  and  CdZnSe  QW  laser 
structures.  Next,  we  consider  an  active  region  model  to  be  used  in  the  theoretical  simula- 
tion for  the  QW  laser  structures.  It  consists  of  a  QW  layer  and  two  barrier  layers  of  infmite 
thickness  on  each  side.  A  flat-band  approximation  is  assumed  in  constructing  the  potential 
well  diagram  for  the  active  region  model.  A  typical  quantum  well  problem  is  then  solved 
for  the  quantized  energy  levels  in  the  CB  and  VB  (including  HH  and  LH  bands)  potential 
wells.  Numerical  results  from  the  calculation  for  the  two  material  systems  will  be  shown 
and  summarized.  Finally,  from  the  active  region  model  developed,  we  will  be  able  to  ex- 
plore the  optical  gain  and  spontaneous  emission  functions  using  the  single-particle  model. 
A  set  of  calculated  gain  and  spontaneous  emission  spectra  as  a  function  of  carrier  density 
in  the  QW  for  both  InGaAs  and  CdZnSe  QW  lasers  will  be  presented.  Comparison  between 
two  material  systems  will  be  made. 

Several  appendices  will  serve  as  the  supplements  for  more  in-depth  discussions  of 


many  relevant  subjects  in  this  chapter.  In  brief,  derivation  of  material  parameters  for 
CdZnSe  and  ZnSSe  ternary  alloys  from  the  material  parameters  of  the  related  binary  alloys, 
based  on  an  interpolation  scheme,  can  be  found  in  Appendix  A.  In  Appendix  B,  the  tran- 
sition matrix  element,  which  determines  the  strength  of  interaction  between  two  states,  and 
its  polarization-dependent  effects  are  discussed. 

2.2  Laser  Structures 

2.2.1  InGaAs  OW  Laser 

The  InGaAs/GaAs  graded-index  separate-confinement  heterostructure  (GRIN-SCH) 
single  QW  diode  laser  structure  studied  in  this  work,  for  which  the  schematic  sketch  of 
cross-section  epitaxial  layers  and  the  corresponding  diagram  of  Al  and  In  composition  pro- 
file are  shown  in  Fig.  2.1,  consists  of  a  25  nm  graded  Al^Gai.^As  layer  (x  =  0.05  -  0.6) 
grown  on  a  n-GaAs  substrate,  a  1400  nm  n-Alo.6Gao.4As  cladding  layer,  a  200  nm  graded 
Alj,Ga,.yAs  (undoped)  guiding  layer  (y  =  0.6  -  0.3),  an  8  nm  undoped  In^Gai.^As  (z  -  0.15) 
compressively-strained  quantum  well  centered  between  two  7  nm  GaAs  layers,  a  200  nm 
graded  AlyGa,.j,As  (undoped)  guiding  layer  (y  =  0.3  -  0.6),  a  1300  nm  p-Alo6Gao.4As  clad- 
ding layer,  a  25  nm  graded  Al^Ga,.:(As  layer  (x  =  0.6  -  0.05)  and  a  100  nm  p*-GaAs  cap 
layer.  The  guiding  and  cladding  layers  provide  for  electrical  and  optical  confinement  re- 
spectively in  the  growth  direction. 

2.2.2  CdZnSe  OW  Laser 

A  cross-section  diagram  of  the  epitaxial  layers  of  the  CdZnSe/ZnSSe  SCH  SQW  di- 
ode laser  studied  in  this  work  is  shown  in  Fig.  2.2.  It  consists  of  a  GaAs  buffer  layer  grown 
by  MBE  on  an  n-type  GaAs  (001)  substrate  [Gai93],  a  ZnSe  buffer  layer,  a  quaternary  n- 


ZnMgSSe  cladding  layer,  a  4  nm  Cd^Zrii.jSe  (x  ~  0.3)  compressively-strained  quantum 
well  (QW)  centered  between  n-  and  p-ZnSo.06Seo.94  guiding  layers,  a  p-ZnMgSSe  cladding 
layer  and  a  graded  p*-ZnSe,Tei.,  (x  =  1  to  0)  contact  layer.  This  contact  scheme  provides 
a  reasonably  low  resistance  p-contact  to  the  laser  structure  [Fan92]. 

2.3  Strained  Quantum  Wells 

2.3.1  Strained  Bandgaps 

When  a  very  thin  layer  of  QW  is  grown  out  of  a  material  with  a  larger  native  lattice 
constant  than  that  of  the  surrounding  barrier  layers,  the  lattice  compresses  in  the  plane  of 
the  well  (in-plane  direction  denoted  as  II)  to  match  that  of  the  barrier  layers  and  elongates 
in  the  growth  direction  (normal  to  the  plane  of  well  and  denoted  as  ±)  to  keep  the  volume 
of  each  unit  cell  the  same,  as  shown  in  Fig.  2.3.  In  this  case,  QW  is  said  to  be  under  biaxial 
compressive  strain.  If  the  native  lattice  constant  of  QW  is  smaller  than  that  of  the  barrier 
layers,  QW  is  then  under  biaxial  tensile  strain.  Degree  of  lattice-mismatch  can  be  described 
with  a  lattice-mismatch  parameter  e  defined  as  follows  [Col95]: 

e  =  l--^.  (2.1) 

where  «;,  is  the  in-plane  lattice  constant  of  the  barrier  material  and  a„  is  the  native  un- 
strained in-plane  lattice  constant  of  the  strained  QW  material.  From  Eq.  (2.1),  we  have  e  > 
0  for  biaxial  compressive  strain  (a„  >  aj)  and  e  <  0  for  biaxial  tensile  strain  (a„  <  Qj,).  The 
lattice  distortion  of  the  crystal  due  to  the  stress  resulting  from  lattice-mismatch  is  usually 
defined  mathematically  by  a  strain  tensor  E;,.  Assuming  the  shear  component  of  stress  (i  * 
J)  can  be  ignored,  which  is  valid  in  typical  semiconductor  application,  we  thus  only  need  to 
consider  the  three  diagonal  components  of  the  strain  tensor:  e„,  Eyy  and  e,j,  where  x  and  y 
subscripts  indicate  the  in-plane  directions  (II)  and  z  subscript  is  the  growth  direction  (X). 
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Figure  2.1  Schematic  diagrams  of  the  cross-section  epitaxial  layers  and  the  corre- 
sponding aluminum  and  indium  composition  profile  for  the  InoisGao.gsAs/ 
GaAs  GRIN-SCH  single  QW  laser. 
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Figure  2.2        Schematic  diagram  of  the  cross-section  epitaxial  layers  for  the  CdojZrio  7Se/ 
ZnSo.o6Seo.94  SCH  single  QW  laser. 
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By  symmetry,  the  strain  in  both  x  and  y  directions  must  be  equal  and  we  obtain 

^\\-^.x  =  h,  = -^-  (2.2) 

Using  Eq.  (2. 1)  in  Eq.  (2.2),  we  find  that  e„,  e„  <  0  for  biaxial  compressive  strain  and  e„ 
e^  >  0  for  biaxial  tensile  strain.  With  no  stress  applied  to  the  z  direction  and  the  cubic  sym- 
metry of  the  crystal,  the  strain  in  the  growth  direction  is  related  to  the  strain  in  the  in-plane 
direction  by  [Col95] 

ex-e.z  =  -2F^e||.  (2.3) 

where  Cn  and  C^  are  referred  to  as  the  Young's  elastic  moduli  (usually  described  in  the 
units  of  10'°  N/m^).  Since  Cn  and  C^  are  both  positive  in  common  semiconductors  (usu- 
ally Cii  >  Cn),  the  lattice  deformation  along  ±  direction  will  be  opposite  to  the  deformation 
along  II  directions,  as  depicted  in  Fig.  2.3. 

Because  die  bandgap  of  a  semiconductor  material  is  related  to  its  lattice  spacing,  the 
distortion  in  crystal  lattice  due  to  the  strain  should  lead  to  modifications  in  the  bandgap  of 
the  strained  QW  layer.  Putting  aside  the  bandgap  changes  due  to  the  quantum  confinement 
for  the  present  moment,  which  is  left  to  be  discussed  in  Section  2.4.1,  there  are  two  modi- 
fications in  die  strained  bandgap.  The  first  modification  originates  from  the  hydrostatic 
component  of  the  strain  and  under  biaxial  compressive  strain,  it  produces  an  upward  shift 
in  die  conduction  band  (CB)  as  well  as  a  downward  shift  in  both  heavy  hole  (HH)  and  light 
hole  (LH)  valence  bands  (VB),  as  sketched  in  Fig.  2.4(a).  All  the  energy  shifts  are  in  die 
opposite  directions  under  biaxial  tensile  strain.  As  a  result  of  the  hydrostatic  strain,  the 
bandgap  is  changed  by  an  amount  8e;/  which  is  given  by  [Col95] 

8e„  =  2az,^^f^,  (2.4) 

where  a  is  the  hydrostatic  deformation  potential  and  a  <  0  in  common  semiconductors.  We 
find  that  the  bandgap  is  increased  (Se,,  >  0)  for  biaxial  compressive  strain  and  decreased 
(SE;/  <  0)  for  biaxial  tensile  strain.  The  second  modification  originates  from  the  shear  com- 
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Figure  2.3  Illustration  of  the  crystal  lattice  deformation  resulting  from  the  epitaxial 
growth  of  a  thin  layer  of  QW  with  a  native  lattice  constant  a„  between  two 
thick  barrier  layers  with  lattice  constant  (a)  Oj,  <  a„  (biaxial  compressive 
strain);  (b)  aj  >  a„  (biaxial  tensile  strain). 


ponent  of  the  strain  (the  shear  strain  should  not  be  confused  with  the  shear  stress  which  is 
zero  in  this  case  [Col95])  and  separates  the  LH  band  from  the  HH  band,  as  sketched  in  Fig. 
2.4(a).  Under  biaxial  compressive  strain,  the  shear  strain  produces  an  upward  shift  in  HH 
band  by  an  amount  Ses,  a*  ^d  ^  downward  shift  in  LH  band  by  an  amount  8e^  ,;,  which  can 
be  expressed  as 

8e..*A  =  ^E|i^"c^^".  (2-5) 

5'=5,M  =  5es.AA(l-^0-  (2.6) 

where  b  is  the  shear  deformation  potential  and  ii<  0  in  common  semiconductors  and  A  is 
the  spin-orbit  energy  which  separates  the  split-off  (SO)  band  from  both  the  HH  and  LH 
bands.  Under  biaxial  tensile  strain,  the  shear  strain  produces  an  downward  shift  in  HH  band 
(Sej,  wi  <  0)  and  an  upward  shift  in  LH  band  (Sej.  ,*  <  0).  As  a  consequence  of  the  shear 
strain,  the  band  edge  degeneracy  of  the  HH  and  LH  bands  is  removed  and  two  strained 
bandgaps  are  now  needed  to  be  defined  as  follows: 

£^(HH)  =£^(bulk)+5e^-5e,_,,  ,-r^.^ 
Eg  (LH)  =  Eg  (bulk)  +  5e„  +  8es,  ,^ ' 

where  Eg  (bulk)  is  the  unstrained  bulk  bandgap,  E^  (HH)  is  the  C-HH  strained  bandgap  and 
Eg  (LH)  is  the  C-LH  strained  bandgap.  The  splitting  of  the  HH  and  LH  bands  is  usually 
quantified  by  the  splitting  energy  5  which  can  be  defined  as 
5^£^(LH)-£^(HH) 

where  Eqs.  (2.7)  and  (2.6)  are  utilized  in  the  last  two  expressions,  respectively. 
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Figure  2.4       Effects  of  the  biaxial  tensile  and  compressive  strains  on  (a)  the  bulk  band- 
gap  of  QW  material  and  (b)  the  bandgap  lineup  with  two  barrier  layers. 


15 

2.3.2  Bandgap  Lineup 

The  bandgap  difference  between  QW  and  barrier  materials  causes  band  discontinui- 
ties in  both  CB  and  VB  at  the  interfaces  of  the  heterostracture,  as  sketched  in  Fig.  2.4(b). 
The  CB  offset  is  defined  as  A£^  and  the  VB  offsets  are  defined  as  A£**  and  6£u,  for  HH  and 
LH  bands  since  they  are  split  by  the  shear  strain.  Defining  A£g  as  the  bandgap  difference 
between  the  C-HH  strained  bandgap  of  QW  and  the  bandgap  of  barrier  layers,  we  obtain 
A£j  =  £g  (barrier) -£:^(HH),  (2.9) 

where  £j  (barrier)  is  the  bandgap  of  barrier  material.  As  shown  in  Fig.  2.4(a),  CB  edge  of 
the  strained  QW,  as  well  as  VB  edge,  experiences  a  shift  due  to  the  hydrostatic  strain.  In 
reality,  it  is  difficult  to  experimentally  separate  the  CB  shift  from  the  total  bandgap  shift 
[Cor93].  A  common  approach  is,  first,  not  to  worry  about  how  the  total  shift  is  divided  up 
between  the  CB  and  VB,  and  then  for  heterostructures,  the  lineup  of  the  strained  QW  band- 
gap  with  the  barrier  bandgaps  is  typically  determined  by  a  CB  offset  firaction  Q^  which  can 
be  experimentally  measured.  In  this  approach,  the  CB  offset  can  therefore  be  expressed  as 

^E^  =  Q^d^Eg  (2.10) 

and  the  VB  offset  for  HH  band  is  given  by 

A£,,  =   (\-Q^)^E^.  (2.11) 

Since  the  LH  band  is  pushed  away  from  the  HH  band  by  an  amount  S  as  defined  in  Eq. 
(2.8),  the  VB  offset  for  LH  band  can  be  shown  to  be 

A£,^  =  A£^^-S,  (2.12) 

which  means  Af,),  <  Af^j  under  biaxial  compressive  strain  (S  >  0)  and  A£,j  >  tiE^t  under 
biaxial  tensile  strain  (S  <  0),  as  depicted  in  Fig.  2.4(b). 
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2.3.3  Numerical  Results 

InGaAs  OW.  The  material  parameters  as  functions  of  In  and  Al  mole  fractions  for 
InGaAs/AlGaAs  system  have  been  summarized  in  Table  III  of  [Cot93].  As  shown  in  Fig. 
2.5(a),  for  the  unstrained  InoisGao.gsAs,  we  have  £,  (bulk)  s  1 .  194  eV  and  for  the  strained 
Ino.15Gao.s5As  QW/GaAs  barrier,  we  obtain  £,  (barrier)  =  1.424  eV,  Eg  (HH)  s  1.267  eV,  £, 
(LH)  =  1.329  eV  and  S  =  0.062  eV.  The  CB  offset  in  the  strained  InGaAs/GaAs  system  is 
uncertain  at  present.  Various  values  for  Q^  have  been  reported  anywhere  in  the  range  from 
0.4  to  0.8  [Cor93].  We  will  assume  that  Q^  =  0.55  for  the  InGaAs/GaAs  system  in  this  work 
[Wu94,  Hsu97a].  As  a  consequence,  with  AEg  =  0.157  eV,  we  obtain  M,  =  0.087  eV,  AE^, 
=  0.07 1  eV  and  AE,,,  =  0.025  eV,  as  sketched  in  Fig.  2.5(b).  Numerical  values  of  the  mate- 
rial parameters  used  in  this  work  for  the  strained  Ino.15Gao.g5As  QW/GaAs  barrier  system 
are  summarized  in  Table  2.1. 

CdZnSe  OW.  For  II- VI  CdZnSe/ZnSSe  system,  the  ternary  material  parameters  are 
less  well  known.  In  this  work,  we  derive  some  material  parameters  used  in  the  calculation 
for  CdZnSe  and  ZnSSe  ternary  alloys  from  the  material  parameters  of  the  related  binary  al- 
loys based  on  an  interpolation  scheme.  Given  the  present  focus  on  how  strain  modifies  the 
bandgap  of  CdZnSe  QW,  we  leave  the  derivation  of  the  ternary  material  parameters  to  Ap- 
pendix A.  Then  plugging  the  II-VI  ternary  material  parameters  derived  in  Appendix  A  into 
Eqs.  (2.1)  through  (2.8),  we  have  Eg  (bulk)  =  2.373  eV  for  the  unstrained  CdosZno.vSe  and 
for  the  strained  Cdo.3Zno.7Se  QW/ZnSo  0586094  barrier,  we  obtain  Eg  (barrier)  =  2.729  eV, 
5e„  =  0.08 1  eV,  8ej,  u,  =  0.057  eV,  5ej,  „  =  0.04 1  eV,  £,  (HH)  =  2.396  eV,  Eg  (LH)  =  2.494 
eV  and  S  =  0.098  eV,  as  depicted  in  Fig.  2.6(a).  Assuming  Q,  =  0.6  for  the  CdZnSe/ZnSSe 
heterostructure  [Hsu97b]  and  with  Afg  =  0.333  eV,  we  obtain  AE^  =  0.200  eV,  AE^,  =  0. 133 
eV  and  A£,^  =  0.035  eV,  as  sketched  in  Fig.  2.6(b).  Numerical  values  of  the  material  pa- 
rameters used  in  this  work  for  the  strained  Cdo  sZnojSe  QW/ZuSq  06^60.94  barrier  system  are 
summarized  in  Table  2.2. 
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Figure  2.5        The  strain  modified  (a)  bulk  bandgap  and  (b)  potential  wells  of  the  8  nm 
Ino.15Gao.85As/GaAs  QW  under  biaxial  compressive  strain. 
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Parameters  Ino.15Gao.85As  GaAs 


bulk  bandgap  Eg  (eV) 
Luttinger  Parameter  Yi  (x  f?llm^ 
Luttinger  Parameter  Y2  (x  S^/2m„) 
electron  effective  mass  m^  (x  mj 
growth-direction  HH  effective  mass  mjjj.  (x  mj 
growth-direction  LH  effective  mass  rtiu^  (x  m„) 
background  refractive  index  n 
static  dielectric  constant  e^ 


1.194 

IA2A 

8.773 

6.85 

3.041 

2.1 

0.061 

0.067 

0.372 

0.377 

0.067 

0.091 

3.6 

13.33 

Ino.isGaojsAs/GaAs  strained  parameters 


Ino  isGao.gsAs  HH  strained  bandgap  £g(HH)  =  1.267  eV 

Ino.15Gao.g5As  LH  strained  bandgap  £j(HH)  =  1.329  eV 

Ino.15Gao.85As  splitting  energy  5  =  Eg(Ui)  -  Eg(HH)  =  0.062  eV 

Bandgap  difference  A£j  =  0.157  eV  '■   - 

CB  offset  fraction  Q,  =  AEJAE^  =  0.55  .  L  '  '     '"' 

CB  offset  A£^  =  0.087  eV 

HH  band  offset  AE^,  =  0.07 1  eV 

LH  band  offset  A£,^  =  0.025  eV 

Ino.15Gao.85As  in-plane  HH  effective  mass  m/,/^  =  0.085  m„ 

Ino.i5Gao85As  in-plane  LH  effective  mass  Wj^n  =  0.174  m„ 


Table  2. 1         List  of  material  parameters  for  the  Ino,i5Gafl  gsAs/GaAs  strained  QW. 
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Figure  2.6       The  strain  modified  (a)  bulk  bandgap  and  (b)  potential  wells  of  the  4  nm 
Cdo.3Zno.7Se/ZnSo.06Seo.94  QW  under  biaxial  compressive  strain. 
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Parameters  Cdo.3Zno.7Se         ZnSo.o6Sea94 


bulk  bandgap  Eg  (eV) 
Luttinger  Parameter  Yi  (x  S^/2m„) 
Luttinger  Parameter  72  (x  hVlmJ 
electron  effective  mass  ntc  (x  m^) 
growth-direction  HH  effective  mass  m^^j^  (x  m„) 
growth-direction  LH  effective  mass  mn^  (x  m„) 
background  refractive  index  n 
static  dielectric  constant  e. 


2.373 

2.729 

4.495 

4.3 

1.206 

1.14 

0.139 

0.152 

0.480 

0.571 

0.145 

0.157 

2.947 

2.726 

9.524 

9.23 

Cdo.3Zno.7Se/ZnS006Seo.94  strained  parameters 


Cdo.3Zno.7Se  HH  strained  bandgap  £g(HH)  =  2.396  eV 

Cdo.3Zno.7Se  LH  strained  bandgap  £j(HH)  =  2.494  eV 

Cdo.3Zno.7Se  splitting  energy  S  =  £j(LH)  -  £g(HH)  =  0.098  eV 

Bandgap  difference  AEg  =  0.333  eV 

CB  offset  fraction  Q,  =  A£yA£g  =  0.6 

CB  offset  A£,  =  0.200  eV 

HH  band  offset  AEy,  =  0.133  eV 

LH  band  offset  A£,*  =  0.035  eV 

Cdo.3Zno7Se  in-plane  HH  effective  mass  mj^,  =  0. 175  m^ 

Cdo.3Zno.7Se  in-plane  LH  effective  mass  m,;,||  =  0.304  m„ 


Table  2.2         List  of  material  parameters  for  the  Cdo.3Zno.7Se/ZnSo.06Seo.94  strained  QW. 
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2.4  Active  Region  Model 

So  far,  we  have  determined  the  bandgap  lineups  at  interfaces  of  the  strained  hetero- 
structures.  The  next  step,  before  going  ahead  to  intnxiuce  the  gain  and  spontaneous  emis- 
sion functions,  is  to  define  an  active  region  model  for  the  laser  structures  to  be  simulated. 
For  the  QW  laser  structures  studied,  as  shown  earlier  in  Figs.  2.1  and  2.2  for  the  InGaAs 
QW  and  CdZnSe  QW  respectively,  they  both  have  an  active  region  structure  consisting  of 
a  strained  QW  material  sandwiched  between  two  barrier  layers.  It  is  conventional  to  model 
the  actual  active  region  structure  with  the  barrier  layers  on  each  side  of  QW  extending  to 
infinity.  It  is  also  assumed  that  the  potential  wells  associated  with  the  active  region  are 
symmetric  and  rectangular,  like  the  ones  shown  in  Figs.  2.5(a)  and  2.6(a).  As  discussed  in 
[Blo96],  there  is  good  reason  to  believe  that  this  so-called  flat-band  approximation  is  valid 
for  QW  lasers  under  typical  operating  conditions. 

With  the  active  region  model  described  above,  in  the  following  subsections  we  first 
will  solve  the  quantized  electron  and  hole  states  for  the  one-dimensional  potential  wells  in 
CB  and  VB  respectively.  Then  we  will  discuss  the  in-plane  subband  structures  in  a  para- 
bolic band  approximation. 

2.4.1  Potential  Wells  and  Quantized  States 

Fig.  2.7(a)  shows  a  generic  potential  well  diagram  along  the  growdi  direction  of  the 
active  region  model  used  in  the  simulation.  For  the  CB  electrons  and  VB  holes  moving 
along  the  growth  direction,  they  will  experience  potential  wells  with  barrier  heights  equal 
to  the  CB  offset  AE^  and  VB  offset  A£„,  respectively.  As  a  result,  energies  are  quantized 
in  the  potential  wells.  We  define  E^j  and  E^j  as  the  energies  of  theyth  quantized  states  in  the 
CB  and  VB  potential  wells,  respectively,  and  are  measured  relative  to  the  bottoms  of  their 
corresponding  potential  wells.  In  the  VB,  there  are  HH  and  LH  bands  in  our  active  region 
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model.  Therefore,  there  is  the  need  to  clarify  that  the  HH  quantized  energy  Eif^j  and  LH 
quantized  energy  Eujj  are  measured  from  the  bottoms  of  the  HH  and  LH  potential  wells  of 
barrier  heights  A£)^  and  A£,j,  respectively.  Finding  E^  and  £,y  now  becomes  solving  a 
quantum  well  problem. 

Details  of  solving  this  quantum  well  problem  can  be  found  in  many  basic  quantum 
mechanics  textbooks,  here  we  refer  it  to  the  treatment  in  [Cor93].  Basically,  first  it  is  to 
solve  the  Schrodinger's  equation  for  each  section  of  the  active  region  (two  barrier  and  one 
QW).  Then  by  matching  proper  boundary  conditions  across  two  well-barrier  interfaces,  a 
pair  of  characteristic  equations  are  obtained 

HVyJ-^jt;;^  =  °    for  y  is  odd.  an) 

'=°'(vt)  +  ?;r:=0      fow'iseven,  (2.14) 

no    njl. 

with 


(2.15) 


%  =  - ^ -'  •  (2.16) 

where  (f,  is  the  thickness  of  QW,  m„j^  is  the  growth-direction  effective  mass  of  the  QW  ma- 
terial, m„i,  is  the  effective  mass  of  the  barrier  material,  k„jx  is  referred  later  in  Chapter  3  to 
as  the  equivalent  growth-direction  component  of  the  wavevector  in  QW  associated  with 
£„y,  a„j  is  the  magnitude  of  the  complex  wavevector  in  the  barrier  layer  associated  with  E„j 
and  A£,  is  the  barrier  height.  The  n  subscript  indicates  CB  for  electrons  and  HH  and  LH 
bands  in  VB  for  heavy  and  light  holes,  respectively.  Eqs.  (2.13)  and  (2. 14)  can  only  be  sat- 
isfied by  a  set  of  discrete  energies  E„j  which  represent  the  bound  solutions  of  the  corre- 
sponding quantum  well. 


*„A  = 

h 

_  j2m„,(AE„-E^.) 
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Figure  2.7        Generic  diagrams  of  (a)  the  flat-band  potential  wells  and  (b)  the  parabolic 
subband  structure  of  the  QW  active  region  model. 
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2.4.2  Parabolic  In-plane  Subband  Structures 

As  discussed  above,  the  potential  wells  are  one-dimensional  along  the  growth  direc- 
tion. Within  the  plane  of  the  well,  carriers  (electrons  and  holes)  are  not  confined  and  still 
behave  like  "free."  Thus,  for  each  quantized  energy  level  in  the  QW,  an  in-plane  energy 
subband  exists,  as  shown  in  Fig.  2.7(b).  The  subband  structure  is  usually  characterized  by 
defining  an  in-plane  effective  mass  as  [Kit86] 

d^E  „"!-' 

njk 


(2.17) 


where  m„^i  is  the  in-plane  effective  mass,  E^^^  is  the  in-plane  energy  measured  relative  to 
the  bottom  of  the  well  and  k,.  is  the  in-plane  wavevector  associated  with  E^^  .  They  sub- 
script indicates  they'th  subband  in  the  well  and  the  n  subscript  indicates  CB  for  electrons, 
HH  and  LH  bands  for  heavy  and  light  holes,  respectively.  If  the  parameter  m„f\  is  a  constant 
with  respect  to  E-^.  ,  the  subband  structure  is  said  to  be  parabolic.  In  a  nondegenerate 
band,  such  as  the  CB  (excluding  spin  degeneracy),  all  the  subbands  are  assumed  to  be  par- 
abolic and  the  electron  effective  masses  along  the  1  and  II  directions  are  identical  (m^x  = 
nicf^.  However,  the  in-plane  subband  structure  in  a  degenerate  band,  such  as  the  VB,  is 
more  complicated,  since  band  coupling  between  the  HH  and  LH  bands  can  be  very  strong. 
In  that  case,  the  VB  subbands  become  nonparabolic  band  due  to  this  so-called  valence  band 
mixing  effect.  Since  the  quantum  wells  are  under  biaxial  compressive  strain  for  the  laser 
structures  studied  in  this  work,  the  LH  band  is  pushed  away  from  the  HH  bands  due  to  the 
shear  component  of  the  strain.  This  greatly  reduces  the  valence  band  mixing  effect  between 
the  HH  and  LH  bands  and  allows  us  to  assume  that  the  subband  structures  in  the  VB  are 
also  parabolic.  In  any  event,  the  growth-direction  HH  and  LH  effective  masses  are  not  af- 
fected by  the  strain  and  are  given  by  [Cor93] 

1      h'^  1      h?-  ,.  ,„, 

m,,L  I   =  • r — —  ,       m,.  I  = — — ,  (2.16) 

hhL      Yj_2y2  2  ""-L     7i-i-2Yj2 
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where  the  material  constants  Yi  and  Y2  are  referred  to  as  the  Luttinger  parameters  [Lut56]. 
In  the  high  strain  approximation,  the  in-plane  HH  and  LH  effective  masses  can  be  ex- 
pressed in  close  foims  with  in  terms  of  Luttinger  parameters  y,  and  Y2  (in  the  units  of 
h'^/2mj  [Cor93] 

1      h^  l_ft2 

'"**ii  =  :^^j  ■     '"Mil  -  :^-::^j  ■  (2.19) 

From  the  preceding  discussion  of  strain  and  by  comparing  Eqs.  (2.18)  and  (2.19),  we  find 
that  not  only  do  the  HH  and  LH  bands  split  apart,  but  also  the  HH  band  becomes  "lighter" 
and  the  LH  band  becomes  "heavier"  in  the  plane  of  compression.  Thus,  this  is  contrary  to 
the  bulk  material  where  the  effective  masses  along  1  and  II  directions  are  identical. 

2.4.3  Numerical  Results 

InGaAs  OW.  For  the  IngisGao  ssAs  QW,  we  obtain  electron  effective  mass  m^  =  m^^ 
=  m^ii  =  0.061m„,  Luttinger  parameters  Yi  =  8.773,  Y:  =  3.041  and  for  the  GaAs  barrier,  we 
have  m^  =  rtici,  =  0.067m„  Yi  =  6.85  and  Y2  =  2.1  using  the  material  parameters  for  the  In- 
GaAs/AlGaAs  system  summarized  in  Table  III  of  [Cor93].  Plugging  the  numerical  values 
for  Yi  and  Y2  in  Eqs.  (2.18)  and  (2.19),  we  obtain  m^^x  =  0.372m„,  m,^  =  0.067m,,,  m«,||  = 
0.085m„  and  m;j||  =  0.174m„  for  the  Ino.15Gao.g5As  QW  and  mi,i^  =  0.377m„  and  m,!^  - 
0.09  lm„  for  the  GaAs  barrier.  Having  all  the  values  of  effective  masses  plugged  into  Eqs. 
(2.13)  through  (2.16),  we  solve  the  quantized  levels  for  the  CB  and  VB  (HH  and  LH)  po- 
tential wells  of  width  d,  =  S  nm,  as  shown  in  Fig.  2.5(b).  The  calculation  yields  one  electron 
subband  in  the  CB  (£cj  =  31-6  meV  relative  to  the  bottom  of  CB  potential  well),  two  HH 
subbands  (Eum  =  9.2  meV  and  Ei,[j2  =  35.4  meV  relative  to  the  bottom  of  HH  potential 
well)  and  one  LH  subband  (E/^i  =  6.6  meV  relative  to  the  bottom  of  LH  potential  well)  in 
the  VB.  Numerical  values  of  the  material  parameters  used  in  this  work  for  the  strained 
Inoi5Gaog5As  QW/GaAs  barrier  system  are  summarized  in  Table  2.1. 
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CdZnSe  OW.  For  the  CdojZno jSe  QW,  we  obtain  electron  effective  mass  m^  =  m^^^ 
^m^fl  =  0.139/n„,  Luttinger  parameters  7i  =  4.495,  Yj  =  1-206  and  for  the  ZnSo.06Seo.94  bar- 
rier, we  have  m^  =  m^t  =  0.152m<„  Yi  =  4.3  and  Y2  =  114  using  the  ternary  material  param- 
eters for  CdZnSe/ZnSSe  system  derived  in  Appendix  A.  Using  the  numerical  values  for  Yi 
and  Y2  in  Eqs.  (2.18)  and  (2.19),  we  obtain  m/jj.  =  0.480m„,  nii^  =  0. 145m„,  nii^  =  0.175m<, 
and  m,*||  =  0.304m„  for  the  Cdo.3Zno.7Se  QW  and  nt/j^  =  0.571m„  and  m,i^  =  0.157m„  for 
the  ZnSo.06Seo.94  barrier.  Plugging  all  the  values  of  effective  masses  into  Eqs.  (2.13) 
through  (2.16),  we  solve  the  quantized  levels  for  the  CB,  HH  and  LH  potential  wells  of 
width  d,  =  4  nm,  as  shown  in  Fig.  2.6(b).  The  calculation  yields  one  electron  subband  in 
the  CB  (Ec,  =  62.6  meV),  two  HH  subbands  (£hhi  =  23.9  meV,  Ehhi  =  90.3  meV)  and  one 
LH  subband  (Eu/^  =  23.1  meV)  in  the  VB.  Numerical  values  of  the  material  parameters 
used  in  this  work  for  the  strained  (Mo.3Zno.7Se  QW/ZnSo.o6Seo.94  barrier  system  are  summa- 
rized in  Table  2.2. 

Comparison.  As  shown  above,  electron  and  hole  effective  masses  of  the  CdZnSe/ 
ZnSSe  system  are  much  larger  than  those  of  the  InGaAs/GaAs  system.  In  Section  2.6.3, 
the  difference  in  the  effective  masses  will  be  shown  to  play  an  important  role  in  the  differ- 
ence of  the  transparency  carrier  densities  between  the  two  QW  diode  laser  systems. 

2.5  Single-Panicle  Gain  and  Spontaneous  Emission  Functions 

2.5.1  Transition  Energy 

In  Section  2.4,  we  have  shown  how  the  strain  modifies  the  bandgap  of  a  QW  and  how 
the  quantum  effect  modifies  its  strained  bandgap.  As  illustrated  in  Fig.  2.7(a),  we  define 
E^  as  the  "modified"  energy  gap  between  the/th  conduction  and  valence  subbands 

^.y  =  ^.  +  ^cy  +  £v;.  (2.20) 

where  E.  is  E  (HH)  for  the  C-HH  strained  bandgap  or  £.(LH)  for  the  C-LH  strained  band- 
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gap.  Eg  and  E,j  are  the  quantized  energy  levels  in  the  QW.  For  an  electron-hole-pair  par- 
ticipating in  a  transition  with  energy  £,.„  as  depicted  in  Fig.  2.7(b),  we  obtain 

£cv    =    f,  +  ^cA  +  £v;.,..  (2-21) 

where  E  ^  and  £  ,,  are  the  energies  of  electrons  and  holes  respectively  participating  in 
this  transition  and  are  measured  relative  to  the  bottoms  of  CB  and  VB  potential  wells.  In 
terms  of  the  in-plane  wavevector  kn,  the  expressions  for  £  -^  and  E-^^  arc  given  by 

where  m^jn  and  m,yii  are  the  in-plane  effective  masses  and  both  conduction  and  valence  sub- 
bands  are  assumed  to  be  parabolic  as  discussed  in  Section  2.4.2. 

2..5.2  Gain  and  Spontaneous  Emission 

In  simulating  the  operating  characteristics  of  standard  QW  diode  lasers,  it  is  custom- 
ary to  utilize  two  functions  which  describe  the  dependence  of  QW  gain  and  spontaneous 
emission  on  photon  energy  tuo.  In  deriving  these  functions,  defining  g  (Sco)  as  the  spectral 
gain  function  and  R^p  (Sw)  as  the  spontaneous  emission  function,  approximations  with  var- 
ious levels  of  accuracy  can  be  made.  In  the  first  level  approximation,  one  ignores  many- 
body  effects  and  assumes  that  electrons  trapped  in  the  QW  interact  only  with  the  electro- 
magnetic field.  In  this  "single-particle"  model,  where  fto)  is  equal  to  the  transition  energy 
£„,  we  have  g  (ha)  =  g  (£„)  and  R,p  (Sco)  =  «^^  (£„).  g  (£„)  and  R^p  (£„)  are  calculated  by 
summing  over  tiie  contributions  from  all  the  allowed  subband  transition  pairs  as  given  by 

«(£cv)   =Z«;(^cv).  (2-23) 

;■ 

«.p(£cv)  =!«,,;(£  J.  (2-24) 

J 

where  gj  (£„)  and  R^j  (£„)  are  the  spectral  gain  function  and  spontaneous  emission  function 
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for  an  allowed  subband  transition  pair,  respectively.  Since  there  is  only  one  subband  in  the 
CB  in  our  structure,  there  are  two  allowed  transition  pairs,  viz.  Cl-HHl  and  Cl-LHl.  In 
the  single-particle  model,  one  ignores  many-body  effects  and  assumes  that  momentum  is 
conserved  (k-selection  holds)  in  the  photon  emission  and  absorption  processes.  In  this 
case,  it  can  be  shown  that  gj  (£„)  and  R^pj  (£„)  are  given  by  [Cor93] 

«/^"^  =  ii.)fBk\''r\'PreAfc(^cj.)  -A(£vW  -1]  .  (2.25) 


R,pj(EJ    =    (»»);^^2H5;72N'p...[/c(£c;K,)//.(V.)] 


(2.26) 


where  e  is  the  electron  charge,  E„  is  the  vacuum  permittivity,  c  is  the  vacuum  light  speed, 
m„  is  the  electron  rest  mass,  n  is  the  background  refractive  index  of  the  QW  material, 
|A/j.P  is  the  transition  matrix  element,  \M-j\^  is  the  average  transition  matrix  element,  pr^j 
is  the  reduced  density  of  states,/^  is  the  Fermi-Dirac  distribution  for  CB  electrons  and/*  = 
1  -/v  is  the  distribution  for  VB  holes.  The  transition  matrix  element  \Mj\^  is  given  by 

^  +  i  cos^e^       for  C-HH  TE  gain 


IMP 


2     2 
l-cos^e^ 

--^cos^e^. 
i  +  cos^e. 


for  C-HH  TM  gain 
for  C-LH  TE  gain 

for  C-LH  TM  gain 


(2.27) 


where 


IMP 


]    (E+A) 


m„E_ 


-^^A')' 


(2.28) 


C0S^6  ; 


(^.■^^vy)^ 


rjW 


(£cv-£,-£c;-£v;)  +  (^c;  +  V;„ 


(2.29) 


For  the  In,Gai,jAs  QW,  |A/p  can  be  simplified  and  expressed  as  a  function  of  indium  mole 


fraction  [Cor93] 
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IMP  =  I  (28.8  -  6.6x)  (em J  .  (2.30) 


The  parameters  m^jj.  and  m^fl  are  the  reduced  effective  masses  along  the  growth  and  in-plane 
directions,  respectively  and  are  defined  as  follows: 

-L  =  -L  +  _L,    J L  +  _L.  (2.31) 

"'rjl       '"cjl      '"v;X  '"rjW         "'cjW      '"vyil 

The  average  transition  matrix  element  |A/j-p  can  be  shown  to  be 

IMtP      2 

'      '    =  -      for  both  C-HH  and  C-LH  spontaneous  emission .  (2.32) 

Derivations  of  Eqs.  (2.27)  and  (2.32),  as  well  as  Eq.  (2.29)  which  shows  that  9y  is  a  function 
of  transition  energy,  are  given  in  Appendix  B.  The  reduced  density  of  states  p,,^  (spin  not 
included)  is  given  by  [Cor93] 


(2.33) 


^^ '^^^J^}  -  1  +  exp  [  (£^.,^^ - E^^)  /k,T]  '  '•^■^'^^ 

^"^^v;".,,)   =  l  +  exp[(£„.,^^-£p/*,r]'  _    ^^-3^) 

where  E^^  and  Ef,  are  the  quasi-Fermi  levels  in  the  CB  and  VB  respectively  and  are  mea- 
sured from  the  CB  and  VB  edges  into  the  bands,  kg  is  the  Boltzmann  constant  and  Tis  the 
temperature  of  the  structure.  For  parabolic  conduction  and  valence  subbands,  we  can  write 
the  electron  concentration  N  and  hole  concentration  P  in  the  quantum  well  as 

^  =  r^I'"c;«ln  { 1  +  exp  [- (£,y  -  £^,)  /k,T]  }  (2.36) 

Tin  a^  j 
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P  =  ;^I'"v;i|ln  { 1  +  «P  [- (£v; -  £/v)  /*fl^]  }  (2.37) 

The  summation  is  over  all  subbands  which  includes  HH  and  LH  subbands  for  P  in  the  VB. 
Under  the  condition  of  overall  charge  neutrality  in  the  quantum  well,  we  have  N  =  P.  For 
a  given  value  of  N,  we  can  solve  for  the  values  of  Ef^  and  Ej,  which  satisfy  Eqs.  (2.36)  and 
(2.37)  at  a  given  T. 

2.6  Simulation  Results 

The  general  procedure  for  calculating  the  single-particle  g  (.E„)  and  R^p  (£„)  func- 
tions is  summarized  as  follows:  (1)  Specify  A'  and  T;  (2)  Input  numerical  values  of  material 
parameters  discussed  in  Section  2.4.3;  (3)  Input  the  quantized  energy  levels  solved  in  Sec- 
tion 2.4.3;  (4)  Solve  E^.  and  £^  from  a  plot  ofEfi.  and  Ef,  as  a  function  ofN  using  Eqs.  (2.36) 
and  (2.37);  (5)  Calculate  g  (£„)  and  R^  (£„)  using  Eqs.  (2.23)  through  (2.26). 

2.6.1  InGaAsOW 

Fig.  2.8  shows  plots  of  Ef^  and  Ef^vs.N  for  the  8  nm  Ino.isGao  ^As  strained  QW  with 
GaAs  barriers.  Because  the  QW  is  under  biaxial  compressive  strain,  the  first  VB  subband 
to  be  filled  up  with  holes  is  a  HH  subband.  From  Eq.  (2.27),  it  indicates  that  the  TE  tran- 
sition matrix  element  is  stronger  than  its  TM  counterpart  for  a  C-HH  transition.  Asa  result, 
TE  gain  dominates  over  TM  gain  for  the  InGaAs  strained  QW  lasers.  Typical  TE  gain 
spectra  for  A^  =  1 .3  x  10"  and  1.7  x  10"  cm''  at  r  =  300  K,  calculated  from  Eqs.  (2.23)  and 
(2.25),  are  shown  in  Fig.  2.9(a).  It  shows  that  gain  increases  with  increasing  N.  The  spec- 
tral position  of  the  peak  gain  occurs  at  E^j  u,<  which  corresponds  to  the  energy  gap  between 
CI  and  HHl  subbands  defined  in  Eq.  (2.20),  and  remains  there  as  N  increases.  Since  no 
transition  is  allowed  for  £„  smaller  than  E^,^  ^,  the  gain  for  photons  with  energies  below 
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Eg,  ^  is  zero.  Fig.  2.9(b)  shows  the  spontaneous  emission  spectra  for  the  corresponding  N, 
calculated  from  Eqs.  (2.24)  and  (2.26).  Two  unrealistic  discontinuous  peaks  occurring  at 
£j/  u,  and  £,;,  //,  in  the  spontaneous  emission  spectra  are  due  to  the  step-like  density  of  states 
in  a  QW  structui*  [Coi93].  For  the  gain  spectra  shown  in  Fig.  2.9(a),  only  the  discontinuity 
at  £,;,  ij,  is  obvious.  The  one  at  Eg,^  u,  should  also  show  up,  but  it  is  in  the  absorption  region 
which  is  beyond  the  scale  of  the  plot. 

■.'■■■  v\ 

2.6.2  CdZnSe  OW  |  pf  f.  , ;-.  ■■  \  ;■  n  f 

Fig.  2. 10  shows  plots  of  E^^  and  Efi,\s.N  for  the  4  nm  Cdo  sZnojSe  strained  QW  with 
ZnSo.o6Seo.94  barriers.  Typical  TE  gain  spectra  for N  =  4.6  x  10"  and  5.4  x  10'*  cm'  at  T 
=  300  K,  calculated  from  Eqs.  (2.23)  and  (2.25),  are  shown  in  Fig.  2.1 1(a).  SimUar  to  the 
gain  spectra  in  Fig.  2.9(a),  only  the  discontinuity  at  £,;,  w,  is  obvious.  In  Chapter  5,  when 
we  discuss  the  absorption  spectra  in  the  CdZnSe  QW  lasers,  we  will  show  clearly  that  both 
discontinuities  exist  within  single-particle  model.  Fig.  2.11(b)  shows  the  spontaneous 
emission  spectra  for  the  corresponding  N,  calculated  from  Eqs.  (2.24)  and  (2.26).  They  also 
show  unrealistic  spectral  features  of  discontinuity  due  to  the  step-like  density  of  states  for 
C-HHl  andC-LHl  transitions. 

2.6.3  Comparison 

From  Eq.  (2.25),  we  see  that  when 

/c(£W+/.(£vyK.)-l>0.  (2.38) 

gj(,E„)  is  positive  and  an  incoming  light  wave  with  photon  energy  £„  will  be  amplified  by 
the  material.  Using  Eqs.  (2.34)  and  (2.35)  in  Eq.  (2.38),  one  can  show  that  this  inequality 
is  equivalent  to  saying  that 
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Figure  2.8        Quasi-Fermi  level  as  a  function  of  carrier  density  for  (a)  CB  electrons  and 
(b)  VB  holes  of  the  8  nm  Ino.1sGao.g5As/GaAs  QW. 
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Figure  2.9  Single-particle  (a)  TE  gain  and  (b)  spontaneous  emission  spectra  for  the  8 
nm  Ino.15Gao.g5As/GaAs  QW  at  various  indicated  carrier  densities  (in  units 
of  10"  cm')  at  room  temperature. 
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Figure  2.10      Quasi-Fermi  level  as  a  function  of  carrier  density  for  (a)  CB  electrons  and 
(b)  VB  holes  of  the  4  nra  Cdo.3Zno.7Se/ZnSo.06Seo.94  QW. 
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Figure  2.10  Single-particle  (a)  TE  gain  and  (b)  spontaneous  emission  spectra  for  the  4 
nm  Cdo.3Zno.7Se/ZnSo.06Seo.94  QW  at  various  indicated  carrier  densities  (in 
units  of  10"  cm'')  at  room  temperature. 
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£^<^cv<£w  +  ^/c  +  ^/v  =  ^^/  (2.39) 


is  the  requirement  for  gain  at  a  photon  energy  E„.  In  other  words,  the  quasi-Fermi  level 
separation  A£f  must  be  greater  than  the  QW  subband  energy  gap  E^  (or  A£y-  E^j  >  0)  to 
achieve  optical  gain  in  the  material.  Eq.  (2.39)  is  similar  to  the  one  for  bulk  semiconductors 
given  in  [Cas78,  Cor93]  where  E^j  is  replaced  by  bulk  Eg.  For  a  QW  under  biaxial  com- 
pressive strain,  E^j  indicates  the  smallest  Cl-HHl  subband  energy  gap  £j;,  u,-  P'o's  of  (A£^ 
-  Eg,^  ^vs.Naie  shown  in  Fig.  2. 12  for  both  the  InGaAs  QW  and  CdZnSe  QW  structures 
studied.  The  cross-over  of  the  (A£y-  £,,«,)  vs.  N  curve  on  the  W-axis  is  defined  as  the  trans- 
parency carrier  density  N„  which  is  the  electron-hole  plasma  density  required  to  provide  the 
quasi-Fermi  level  separation  so  that  the  material  will  become  transparent  for  photon  energy 
equal  to  Egj  according  to  Eq.  (2.39).  Optical  gain  in  the  material  is  attained  when  we  inject 
a  carrier  density  N  >  N,,  such  that  A£}-  >  £j;  j^.  From  the  curves,  we  find  that  N„  =  0.9  x 
10"  cm-^  for  the  strained  InGaAs  QW  which  is  smaller  than  N,^  =  3.1  x  10'*  cm"'  for  the 
strained  CdZnSe  QW.  This  is  mainly  due  to  that  the  density  of  states  for  the  CdZnSe  QW 
is  larger  than  that  for  the  InGaAs  QW.  From  Eq.  (2.33),  the  density  of  states  is  found  to  be 
proportional  to  the  carrier  effective  masses  of  the  material  system  and,  as  discussed  in  Sec- 
tion 2.4.3,  the  effective  masses  of  the  CdZnSe  QW  are  much  larger  than  those  of  the  In- 
GaAs QW.  Another  observation  is  that  the  slope  of  the  curve  for  the  InGaAs  QW  is  larger 
than  that  for  the  CdZnSe  QW.  This  means  that  carrier  fills  up  the  subbands  in  the  InGaAs 
QW  faster  than  in  the  CdZnSe  QW.  Effect  of  density  of  states  on  the  transparency  condi- 
tion is  discussed  in  great  detail  in  [Cor93].  In  Chapter  4,  we  will  see  how  the  difference  in 
A',,  contributes  to  the  difference  in  transparency  current  densities  for  the  two  material  sys- 
tems. 

2.7  Summary 

Optical  gain  and  spontaneous  emission  spectra  are  calculated  for  both  InGaAs  and 
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CdZnSe  QW  laser  structures  using  the  single-particle  model.  Since  the  measured  spectra 
do  not  have  sharp  features,  it  is  clear  that  this  single-particle  model  is  incomplete.  In  Chap- 
ter 3,  we  show  that  when  many-body  effects  are  incorporated  into  the  single-particle  model, 
realistic  gain  and  spontaneous  spectra  are  obtained.  However,  it  is  interesting  to  note  that 
in  the  limit  of  very  thick  QWs,  the  single-particle  model  gives  realistic-looking  spectra. 
This  is  due  to  the  fact  that  the  stair-step  density  of  states  function  gradually  smooths  out  as 
QW  thickness  increases  [Cor93].  As  a  consequence,  the  single-particle  model  applied  to 
thick  active  layer  lasers  generates  spectra  which  look  realistic  and  it  is  not  obvious  that  the 
model  is  incomplete. 
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CHAPTERS 
MANY-BODY  EFFECTS 


3.1  Introduction 


In  Chapter  2,  we  have  discussed  the  optical  gain  and  spontaneous  emission  functions 
based  on  a  single-particle  model  where  mutual  interactions  among  carriers  are  ignored. 
However,  upon  current  injection  into  the  semiconductor  QW  diode  lasers,  the  carrier  den- 
sities in  the  QW  are  typically  on  the  order  of  10"  ~  10"  cm'.  In  the  presence  of  such  high 
density  of  electron-hole  plasma,  mutual  interactions  among  the  carriers  are  expected  to  be 
important  These  mutual  interactions  are  usually  referred  as  "many-body  effects"  and  in- 
clude carrier  scattering  (CS),  bandgap  renormalization  (BGR)  and  Coulomb  enhancement 
(CE).  In  this  chapter,  effects  of  CS,  BGR  and  CE  on  the  single-particle  g  (E^y)  and  R,p  (£„) 
functions  will  be  discussed  sequentially. 

-3.2  Carrier  Scattering  (CS^ 

3.2.1  Gain  and  Spontaneous  Emission  Functions 

In  high  quality  QW  laser  material,  CS  is  the  main  contributor  to  spectral  broadening 
and  is  taken  into  account  by  convolving  the  single-particle  gj  (£„)  and  R,^  (£„)  functions 
for  a  given  subband  transition  pair  with  a  spectral-broadening  lineshape  function.  The  con- 
volved gain  function  can  be  expressed  as  [Asa93] 
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gg{n(a)  =lteat['£f  gjiEJ 


xLj(.hu>-EJ 


-ha 


l-h 


and  the  convolved  spontaneous  emission  is  given  by 
xLj(ha>-EJ 


dE  , 


(3.1) 


£„.  -  fiw 


1-1= 


dE^, 


(3.2) 


where  the  notation  Real  refers  to  the  real  pan  of  a  complex  quantity.  The  subscript  B  indi- 
cates the  spectral  broadening  due  to  carrier  scattering.  Since  gj  (£„),  R^pj  (E„)  and,  as 
shown  later,  Lj{ha  -  E„)  are  all  real  quantities,  Eqs.  (3.1)  and  (3.2)  are  often  rewritten  as 

«b(»»)     =lj^«/£cv)^;(^"-£cv)^£.v  (3-3) 


V('*")    =  %Z'^sp^^c.)h^^^-EJdE, 


(3.4) 


by  dropping  the  square  brackets  which  are  the  only  complex  terms  in  these  equations  and 
thus  dropping  the  notations  Real.  Eqs.  (3.3)  and  (3.4)  mean  that,  for  a  given  subband  tran- 
sition pair,  the  transitions  with  energy  £„  /  Sco  contribute  to  optical  gain  and  spontaneous 
emission  at  tim  with  the  weight  Lj  (ha  -  E„)  and  the  total  optical  gain  and  spontaneous  emis- 
sion are  the  summation  over  the  contributions  from  all  the  allowed  subband  transition  pairs. 
The  width  of  the  spectral-broadening  lineshape  function  is  closely  related  to  the  time 
constants  associated  with  various  carrier  scattering  processes  [Asa93,  Yam87].  The  param- 
eters Xcj  and  Tv,  denote  the  electron  scattering  time  and  hole  scattering  time  in  theyth  sub- 
bands  of  CB  and  VB  respectively.  Various  CS  processes  contributing  to  x^  are  electron- 
electron  (e-e),  electron-hole  (e-h)  and  electron-LO  phonon  (e-LO)  and  those  contributing 
to  Tv;  are  hole-hole  (h-h),  hole-electron  (h-e)  and  hole-LO  phonon  (h-LO).  Details  concern- 
ing the  calculation  of  the  these  scattering  times  can  be  found  in  [Asa93]. 
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The  non-Lorentzian  fomi  used  for  the  lineshape  function  is  given  by  [Asa93,  Hsu97a, 
Yam87] 

r    (h(a,E    ) 


I     ffk  F     \  cv  ^        '      cv 

Lj  (RO)  -  is„)  =  -  ^^^  _  2  +  rl  (/iCO,  EJ  ' 


(3.5) 


where 


r„  (««a.  £cv)    =   r,,k,  Ci"  -  £.v  +  ^cyk,)  +  ^.jV.,  ('*'«  -  ^cv  +  £v;k„)  •  (3.6) 

This  lineshape  function  differs  from  a  Lorentzian  function  in  that  r„  depends  on  photon 
energy  fico.  T-y^  for  the/th  VB  subband  can  be  approximated  as  [Asa93] 


^vK^^^-2i 


n  ,,e''-'-^'-'-"" 


l+e-^A''*'" 


where 


with 


<^\i         T(0,k^J^,k^J^)r2m^J„kT^2 


4SKEle}hkTdl     1  +  g-^A^*^ 


no.v.V)  =  ^il'x[7^,-^^^'-^'''^-W^ 


(3.7) 


(3.8) 


1         ^?  +  8*v', 


v;l 


^}    C^^^^.j^ 


(3.9) 


4(^^  +  «?;i)^[exp(rf„^)-l]. 


The  parameter  Xj  is  the  inverse  screening  length  and  is  given  by  [Asa93] 

^2 


X?  = 


OS) 


(3.10) 


c;  vj 

where  e,  is  the  static  dielectric  constant  of  the  quantum  well  material,  d^  is  the  effective 
well  width  for  theyth  subband  in  the  CB  and  VB  (n  =c  and  v,  respectively).  The  summa- 
tion over_/  means  that  the  screening  effect  is  considered  with  inclusion  of  all  the  subbands 
in  CB  and  VB.  d„j  is  defined  by  [Hsu97a] 
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d  .  =  j-^M=  =  j^,  (3.11) 

where  the  parameter  it  j^  is  the  equivalent  growth-direction  component  of  the  wavevector 
associated  with  the  quantized  energy  level  E„j.  For  they'th  CB  subband,  V^j^  can  be  cal- 
culated by  replacing  all  the  subscripts  v  with  c  in  Eqs.  (3.7)  through  (3.9).  Details  concern- 
ing the  derivation  of  Eq.  (3.7)  and  the  associated  parameters  are  given  in  [Asa93]. 

3.2.2  Simulation  Results  ; .    v  ;  '  -.      ■• 

InGaAs  OW.  The  convolved  TE  gain  and  spontaneous  emission  spectra  with  inclu- 
sion of  spectral  broadening  due  to  carrier  scattering  are  plotted  as  the  solid  curves  in  Figs. 
3.1(a)  and  (b),  respectively,  at  W  =  1.7  x  10'*  cm'  and  7  =  3(X)  K.  It  is  evident  that  these 
spectra  are  broadened  and  smoothed  relative  to  their  single-particle  counterparts  which  are 
replotted  as  the  dotted  curves  in  Figs.  3. 1  (a)  and  (b)  for  easy  comparison.  The  peak  gain  is 
reduced  significantly  and  its  spectral  location  is  shifted  away  from  the  Cl-HHl  subband 
energy  gap  £,;,  u,  to  higher  photon  energy  (in  this  case  by  about  1 1  meV).  As  discussed 
later  in  Chapter  4,  this  shift  is  important  in  determining  the  dependence  of  lasing  wave- 
length on  cavity  length. 

CdZnSe  OW.  The  solid  curves  in  Figs.  3.2(a)  and  (b)  are  the  convolved  TE  gain  and 
spontaneous  emission  spectra,  respectively,  of  their  single-particle  counterparts  plotted  as 
the  dotted  curves  for  the  4  nm  Cdo.3Zno,7Se/ZnSo.o6Seo.94  QW  at  A'  =  5.4  x  10"  cm''  and  T 
=300  K.  In  this  case,  peak  gain  is  also  reduced  and  its  spectral  location  is  blue-shifted  from 
Egi,  uii  to  higher  photon  energy  by  about  1 1.5  meV. 
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Figure  3.1  The  convolved  (a)  TE  gain  and  (b)  spontaneous  emission  spectra  (solid 
curves)  of  their  single-panicle  counterparts  (dotted  curves)  for  the  8  nm 
Ino.15Gao.85As/GaAs  QW  at  N  =  1.7  X  10"  cm-'  and  T  =  300  K.  The  convo- 
lution takes  into  account  the  spectral  broadening  due  to  carrier  scattering. 
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Figure  3.2  The  convolved  (a)  TE  gain  and  (b)  spontaneous  emission  spectra  (solid 
curves)  of  their  single-particle  counterparts  (dotted  curves)  for  the  4  nm 
Cdo.3Zno.7Se/ZnSo.o6Seo,94  QW  at  W  =  5.4  x  10'*  cm''  and  7  =  300  K.  The 
convolution  takes  into  account  the  spectral  broadening  due  to  carrier  scat- 
tering. 
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3.3  Bandgap  Renormalization  (BGR) 

3.3.1  Plasma  Screening 

One  of  the  most  important  many-body  interactions  in  an  electron-hole  plasma  (EHP) 
is  the  screening  of  the  Coulomb  potential.  In  a  two-dimensional  (2D)  quantum  well  sys- 
tem, the  Fourier  transformation  of  the  unscreened  Coulomb  potential  in  the  real  space  to 
the  momentum  space  (or  wavevector  k  space)  is  given  by  [Cho94,  Hau90] 

V„   =  -— ^ (3.12) 

"•      2e„e,<?„A 

where  A  is  the  in-plane  area  of  the  quantum  well.  The  parameter  qu  is  the  magnitude  of  a 
qil  vector  which  is  defined  as  the  difference  between  in-plane  wavevector  k,.  and  k,, ,  that 
is  q,.  =  k,'|  -  k|, .  Given  the  angle  6  between  kl  and  k,, ,  qll  can  be  given  from  the  law  of 
cosine 

q^  =  it|'|2  +  kj  -  2A|'|jt||Cose .  (3.13) 

Defining  the  exciton  Bohr  radius  fl„  and  exciton  Rydberg  constant  Eg  as  (in  M.K.S.  units) 
[Cho94,  Hau90] 


4jcS^e  E, 

o  s 


(3.14) 


Er  =  7"^.  (3-15) 


we  obtain 


e^ 


87ce„e, 


(3.16) 


where  a^g  depends  only  on  the  static  dielectric  constant  e,  of  the  material.  Unlike  in 
[Cho94,  Hau90],  we  use  the  in-plane  reduced  effective  mass  discussed  in  Chapter  2  to  take 
the  effective-mass  anisotropy  into  account  in  Eqs.  (3.14)  and  (3.15).  Using  (3.16)  in  Eq. 


46 


(3.12),  the  unscreened  Coulomb  potential  can  be  expressed  in  terras  of  a„  and  E^  as 
The  screened  Coulomb  potential  is  given  by  [Cho94,  Hau90] 


(3.17) 


'I'      e|,(?||,co) 


(3.18) 


where  e,,  (i?,,,  co)  is  the  longitudinal  dielectric  function  and  is  frequency  dependent  (dy- 
namic). At  high  temperature  where  screening  is  governed  by  free  carriers  (not  excitons), 
screening  effect  can  be  described  within  the  so-called  random  phase  approximation  (RPA) 
[Cho94,  Hau90].  Replacing  a  continuum  of  poles  in  RPA  by  a  single  effective  plasmon- 
pole  and  neglecting  the  frequency  dependence  of  screening  effect,  one  obtains  the  much 
simpler  static  plasmon-pole  approximation  (SPPA)  [Cho94,  Hau90].  In  this  case,  the  lon- 
gitudinal dielectric  function  can  be  expressed  as 
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eiiC^ii-O) 


1- 


1' 

„2   ' 


(3.19) 


where  O)^/  is  the  plasma  frequency  and  co    is  the  effective  plasmon  frequency.  For  a  gen- 
eral multi-subband  QW  system,  cOp;  can  be  expressed  as  [Cam94] 


J  V 


{'"cjW      ""vjllj 

V 

l^.'A 

;■ 

['"cjn    '"v 

fl) 

(3.20) 


where  the  summation  is  overall  all  the  conduction  and  valence  subbands.  In  the  last  form 
of  (i>l, ,  the  expression  for  V    in  Eq.  (3.17)  is  used.  The  effective  plasmon  frequency  co 
is  related  to  co^,  by  [Cho94,  Hau90] 


1  + 


<?!! 


^jl  + 


(3.21) 


where  the  inverse  screening  length  X,  is  given  in  Eq.  (3.10)  with  inclusion  of  all  the  sub- 
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bands.  Dividing  both  sides  of  Eq.  (3.21)  by  wj^  and  taking  the  inverse,  we  obtain 


^ 


1  + 


^li 


"pi 


(3.22) 


and  therefore, 


1-- 


•!  S 


2-1 

2n-l 


(3.23) 


Before  further  expanding  Eq.  (3.23),  we  introduce  several  new  parameters  which,  in  anal- 
ogy to  a„  and  E^  defined  Eqs.  (3.14)  and  (3.15),  are  defined  as 


a„  =  — 7. ,     a„„  =  — ^ , 


"cjll 


E ^!_ 

*^"2m,,„a?' 


"v;ll 


E,..= 


h^ 


2'"vjll«Jv  ' 


(3.24) 


(3.25) 


cji\    oc 

These  parameters  a„c>  i<n>  ^rc  and  ^kv  don't  have  any  physical  meaning  except  for  the  con- 
venience of  mathematical  manipulation.  Similar  to  Eq.  (3.16),  we  obtain 

p2 


O    D 


=  '^o^R 


Using  Eq.  (3.26)  in  Eq.  (3.24)  to  solve  for  the  effective  masses,  we  have 


J_ 

"c;ll 


::^  =  ^  ^"o^R^ " 


1       2  ,    _  , 


(3.26) 


(3.27) 


For  the  terms  containing  cOp,  in  both  square  braclcets  of  Eq.  (3.23),  we  use  Eq.  (3.20)  for 
co^(  and  Eq.  (3.27)  for  m^n  and  m„^i  to  obtain 


1  I    »?ii 


®p/l2m,^.|l 


(3.28) 


S^E^Voc  +  Vov) 
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We  have  generalized  the  expression  of  Eq.  (3.28)  for  a  multi-subband  QW  system.  But  for 
a  single-subband  QW  system  in  both  the  CB  and  VB,  which  is  an  assumption  always  made 
in  [Cho94,  Hau90],  we  have  nj=Pj  =  1x20  and  therefore  the  summation  over;  can  be  sim- 
plified as 


"iD^o 


X  ("foe  +  Py^ov)    =  «2D  ("oc  +  "oJ 
J 

where  it  can  be  shown  that 


a„„  +  a„„  = 


e-      \"'cj\\ 
47tS2e  E     1 
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(3.29) 


(3.30) 


Finally  with  the  use  of  Eq.  (3.29),  Eq.  (3.28)  reduces  to 


"pi 


87C/I.,„ 


(3.31) 


For  a  multi-subband  QW  system,  which  is  the  real  case  of  the  two  QW  laser  systems  dis- 
cussed in  Chapter  2,  it  can  be  shown  that  Eq.  (3.31)  still  gives  a  close  numerical  result  to 
that  Eq.  (3.28)  gives.  To  simplify  the  calculation,  Eq.  (3.31)  will  be  used  hereinafter  for 
the  multi-subband  system  as  well.  Plugging  Eq.  (3.31)  into  Eq.  (3.23),  we  have 


1  +  ^  +  ^11 
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(3.32) 


Finally,  using  Eq.  (3.32)  in  Eq.  (3.19)  for  e,,  and  then  plugging  e,,  into  Eq.  (3.18),  the 
screened  Coulomb  potential  can  be  expressed  as 
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With  this  expression  for  V    ,  we  can  discuss  the  effects  of  bandgap  renormalization  due  to 
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plasma  screening  next. 

3.3.2  Bandgap  Renormalization 

As  carrier  density  increases,  more  unoccupied  states  become  available  to  the  VB 
electrons.  As  a  consequence,  the  VB  electron  distribution  can  change  in  such  a  way  that 
the  Coulomb  repulsion  with  the  CB  electrons  is  more  effectively  screened.  This  slight  re- 
arrangement of  the  charge  carriers  results  in  decreases  in  the  single-particle  energies  and, 
as  a  consequence,  a  reduction  in  bandgap.  In  this  case,  the  total  bandgap  renormalization 
(BGR)  Afgc^  is  the  sum  of  Coulomb-hole  (CH)  self-energy  AE^h  and  screened-exchange 
(SX)  self-energy  AE^x-  The  first  contribution  AEqh  describes  the  energy  reduction  of  elec- 
trons (holes)  with  different  spin  by  avoiding  each  other  because  of  the  mutual  Coulomb  re- 
pulsion and  is  the  difference  between  screened  and  unscreened  Coulomb  potential  [Cho94, 
Hau90],  ■•  .,       ,,  . 

111 
Using  Eqs.  (3.17)  and  (3.33),  we  obtain 

4KaE 
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Eq.  (3.34)  can  be  then  rewritten  as 
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(3.35) 


(3.36) 


To  evaluate  the  summation  over  q||,  we  first  must  be  able  to  count  the  available  qg  states. 
For  a  QW,  within  the  plane  of  QW,  a  2D  vector  qn  sweeps  out  a  circle  of  radius  q^  in  q,, 
space,  as  shown  in  Fig.  3.3.  One  q^  state  occupies  a  square  of  area  in  ^n  space  equal  to 
4n^/A ,  which  is  just  like  that  a  kj  state  takes  up  a  square  area  of  4n'^/A  in  /t,,  space 
[Cor93].  Therefore,  number  of  qy  states  within  the  area  defined  by  increments  in  q,,  and 


50 


Figure  3.3  Distribution  of  qg  vector  states  in  q,i  space  for  a  QW.  One  qq  state  takes  up 
a  square  area  of  47tVA  indicated  by  a  mesh  of  daslied  lines.  The  shaded  area 
of  ?i|d(|)rf?ii  represents  the  infinitesimal  area  defined  by  an  increment  of  dtjn 
in  the  magnitude  of  qj  and  an  increment  of  d0  in  its  sweep  angle  0.  The 
number  N,  of  qy  states  within  this  shaded  area  is  just  the  number  of  the 
squares  within  it.  Every  qy  state  can  be  taken  into  account  by  an  integration 
of  Wj  over  the  whole  (7||  space. 


sweep  angle  (|>  (shown  as  the  shaded  area  in  Fig.  3.3)  can  be  written  as 

?iid(t)rfi?|| 


N. 


47C^ 

A 


(3.37) 


where  the  spin  degeneracy  is  not  included.  The  summation  over  qy  can  be  replaced  by  an 
integration  of  N^  over  the  whole  q^  space 

p|.2n?||d(t)d<ir|| 
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4n^ 
A 


(3.38) 


Since  the  function  to  be  summed  in  Eq.  (3.36)  is  independent  of  <[),  its  integration  overij)  is 
equal  to  2e  and  Eq.  (3.36)  becomes 
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To  perform  the  integration  in  Eq.  (3.39),  we  evaluate  the  integrand  at  two  extreme  cases: 
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Using  the  evaluation  in  Eq.  (3.40),  Eq.  (3.39)  can  be  approximated  as 
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The  second  contribution  AE^x  results  from  the  change  in  electron  (hole)  self-energies  due 
to  the  fact  that  electrons  (holes)  with  equal  spin  avoid  each  other  (Pauli  exclusion  principle) 
and  is  given  by  [Cho94,  Hau90] 


SX  Zj     sq,  ^J  cq,       J  hq,> 
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(3.42) 
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We  substitute  Eq.  (3.33)  for  the  screened  Coulomb  potential  V^^_  and  Eq.  (3.38)  for  the 
summation  over  qg  in  Eq.  (3.42)  to  obtain 
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where  a  factor  of  q^/X^  is  pulled  out  of  the  first  bracket  in  last  form  of  AE^x-  Often,  the 
weak-momentum  dependence  of  AE^x  is  neglected  and  at  the  subband  edge  q^^  =  k[^  with 
it|i  =  0  in  Eq.  (3.13),  then  we  have 
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3.3.3  Simulation  Results 


InGaAs  OW.  Plots  of  AE^h,  AE^  and  AE^q^  as  a  function  of  volume  carrier  density 
A? for  Cl-HHl  and  Cl-LHl  transitions  are  shown  as  the  dash-dotted,  dashed  and  solid  lines 
in  Fig.  3.4,  respectively.  Often,  if  the  weak  momentum  dependence  of  BGR  is  neglected, 
BGR  can  be  taken  into  account  by  rigidly  shifting  the  gain  and  spontaneous  emission  by  an 
amount  of  AEbgr,  where  the  negative  sign  on  the  values  of  AEgQ^  means  a  red  shift  of  spec- 
ti'a  to  lower  photon  energy.  The  corresponding  BGR-shifted  TE  gain  and  spontaneous 
emission  spectra  of  Figs.  3.1(a)  and  (b)  are  shown  in  Figs.  3.5(a)  and  (b).  As  discussed  later 
in  Chapter  4,  the  red  shift  in  gain  spectra  gives  an  opposite  direction  in  predicting  the  de- 
pendence of  lasing  wavelength  on  cavity  length. 

CdZnSe  OW.  Plots  of  A£c;/.  AEsx  and  AEbcr  as  a  function  ofN  for  Cl-HHl  and  Cl- 
LHl  transitions  are  shown  as  the  dash-dotted,  dashed  and  solid  lines  in  Fig.  3.6,  respective- 
ly. The  coiresponding  BGR-shifted  spectra  of  Figs.  3.2(a)  and  (b)  are  shown  in  Fig.  3.7(a) 
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and  (b)  where  the  resulting  peak  gain  location  is  well  below  E^,  ^  by  about  55  meV.  This 
situation,  which  will  be  discussed  later  in  Chapter  5,  was  often  mistakenly  elucidated  as  the 
evidence  for  existence  of  excitonic  gain  in  CdZnSe  QW  lasers  even  at  room  temperature 
[Ding94]. 

Comparison.  By  comparing  Fig.  3.4  with  Fig.  3.6,  we  find  that  A£jc«  vs.  N  curves 
are  about  the  same  for  the  InGaAs  QW  and  CdZnSe  QW,  although  their  AEch  and  Afij^ 
components  are  different.  This  could  be  an  coincidence,  however,  it  agrees  with  an  univer- 
sal formula  commonly  used  for  iiEg^g  that  goes  with  one-third  power  ofN  [Chi88,  Col95]. 

3.4  Coulomb  Enhancement  fCE"! 

3.4.1  Coulomb  Enhancement  Factor 

The  Coulomb  attraction  between  electrons  in  the  CB  and  holes  in  the  VB  tends  to 
keep  electrons  and  holes  within  the  vicinity  of  each  other.  This  increases  the  radiative  re- 
combination rate  and  results  in  an  enhancement  in  optical  gain  as  well  as  spontaneous  emis- 
sion rate.  To  investigate  the  effects  of  Coulomb  enhancement  on  optical  gain  g  (ha)  and 
spontaneous  emission  R^^,  (fico),  we  start  with  the  complex  optical  susceptibility  x  (So)), 
since  both  g  (Sco)  and  R^^,  (fico)  are  related  to  %  (Sco).  In  a  multi-subband  system,  the  total 
optical  susceptibility  x(Sco)  is  the  summation  over  the  contributions  Z,(^w)  from  the  al- 
lowed subband  transition  pairs  and  is  expressed  as 

JC(fiw)   =  ^Xy(»co),  (3.45) 

J 

where  the  subscript;  indicates  they'th  subband  transition  pair.  With  inclusion  of  Coulomb 
enhancement,  Xj (^w)  is  given  by  (in  C.G.S.  units)  [Cho94,  Hau90]  ;• 
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Figure  3.4  Bandgap  renomialization  as  a  function  of  carrier  density  W  for  CI -HHl  and 
Cl-LHl  subband  transitions  in  the  8  nm  Ino.isGao  gsAs/GaAs  QW  at  T=  300 
K.  The  solid  curves  are  the  total  bandgap  renormalization  (BGR),  the  dash- 
dotted  curves  are  the  Coulomb-hole  components  (CH)  and  the  dashed 
curves  are  the  screen-exchange  (SX)  components. 
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Figure  3.5  The  rigidly  shifted  (a)  TE  gain  and  (b)  spontaneous  emission  spectra  (solid 
curves)  for  the  8  nm  Ino.15Gao.85As/GaAs  QW  at  N  =  1.7  x  10"  cm'  and  T 
=  300  K  by  the  bandgap  renormalization  due  to  plasma  screening.  The 
dashed  curves  are  the  spectra  without  BGR  shifts. 
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Figure  3.6  Bandgap  renomialization  as  a  function  of  carrier  density  N  for  Cl-HHl  and 
Cl-LHl  subband  transitions  in  the  4  nm  Cdo.3Zno.7Se/ZnSo.06Seo.94  QW  at  T 
=  300  K.  The  solid  curves  are  the  total  bandgap  renormalization  (BGR),  the 
dash-dotted  curves  are  the  Coulomb-hole  components  (CH)  and  the  dashed 
curves  are  the  screen-exchange  (SX)  components. 


57 


'^gl.hh 


4000 

j 

1  ■  II  1  1  1  1  1  1  1  1  1  1  1 

1  iT 1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

r=300K              : 

3000 

'- 

i  > 

'  4nmCdZnSeQW  j 
W  =  5.4  X  10'*  cm-^      : 

-^ 

2000 

'■_ 

1 

a- 

. 

w/CS+BGR 

--     ■    w/cs          : 

1 

So 

1000 

0 

innn 

'- 

\ 

/ 

7\ 

k*           .',     _         

1 — 1 ....  1 ... , 

\    \ 

i\  1 1 1.  ii'ii  1 1 1 1 1 11 1 1 

2.25 


2.35  2.45  2.55 

fico  (eV) 

(a) 


2.65 


> 


8  10^ 


6W    - 


4  10^    - 


<*r      2  10^*    i- 


2.25 


.1111 

1  1  1  1  1  1  1  1  1  1  1  I  1  1  1  1 

-T-I    1    1    1    I    1    1    1    1    11    1    1    1    1 

r=300K 

4  nm  CdZnSe  QW 

-i 

7 

iV  =  5.4xlO"cm-3 

- 

w/CS+BGR 

w/CS 

2.35 


2.45 
ha  (eV) 

(b) 


2.55 


2.65 


Figure  3.7  The  rigidly  shifted  (a)  TE  gain  and  (b)  spontaneous  emission  spectra  (solid 
curves)  for  the  4  nm  Cdo.3Zno.7Se/ZnSo.06Seo  94  QW  at  W  =  5.4  x  10"  cm^^ 
and  T  =  300  K  by  the  bandgap  renormalization  due  to  plasma  screening. 
The  dashed  curves  are  the  spectra  without  BGR  shifts. 
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where  V  is  the  volume  of  QW,  e„  =  n^  is  the  background  dielectric  constant  of  QW,  d„ 
is  the  square  root  of  the  dipole  moment  matrix  element  at  transition  involving  in-plane 
wavevector  k,,,  1/(1  -  Qx^  is  the  Coulomb  enhancement  factor  in  Fade  approximation, 
^s,k,k'  is  the  screened  Coulomb  potential, /,*  and /a*  are  the  Fermi-Dirac  distributions  of 
electrons  and  holes  with  energies  e^t  and  e^  in  the  CB  and  VB  respectively  and  hb  is  the 
dephasing  or  broadening  factor  due  to  carrier  scattering.  Originally,  Eqs.  (3.45)  to  (3.47) 
are  derived  under  the  assumption  of  a  single-subband  system  and  solved  finm  the  Bethe- 
Salpeter  equation  for  the  only  one  existing  subband  transition  pair  [Cho94,  Hau90].  Since 
in  a  multi-subband  system  a  set  of  Bethe-Salpeter  equations  for  various  allowed  subband 
transition  pairs  can  be  solved  independently  for  each  allowed  subband  transition  pair  as 
long  as  the  screening  effect  is  considered  with  inclusion  of  all  the  subbands  (as  we  have 
done  in  Section  3.3.1)  [Cam94],  we  are  allowed  to  generalize  Eqs.  (3.46)  to  (3.48)  to  be 
used  in  Eq.  (3.37)  of  a  multi-subband  system. 

Now  we  want  to  expand  Eqs.  (3.47)  and  (3.48)  in  terms  of  the  parameters  that  we  are 
familiar  with.  We  begin  with  the  dipole  moment  matrix  element.  The  relationship  between 
the  transition  matrix  element  and  the  dipole  moment  matrix  element  in  M.K.S.  units  is  giv- 
en by  [Col95] 

e2|Mj.|2  =  Imla^eM'^,  (3.49) 

where  e^UP  is  the  dipole  moment  matrix  element  and  x  is  the  position  operator.  The  factor 
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of  2  in  the  front  counts  for  the  spin  degeneracy  buried  in  iMj-P  and  is  not  included  in  the 
reduced  density  of  states  p„j  given  in  Eq.  (2.33).  On  the  contrary,  in  [Col95]  the  spin  de- 
generacy is  taken  into  account  in  p„j  instead  of  \MJ  ^ .  From  Eq.  (3.49),  we  can  relate  the 
dipole  moment  matrix  element  d^^,  in  C.G.S.  units  to  iMyP  in  M.K.S.  units  as  follows 


d^ 


e  |Mj.|2 


47:e„ 


(3.50) 


where  the  factor  of  l/(47ce„)  is  the  converting  constant  between  C.G.S.  and  M.K.S.  units. 
By  the  following  substitutions  in  Eqs.  (3.47)  and  (3.48) 
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(3.53) 


where  Eq.  (3.33)  is  used  for  V  .  The  superscript '  indicates  the  parameters  involving  tran- 
sition energy  £^^  and  the  associated  wavevector  k,', .  We  also  use  the  hco  dependent  broad- 
ening factor  r„  defined  in  Eq.  (3.6)  in  the  CE  expression  above  instead  of  the  fio) 
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independent  broadening  factor  hS  usually  used  in  previous  literatures  discussing  many- 
body  effects  [Cho94,  Hau90].  Similar  to  Eq.  (3.38),  we  replace  the  summation  over  kj|  in 
Eq.  (3.53)  with 
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Since  ^n  is  6  dependent  as  defined  in  Eq.  (3.13),  the  integrand  can  not  be  pulled  out  of  in- 
tegration over  6.  So  far,  we  have  successfully  derived  the  Coulomb  enhancement  factor  to 
be  in  a  form  such  that  we  can  use.  Next,  we  will  show  how  it  can  corporate  into  the  optical 
gain  and  spontaneous  emission  functions. 


3.4.2  Gain  and  Spontaneous  Emission  Functions  •  .  r    ; 

The  optical  gain  g  (So))  is  related  to  the  complex  optical  susceptibility  x  (^w)  by 
[Cho94,Hau90] 


g(tia)  = 1  mag  Ix  {fid))] 

=  ^Real[ix(ti(o)] 


(3.57) 


where  the  notations  Imag  and  Real  refer  to  the  imaginary  and  real  parts  of  a  complex  pa- 
rameter, respectively.  Plugging  Eq.  (3.47)  into  Eq.  (3.46)  for  Xy(fiM),  we  can  rewrite  Eq. 
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(4.45)  as 


Applying  Eq.  (3.58)  to  Eq.  (3.57)  and  substituting  e^  with  n^ ,  we  have 
g(»0))   =  "^Real 


ly      d^v     4710)  (fek  +  fhk-'^') 


The  siunmation  over  k  can  be  replaced  by 


(3.58) 


(3.59) 


(3.60) 


where  the  lower  limit  of  the  integral  is  chosen  to  be  Egj  because  p„j  is  zero  for  E„  <  E^  for 
a  given  transition  pair.  Using  the  substitutions  in  Eqs.  (3.51)  and  (3.60),  Eq.  (3.59)  be- 
comes 
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It  can  be  shown  that 
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where  L-  (Sco  -  £„)  is  defined  in  Eq.  (3.5).  Using  Eqs.  (2.25)  and  (3.62)  in  Eq.  (3.61),  fi- 
nally the  optical  gain  with  Coulomb  enhancement  can  be  expressed  as 
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where  the  E  subscript  indicates  the  Coulomb  enhanced  gain.  If  ?y  is  zero,  the  optical  gain 
function  with  CE,  g^  (tua),  in  Eq.  (3.62)  will  reduce  to  its  counterpart  without  CE,  gg  (fifn), 
inEq.  (3.1). 
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Recalling  Eq.  (3.58),  %  (ha)  contains  the  factor  associated  with  the  Fertni-Dirac  dis- 
tributions which  can  be  subdivided  into  two  components  as 

(/.*  +  /«-!)  =  [/.*/m]  -  [  ( 1  -  fek)  ( 1  -  /m)  ]  •  (3-64) 

The  first  square  bracket  is  referred  as  the  emission  component,  since  it  is  the  joint  proba- 
bility of  finding  an  occupied  state  in  the  CB  and  an  unoccupied  state  in  the  VB.  The  second 
square  bracket  is  referred  as  the  absorption  component,  since  it  is  the  joint  probability  of 
finding  an  occupied  state  in  the  VB  and  an  unoccupied  state  in  the  CB.  Using  this  subdi- 
vision in  Eq.  (3.58),  %  (fi^)  can  be  subdivided  into  two  components  as  follows 
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471  Y 


^co) 
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where  x(Sto)  -X<m(*w)  -ZafcC^t^)-  The  subscripts  em  and  aft  indicate  emission  and  ab- 
sorpdon,  respectively.  The  spontaneous  emission  function  is  related  to  the  imaginary  of 
iem  i^^)'  which  can  be  seen  as  an  analogy  to  Eq.  (3.57).  Similar  to  the  derivation  of  Eq. 
(3.63),  the  Coulomb  enhanced  spontaneous  emission  function  R,j^  (ha)  can  be  expressed  as 


R^^(ha)  =«ea/   Xj^- 


%j^EJ 


-qij(ha,EJ 


XL. (ha- E  J 


1-1^ 


(E,,-ha) 


dE. 
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(3.66) 


As  we  recall,  the  emission  component  of  Eq.  (3.64)  appears  in  the  expression  for  single- 
particle  ^^p,  (ftto)  function  given  in  Eq.  (2.26).  For  the  spontaneous  emission  function,  the 
Coulomb  enhancement  factor  1/(1  -  ^,^)  is  slightly  different  from  that  of  the  gain  function 
by  replacing  |Mj|  in  the  denominator  and  |A/^|  in  the  numerator  in  Eq.  (3.55)  with  \M-^ 
and  \Mj\  .  As  a  consequence,  \Mt\  and  |A?j-|  cancel  with  each  otiier  since  they  are  both 
equal  to  J2/3\M\  as  shown  in  Eq.  (2.15).  Again,  the  spontaneous  function  R,pi;(ha)  with 
CE  included  in  Eq.  (3.66)  reduces  to  its  counterpart  R^^b  (^f")  without  CE  included  in  Eq. 


63 

(3.2)if9y  =  0. 

3.4.3  Simulation  Results 

InGaAs  OW.  Fig.  3.8(a)  shows  plots  of  gj  (Sco)  and  ggCftw)  vs.  ftco,  with  inclusion  of 
the  rigid  shift  due  to  BGR,  forW  =  1.7  x  10"  cm'.  In  this  case,  CE  increases  the  peak  ma- 
terial gain  ^  by  a  factor  of  1.4  and  blue  shifts  its  spectral  position  £  by  4  meV.  This  blue 
shift  is  very  important  in  predicting  the  slope  of  cavity  length  dependence  of  lasing  energy 
as  would  be  discussed  in  Chapter  4.  Fig.  3.8(b)  shows  plots  of  R^^g  (h(o)  and  R^^^  (ha>)  vs. 
fid),  with  inclusion  of  the  rigid  shift  due  to  BGR,  for  N  =  1.7  x  10"  cm''.  The  increase  in 
spontaneous  emission  due  to  CE  is  evident  and  is  also  important  in  predicting  the  threshold 
current  of  diode  lasers  (see  Chapter  4). 

CdZnSe  OW.  In  Fig.  3.9(a),  gg  (ha)  and  gg  (ha)  vs.  hay,  which  include  the  rigid  shift 
due  to  BGR,  are  plotted  as  the  dotted  and  solid  curves  for  W  =  5.4  x  10"  cm',  respectively. 
It  is  shown  that  CE  increases  g  by  a  factor  of  2.35  and  blue  shifts  £  by  6  meV.  Fig.  3.9(b) 
shows  plots  of /?j^8(ft(B)  and/Jjp£(Sa))  vs.  Sco,  with  inclusion  of  the  rigid  shift  due  to  BGR, 
forW  =  5.4xlO"cm-3. 

Comparison.  Comparing  the  TE  gain  spectra  in  Figs.  3.8(a)  and  3.9(a),  it  is  clear  that 
CE  is  stronger  in  the  wide  bandgap  II- VI  CdZnSe  QW  than  in  the  III-V  InGaAs  QW.  A 
stronger  Coulomb  attractive  effect  can  be  indicated  by  a  larger  exciton  binding  energy  A£„ 
in  die  semiconductors.  For  a  2D  QW  structure,  A£„  is  equal  to  AE^  where  E„  is  the  Ryd- 
berg  constant  defined  in  Eq.  (3.12).  By  plugging  the  material  parameters  listed  in  Tables 
2. 1  and  2.2,  A£„  of  the  n  =  1  HH  excitons  is  found  to  be  about  1 1  meV  in  an  Ino  jsGao  gjAs 
QW  while  it  is  about  46  meV  (about  a  factor  of  4)  in  a  Cdo.3Zno.7Se  QW. 

Comment.  In  Eq.  (3.46),  the  Coulomb  enhancement  factor  is  derived  in  the  Pad^  ap- 
proximation which  is  a  good  approximation  for  carrier  densities  above  the  so-called  Mott 
density  Nn,„„.  At  the  Mott  density,  the  exciton  binding  energy  vanishes  and  excitons  cease 
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Figure  3.8  The  inclusion  of  Coulomb  enhancement  on  (a)  TE  gain  and  (b)  spontaneous 
emission  spectra  (solid  curves)  for  the  8  nm  Ino.15Gao.85As/GaAs  QW  at  /V  = 
1.7  X  10"  cm-'  and  r  =  300  K.  The  dashed  curves  are  the  spectra  without 
CE. 
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Figure  3.9  The  inclusion  of  Coulomb  enhancement  on  (a)  TE  gain  and  (b)  spontaneous 
emission  spectra  (solid  curves)  for  the  4  nm  Cdo.3Zno,7Se/ZnSo.o6Seo.94  QW 
at  N  =  5.4  X  10'*  cm-^  and  7"  =  300  K.  The  dashed  curves  are  the  spectra 
without  CE. 
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to  exist  due  to  the  screening  of  charged  carriers  in  semiconductors  [Hau90,  Zim88].  The 
value  for  N^m  can  be  simply  found  from  a  plot  of  AEgQ^  vs.  N  at  which  \AEbqd\  is  equal  to 
A£„  [Zim88].  For  an  8  nm  Ino.isGao.gsAs  QW,  with  A£„  =  1 1  meV,  we  thus  find  /V„„„  « 
5  X  10"  cm-3  from  Fig.  3.4.  For  a  4  nm  Cdo.3Zno.7Se  QW,  with  A£„  =  46  meV,  N„^,  is 
found  to  be  about  2  x  10"  cm'  (or  8  x  10"  cm"^  in  sheet  density)  from  Fig.  3.6.  Because 
of  this  large  value  for  Nj^^^  in  II- VI  CdZnSe  QW,  it  has  always  been  a  suspicion  that  lasing 
in  CdZnSe  QW  lasers  is  excitonic  in  namre  [Ding92,  Ding93,  Ding94].  In  Chapter  4,  we 
will  show  that  even  for  the  CdZnSe  QW  lasers,  the  carrier  densities  required  to  achieve  las- 
ing thresholds  at  room  temperature  are  higher  than  the  Mott  density. 

3.5  Summary 

In  this  chapter,  we  have  discussed  many-body  effects,  such  as  carrier  scattering, 
bandgap  renormalization  and  Coulomb  enhancement,  on  the  optical  gain  and  spontaneous 
emission  spectra  of  the  semiconductor  quantum  wells.  We  have  also  demonstrated  that  CE 
is  in  fact  stronger  in  the  CdZnSe  QW  than  in  the  InGaAs  QW  structures. 

In  next  chapter,  we  will  use  these  gain  and  spontaneous  emission  functions  developed 
in  Chapters  2  and  3  to  make  predictions  on  laser  characteristics  at  threshold  and  compare 
them  with  experiments  using  both  InGaAs  and  CdZnSe  QW  lasers. 


CHAPTER  4 
LASER  THRESHOLD  CHARACTERISTICS 


4.1  Introduction 

In  Chapter  2  and  3,  models  of  optical  gain  and  spontaneous  emission  based  on  elec- 
tron-hole plasma  (EHP)  were  described.  In  the  so-called  conventional  EHP  theory,  only 
carrier  scattering  (CS)  and  bandgap  renormalization  (BGR)  are  taken  into  account 
[Asa93,Col95,  Cor93]  and  Coulomb  enhancement  (CE)  is  left  out  mostly  due  to  its  involve- 
ment with  complicated  computation.  In  the  case  of  the  wide  bandgap  lasers,  such  ignorance 
of  CE  can  be  a  mistake  since  it  has  been  reported  that  CE  should  be  more  important  in  blue- 
green  CdZnSe  and  InGaN  quantum  well  (QW)  lasers  [Cho95a,  Cho95b,  Ree95a,  Ree95b, 
Ree96]  than  in  infrared  InGaAs  QW  lasers.  In  fact,  the  calculations  in  Chapter  3  have 
shown  that  CE  has  larger  effects  on  the  shapes  and  magnitudes  of  the  optical  gain  and  spon- 
taneous emission  spectra  for  the  4  nm  CdojZnojSe  QW  than  for  the  8  nm  Ino.15Gao.85As 
QW. 

In  this  chapter,  we  will  investigate  effects  of  CE  on  the  simulation  of  laser  threshold 
characteristics  for  both  the  well-developed  narrow-bandgap  InGaAs  QW  lasers  and  the  rel- 
atively new  wide-bandgap  CdZnSe  QW  lasers.  In  order  to  see  if  the  inclusion  of  CE  in  the 
simulation  yields  improved  agreement  with  experiment,  a  number  of  InGaAs/GaAs  QW 
and  CdZnSe/ZnSSe  QW  lasers  were  fabricated,  characterized  and  compared  to  predictions. 

In  Section  4.2,  the  laser  devices  characterized  are  described.  Comparisons  between 
experiments  and  predictions  on  cavity  length  dependence  of  lasing  energy  and  threshold 
current  are  discussed  in  Section  4.3  and  4.4,  respectively.  Lasing  energy  dependence  on 
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temperature  is  investigated  in  Section  4.5  and  the  work  is  summarized  in  Section  4.6. 

4.2  Laser  Devices 

■"■  '*       '      '■.,■',  '^ 

4.2.1  InGaAsOW 

:.     ■^  'i  •  • 

The  epitaxial  layers  of  the  InGaAs/AlGaAs  GRIN-SCH  single  QW  diode  lasers  are 
described  in  section  2.2.1.  Using  standard  photolithographic  techniques,  100  |J.m  stripes  on 
500  (im  centers  were  defined  on  the  p*-GaAs  cap  layer.  The  p-side  of  the  wafer  was  then 
metalized  using  Au  pulsed  electroplating  technique.  After  thinning  the  n-substrate  down 
to  a  wafer  thickness  of  about  100  nm,  the  n-side  was  also  Au-electroplated,  as  shown  in 
Fig.4.1.  Thewaferwascleavedintobars  with  three  cavity  lengths  (L  =  500, 1000  and  1500 
Hm)  and  then  characterized. 

4.2.2  CdZnSe  OW 

The  epitaxial  layers  of  the  CdZnSe/ZnSSe  SCH  single  QW  diode  laser  are  described 
in  section  2.2.2.  Fig.  4.2  shows  the  5  jtm  wide  buried-ridge  index-guided  devices  fabricat- 
ed from  this  structure  where  the  p-contact  metal  stack  is  Pd/Pt/Au  annealed  at  200°C  for  15 
minutes  in  forming  gas  and  the  n-contact  metal  stack  is  Pd/Ge/Au  annealed  at  200  'C  up  to 
an  hour  in  forming  gas  [Hab97].  These  state-of-the-art  CdZnSe  QW  laser  devices  were 
provided  by  the  3M  Company.  Laser  bars  with  three  cavity  lengths  {L  =  740,  1220  and 
2010  |im)  were  characterized. 

4.3  Cavity  Length  Dependence  of  Lasing  Energy 

As  discussed  in  Section  3.4.3,  CE  has  effects  on  the  spectral  location  and  the  magni- 
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Figure  4. 1  Cross-section  sketch  of  the  InGaAs/AlGaAs  GRIN-SCH  single  QW  diode 
lasers.  Au  stripes  with  widths  of  100  nm  on  500  nm  centers  were  deposited 
on  the  p+-GaAs  cap  layer  using  pulsed  electroplating  technique.  After  thin- 
ning the  n-GaAs  substrate  down  to  a  wafer  thickness  of  about  100  |lm,  Au 
was  deposited  everywhere  on  the  n-GaAs  substrate. 
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Figure  4.2  Cross-section  sketch  of  the  5  nm  wide  buried-ridge  index-guided  CdZnSe/ 
ZnSSe  SCH  single  QW  laser  devices  provided  by  the  3M  Company.  The  p- 
contact  metal  stack  is  Pd/Pt/Au  annealed  at  200°C  for  15  minutes  in  forming 
gas  and  the  n-contact  metal  is  Pd/Ge/Au  annealed  at  200°C  up  to  an  hour  in 
forming  gas  [Hab97]. 
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tude  of  the  peak  of  a  gain  spectrum.  As  a  consequence,  the  lasing  wavelength  X  (or  lasing 
energy  E)  for  conventional  cleaved  facet  lasers  (Fabry-Perot  lasers)  should  depend  on  CE, 
since  they  lase  at  the  peak  of  the  gain  spectra.  Due  to  the  fact  that  an  accurate  prediction 
of  £  depends  on  a  precise  knowledge  of  various  QW  bandgap  detennining  parameters  such 
as  QW  material  composition  and  degree  of  biaxial  compression,  QW  thickness,  band  off- 
sets etc.,  it  is  not  practical  to  test  the  importance  of  CE  by  comparing  measured  and  calcu- 
lated E  values.  However,  it  can  be  shown  that  the  dependence  of  E  on  carrier  density  N  is 
relatively  insensitive  to  the  choice  of  these  bandgap-determining  parameters  but  is  sensitive 
to  CE.  Since  this  dependence  can  be  determined  easily  by  measuring  the  dependence  of  £ 
on  laser  cavity  length  L,  it  is  practical  to  determine  the  importance  of  CE  by  comparing  cal- 
culations of  the  £  vs.  L  dependence  with  experiment. 

4.3.1  Experiment 

InGaAs  QW.  The  arrangement  of  experimental  setup  is  sketched  in  Fig.  4.3.  The 
lasers  were  characterized  in  the  bar  form  on  a  Cu  probing  stage  at  T  =  300  K  using  2  usee 
pulses  at  1  kHz  repetition  rate  to  avoid  heating  effect.  The  laser  beam  from  the  laser  facet 
was  focused  using  a  spherical  lens  on  to  the  front  slit  of  the  monochromator  which  in  addi- 
tion has  two  plane  mirrors,  a  grating  and  a  rear  slit,  as  depicted  in  Fig.  4.4.  The  prime  inir- 
ror  brings  the  laser  beam  to  the  grating  and  the  secondary  mirror  images  the  diffracted  beam 
from  the  grating  to  the  rear  slit.  A  silicon  PIN-IOD  photo-detector  is  placed  after  the  rear 
slit  and  plugged  into  a  digital  voltmeter  to  measure  the  photo  voltage.  The  reading  of  the 
voltmeter  reaches  a  maximum  once  the  laser  beam  is  exactly  imaged  to  the  rear  slit  by  ro- 
tating the  grating  to  a  certain  angle  for  a  given  lasing  wavelength  X.  Since  the  angle  has 
been  calibrated  to  the  corresponding  wavelength  by  the  manufacture  of  the  monochroma- 
tor, the  value  of  X  is  easily  read  out.  At  each  cavity  length  L,  X  was  therefore  recorded  for 
several  lasers  and  converted  to  an  average  emission  energy  £  for  that  L  value  using  the  fol- 


lowing  relation 


72 


E  =  j,  (4.1) 


where  h  is  the  Plank  constant.  The  measured  dependence  of  E  on  1/L  for  the  InGaAs  QW 
lasers  is  plotted  in  Fig.  4.5.  As  shown,  E  increases  by  about  3  raeV  as  L  decreases  from 
1500  to  500  \un  (i.e.  1/L  increases  from  6.7  to  20  cm"'). 

CdZnSe  OW.  The  arrangement  of  experimental  setup  is  sketched  in  Fig.  4.6.  The 
lasers  were  characterized  in  the  bar  form  on  a  Cu  probing  stage  at  T  =  300  K  using  2  usee 
pulses  at  1  kHz  repetition  rate  to  avoid  heating  effect.  The  edge  emission  spectra  from  the 
laser  facet  were  imaged  using  a  microscope  objective  on  to  the  front  slit  of  the  optical  mul- 
tichannel analyzer  (OMA).  Basically  OMA  is  similar  to  a  monochromator  which  has  an 
array  of  1024  photodetectors  (1024  channels)  at  the  rear  image  plane  instead  of  a  rear  slit. 
Photon  counts  received  by  each  channels  are  sent  to  the  computer  and  reconstructed  as  a 
spectrum  with  each  channel  corresponding  to  a  specific  wavelength.  Spectra  obtained  from 
a  typical  CdZnSe  QW  laser,  below  and  just  above  lasing  threshold,  are  shown  in  Fig.  4.7. 
The  lasing  wavelength  X,  defined  at  the  peak  of  the  spectrum  just  above  lasing  threshold, 
was  recorded  for  several  lasers  at  each  cavity  length  L  and  converted  to  average  emission 
energy  E  for  that  L  value  using  Eq.  (3.1).  Fig.  4.8  shows  the  measured  dependence  of  £  on 
1/Z,  for  CdZnSe  QW  lasers.  In  this  case,  E  increases  by  about  6  meV  as  L  decreases  from 
2010  to  740  urn  (i.e.  1/L  increases  from  5  to  13.5  cm'). 


4.3.2  Prediction 


For  typical  cleaved-facet  lasers  (Fabry-Perot  lasers)  with  a  cavity  length  L,  the  fre- 
quency spacing  Av  between  longitudinal  modes  is  given  by  [Ver89] 

where  n  is  the  refractive  index  of  the  cavity.  For  example,  a  GaAs  laser  cavity  with  n  =  3.6 
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Figure  4.3  Illustration  of  experimental  setup  for  measuring  the  lasing  wavelength  for 
the  InGaAs  QW  diode  lasers.  The  lasers  were  tested  in  a  pulsed  mode  and 
edge  emission  from  the  laser  facet  was  imaged  to  a  monochromator  using  a 
spherical  lens. 
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Figure  4.4  Sketch  of  the  components  inside  the  monochromator  used  in  Fig.  4.3.  The 
prime  mirror  brings  the  laser  beam  imaged  at  the  front  slit  to  the  grating  and 
the  secondary  mirror  images  the  diffracted  beam  to  the  rear  slit  by  rotating 
the  grating  to  a  proper  angle  for  a  given  lasing  wavelength. 
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Figure  4.5  Measured  lasing  energies  E  of  various  laser  cavity  lengths  L  for  the  8  nm 
Ino  isGao.gsAs  QW  lasers.  An  increase  of  about  3  meV  in  E  is  found  as  ML 
increases  from  6.7  to  20  cm"'  (i.e.  L  decreases  from  1500  to  500  nm). 
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Figure  4.6  Illustration  of  experimental  setup  for  measuring  the  edge  emission  spectra 
for  the  CdZnSe  QW  diode  lasers.  The  lasers  were  tested  in  a  pulsed  mode 
and  edge  emission  from  the  laser  facet  was  imaged  to  an  optical  multichan- 
nel analyzer  (OMA)  using  a  microscope  objective. 
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Figure  4.7  Edge  emission  spectra  obtained  from  a  typical  4  nm  Cdo.3Zno.7Se  QW  diode 
laser.  The  spectrum  below  the  threshold,  which  is  magnified  10  times  in  the 
plot,  is  wide.  Just  above  the  threshold,  the  spectrum  becomes  very  narrow 
and  the  lasing  wavelength  X  is  defined  at  the  peak  of  the  spectrum. 
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Figure  4.8  Measured  lasing  energies  E  of  various  laser  cavity  lengths  L  for  the  4  nm 
CdojZnojSe  QW  lasers.  In  this  case,  E  increases  by  about  6  meV  as  L  de- 
creases from  2010  to  740  \im  (i.e.  1/Z,  increases  from  5  to  13.5  cm'').  The 
error  bars  show  the  variation  of  E  obtained  from  each  laser  bar  and  the  open 
circles  are  the  corresponding  average  values  of  £. 
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and  L  =  1000  |im,  Av  is  42  GHz  calculated  using  Eq.  (4.2)  which  corresponds  to  an  energy 
spacing  /lAv  of  0.17  meV.  This  spacing  is  relatively  small  compared  to  the  spectral  width 
of  the  gain  function  g(h(D)  which  is  typically  in  the  order  of  100  meV  at  threshold.  As  a 
consequence,  the  Fabiy-Perot  lasers  will  lase  at  the  peak  of  g(ftCD)  and  the  lasing  energy  E 
is  essentially  determined  by  the  spectral  location  of  the  peak  of  gifia).  The  value  of  E  is 
associated  to  L  by  the  fact  that  the  magnitude  of  the  peak  gain  g  ofg(,tUo)  required  for  lasing 
at  L  is  given  by  the  expression  [Chi88,  Col95] 

«  =  «-  =  f[W*«-]-  ^'-'^ 

where  T  is  the  transverse  mode  overlap  factor  which  can  be  calculated  using  the  software 
MODEIG,  R  is  the  modal  reflectivity  at  the  facets  and  a,  is  the  mode  loss  coefficient.  The 
value  of  a,  can  be  determined  from  an  experimental  plot  of  the  inverse  differential  quantum 
efficiency  lAii  vs.  L  [Col95,  Hsu96].  <     ''■■.''./' 

The  procedures  for  obtaining  £  as  a  function  of  L,  depicted  in  Fig.  4.9,  are  summa- 
rized as  follows:  ( 1 )  Specify  N  and  T;  (2)  Calculate  and  plot  the  spectral  gain  function  g(ha>) 
using  Eqs.  (2.3)  and  (2.5)  for  single-particle  model,  Eq.  (3.1)  for  the  inclusion  of  CS,  Eqs. 
(3.41)  and  (3.44)  for  the  addition  of  BGR,  and  Eq.  (3.63)  for  the  inclusion  of  CE;  (3)  De- 
termine the  spectral  location  E  of  the  peak  of  g{h(o)  and  the  peak  gain  g  of  g(fta));  (4)  Use 
this  value  of  g  in  Eq.  (4.3)  to  solve  for  L;  (5)  Plot  this  E  vs.  l/L  pair,  (6)  Vary  A'  and  repeat 
steps  (1)  through  (5).  Finally,  a  plot  of  E  vs.  l/L  curve  is  obtained. 

InGaAs  OW.  For  the  wide-stripe  InGaAs  QW  lasers  used  in  our  experiment,  the  lat- 
eral mode  overlap  factor  is  assumed  to  be  unity.  Using  T  =  0.022,  R  =  0.32  and  a,  =  1 .7 
cm-'  in  Eq.  (4.3),  the  estimated  g,,,  are  422,  595,  1113  cm'  for  lasers  with  L  =  1500,  1000 
and  500  nm,  respectively.  Using  the  CS-modified  gain  function  gg  (ha)  given  in  Eq.  (3.1), 
Fig.  4.10(a)  shows  the  calculated  TE  gain  spectra  for  N  =  1.2,  1.4,  1.6,  1.8  x  10"  cm-' 
which  produce  g  in  the  range  of  interest  (~  300  to  1200  cm').  The  plots  indicate  that  E  blue 
shifts  as  N  increases.  Using  the  procedures  described  in  Fig.  4.9,  we  translate  these  plots 
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to  an  £  vs.  ML  curve  in  Fig.  4.10(b).  An  increase  in  E  (blue  shift)  of  about  5  meV  is  ob- 
tained as L  decreases  from  1500  to  5(X)  |im  (i.e.  1/L  increases  from  6.7  to  20  cm').  The 
corresponding  change  in  the  experimental  data  is  about  3  meV,  a  factor  of  1.7  smaller.  If 
the  BGR  shown  in  Fig.  3.4  is  added  to  the  gain  spectra  of  Fig.  4. 10(a),  the  resulting  rigidly 
shifted  gain  spectra  in  Fig.  4. 1 1  (a)  show  essentially  no  shift  of  peak  gain  position  as  N  in- 
creases. As  a  consequence,  the  corresponding  E  vs.  1/L  curve  in  Fig.  4.1 1(b)  shows  that 
the  BGR  shift  essentially  cancels  the  CS  shift.  It  is  interesting  to  note  that  this  CS/BGR 
cancellation-type  effect  has  been  reported  previously  in  the  modeling  of  the  dependence  of 
£■  on  Z,  in  GaAs  SQW  lasers  [Che93,  Chi88].  It  is  suggested  in  those  works  that  BGR  is 
overestimated  or  some  other  many-body  effect  needs  to  be  taken  into  account.  We  expect 
that  CE  may  be  the  answer  to  this  problem.  Therefore,  we  use  the  CE-modified  gain  func- 
tion ^^(Sto)  given  Eq.  (3.63)  to  calculate  the  TE  gain  spectra  for  A'=  1.1, 1.2, 1.3, 1.4  x  10'* 
cm"'  as  shown  in  Fig.  4. 12(a).  It  clearly  shows  that  E  increases  with  N  and  the  correspond- 
ing E  vs.  \IL  curve  in  Fig.  4.12(b)  shows  a  blue  shift  in  E  of  about  3  meV  as  L  decreases 
from  1500  to  500  )im,  in  good  agreement  with  the  experimental  shift.  Consequently,  we 
conclude  from  this  comparison  that  CE  needs  to  be  included  in  predicting  the  dependence 
of  £  on  L  in  InGaAs  QW  lasers.  The  importance  of  CE  in  determining  the  dependence  of 
operating  wavelength  on  temperature  in  InGaAs/GaAs  vertical-cavity  surface-emitting  la- 
sers (VCSEL)  has  been  reported  previously  [Cho95]. 

CdZnSe  OW.  For  the  CdZnSe  QW  lasers  used  in  our  experiment,  we  have  T  =  0.01, 
/?  =  0.243  and  otiS  lOcm'.  With  these  numbers,  the  estimated  g,A  are  1703, 2159  and  2910 
cm-'  fori  =  2010, 1220  and  740  nm,  respectively.  With  inclusion  of  CS,  the  calculated  TE 
gain  spectra  for  Af  =  6.5, 7.4,  8.3, 9.2  x  10"  cm"'  produce  g  in  the  range  of  interest  (-  1700 
to  3000  cm')  as  shown  in  Fig.  4.13(a).  The  corresponding  E  vs.  1/L  curve  in  Fig.  4.13(b) 
shows  an  increase  in  E  (blue  shift)  of  about  7  meV  as  L  decreases  from  2010  to  740  nm. 
While  this  is  close  to  the  6  meV  blue  shift  obtained  experimentally  (see  Fig.  4.8),  the  pre- 
dicted values  of  £  are  about  85  meV  higher  than  the  measured  ones.  If  the  BGR  shown  in 


Fig.  3.6  is  added  to  the  gain  spectra  of  Fig.  4. 1 3(a),  tiie  resulting  rigidly  shifted  gain  spectra 
in  Fig.  4.14(a)  show  a  red  shift  of  peak  gain  position  as  N  increases.  As  a  consequence,  Fig. 
4. 14(b)  shows  a  red  shift  of  the  predicted  E  with  decreasing  L,  opposite  to  the  observed  blue 
shift.  With  inclusion  of  CE,  the  TE  gain  spectra  for  N  =  5.0, 5.4, 5.7, 6.0  x  10"  cm'  in  Fig. 
4.15(a)  show  that  E  increase  with  N.  As  shown  in  Fig.  4.15(b),  one  obtains  a  blue  shift  in 
E  of  about  4  meV  as  L  decreases  from  2010  to  740  |im,  in  good  agreement  with  the  exper- 
imental shift.  It  also  shows  that  the  predicted  absolute  values  of  £  are  close  to  the  measured 
values,  differing  only  by  about  10  meV.  Consequently,  we  conclude  fttsm  this  comparison 
that  CE  needs  to  be  included  in  predicting  the  dependence  of  £  on  L  in  CdZnSe  QW  lasers. 

Comment  #1.  With  respect  to  the  discrepancy  between  the  measured  absolute  value 
of  E  and  the  prediction  with  CE  included,  we  believe  that  this  is  due  to  a  lack  of  precise 
knowledge  of  various  QW  bandgap  determining  parameters  as  stated  in  the  beginning  of 
this  section.  Since  the  dependence  of  £  on  L  is  relatively  insensitive  to  the  choices  of  these 
parameters,  we  believe  that  the  above  conclusions  about  the  importance  of  CE  in  both  the 
InGaAs  and  CdZnSe  QW  lasers  is  justified. 

Comment  #2.  In  Section  3.5,  we  discussed  a  little  bit  on  the  definition  of  Mott  den- 
sity ^moa  where  bound  excitons  cease  to  exist  and  the  Pad6  approximation  for  CE  factor 
works  better.  For  the  8  nm  Ino.isGaogsAs  QW  lasers,  N^^  is  much  less  than  5  x  10"  cm"' 
and  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers,  it  is  about  2  x  10"  cm'.  As  we  have  just  shown 
in  Figs.  4. 1 2  and  4. 1 5,  the  carrier  densities  required  to  achieve  lasing  action  in  the  InGaAs 
and  CdZnSe  QW  lasers  are  in  the  order  of  1  x  10"  cm''  and  5  x  10"  cm',  respectively, 
which  are  higher  dian  the  their  associated  N^„  values.  This  rules  out  the  possibility  of  ex- 
citonic  gain  at  room  temperature  even  in  CdZnSe  QW  lasers  for  which  excitons  were  said 
to  be  responsible  for  providing  material  gain  [Ding94].  In  Chapter  5,  we  will  show  more 
evidences  that  electron-hole  plasma,  not  excitons,  is  the  source  for  gain/absorption  in  the 
CdZnSe  QW  lasers  by  comparing  the  theoretical  calculations  using  model  ■with  CE  taken 
into  account  to  the  experiments. 
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Figure  4.9       Flow  chart  of  the  procedures  for  obtaining  a  theoretical  E  vs.  1/Z/  curve. 
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Figure  4. 10  With  CS  included,  (a)  the  calculated  TE  gain  spectra  for  the  8  nm 
Ino.15Gao.85As  QW  lasers  at  various  N  (x  10"  cm"')  indicated  by  the  num- 
bers. The  arrow  indicates  that  E  increases  with  N;  and  (b)  the  corresponding 
E  vs.  ML  curve.  An  increase  in  E  (blue  shift)  of  about  5  meV  is  obtained  as 
1/L  increases  from  6.7  to  20  cm'  (i.e.  L  decreases  from  1500  to  500  nm). 
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Figure  4. 1 1  With  CS  and  BGR  included,  (a)  the  calculated  TE  gain  spectra  for  the  8  nm 
Ino.isGao.ssAs  QW  lasers  at  various  N  (x  10"  cm"^)  indicated  by  the  num- 
bers. The  arrow  indicates  essential  no  shift  in  E  with  N\  and  (b)  the  corre- 
sponding E  vs.  ML  curve.  BGR  shift  essentially  cancels  the  CS  shift. 
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Figure  4. 12  With  CS,  BGR  and  CE  included,  (a)  the  calculated  TE  gain  spectra  for  the 
8  nm  Ino.15Gao.g5As  QW  lasers  at  various  N  (x  10"  cm'')  indicated  by  the 
numbers.  The  arrow  indicates  that  E  increases  with  N;  and  (b)  the  corre- 
sponding E  vs.  1/L  curve.  A  blue  shift  in  E  of  about  3  meV  is  obtained  as 
1/L  increases  from  6.7  to  20  cm'  (i.e.  L  decreases  from  1500  to  500  nm). 
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Figure  4.13  With  CS  included,  (a)  the  calculated  TE  gain  spectra  for  the  4  nm 
Cdo.3Zno7Se  QW  lasers  at  various  N  (x  10'*  cm"^)  indicated  by  the  numbers. 
The  arrow  indicates  that  E  increases  with  N;  and  (b)  the  corresponding  E  vs. 
VL  curve.  An  increase  in  E  (blue  shift)  of  about  7  meV  is  obtained  as  1/L 
increases  from  5  to  13.5  cm'  (i.e.  L  decreases  from  2010  to  740  nm). 
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Figure  4. 14  With  CS  and  BGR  included,  (a)  the  calculated  TE  gain  spectra  for  the  4  nm 
Cdo.3Zno.7Se  QW  lasers  at  various  N  (x  10'*  cm'')  indicated  by  the  numbers. 
The  arrow  indicates  that  E  decreases  with  N;  and  (b)  the  corresponding  £  vs. 
1/L  curve  and  the  measured  data.  The  predicted  shift  in  E  with  1/L  is  in  the 
opposite  direction  to  the  measured  shift. 
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Figure  4.15  With  CS,  BGR  and  CE  included,  (a)  the  calculated  TE  gain  spectra  for  the 
4  nm  Cdo.3Zno.7Se  QW  lasers  at  various  N  (x  10"  cm'^)  indicated  by  the 
numbers.  The  arrow  indicates  that  E  increases  with  N\  and  (b)  the  corre- 
sponding E  vs.  1/L  curve  with  the  measured  data.  A  blue  shift  in  E  of  about 
4  meV  is  obtained  as  1/L  increases  from  5  to  13.5  cm'. 
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Comment  #3.  In  previous  work  on  the  theory  of  CdZnSe  QW  lasers  [Cho95b],  it  was 
predicted  that  E  should  shift  to  lower  energies  (longer  wavelengths)  as  N  increases  from  4 
to  5  X  lO'*  cm'-'  in  an  8  nm  thick  QW.  In  this  work,  we  report  on  E  vs.  L  measurements 
using  4  nm  thick  CdZnSe  QW  buried  ridgeguide  lasers  [Hsu97b],  which  indicate  that  E 
shifts  to  higher  energies  (shorter  wavelengths)  as  N  increases.  In  addition,  we  use  a  single 
particle  model  [Cor93]  modified  by  the  addition  of  CE  [Hau90,  Cho94],  carrier  scattering 
[Asa93,  Yam87]  and  bandgap  renormalization  [Cho94,  Hau90]  to  calculate  a  dependence 
of  Eon  N  which  is  in  good  agreement  with  the  experimental  results. 

4.4  Cavity  Length  Dependence  of  Threshold  Current 

Since  the  threshold  gain  increases  with  decreasing  cavity  length,  the  threshold  cur- 
rent should  also  change  as  a  function  of  cavity  length.  By  investigating  the  cavity  length 
dependence  of  threshold  current,  we  are  able  to  demonstrate  the  true  mechanism  behind  the 
optical  gain  and  spontaneous  emission  of  semiconductor  lasers. 

4.4.1  Experiment 

The  experimental  geometry  is  shown  schematically  in  Fig.  4.16.  The  laser  bar  is  put 
on  probing  stage  with  p-side  up.  The  probe  and  stage  are  connected  to  "+"  and  "-",  respec- 
tively, of  a  HP  214A  pulse  generator  via  co-axial  cables  with  a  current-varying  pot  in  the 
circuit.  A  PIN-IOD  photo-detector  is  used  to  measure  the  laser  output  power  from  one  fac- 
et. The  measured  voltage  and  current  pulses  in  the  circuit  alone  with  output  power  pulse 
are  sent  to  three  separate  channels  of  a  Stanford  Research  System  (SRS)  where  a  pulse  sig- 
nal is  averaged  over  a  very  short  sampling  gate  placed  in  the  middle  of  the  pulse  duration 
and  converted  to  a  DC  signal  representing  the  averaged  value.  Then  a  set  of  DC  signals, 
representing  voltage,  current  and  output  power,  are  sent  to  a  personal  computer  (PC)  for 
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further  data  processing  and  a  real-time  output  power  P^  vs.  current  /  curve  is  displayed  on 
the  PC  monitor.  By  extrapolating  the  linear  region  of  P„  vs.  /  curve  back  to  the  /  axis,  the 
cross-over  is  defined  as  the  threshold  current  /,/,.  The  voltage,  current  and  power  pulse  sig- 
nals can  also  be  monitored  simultaneously  on  an  oscilloscope  and  compared  to  the  reading 
shown  on  PC  monitor. 

InGaAs  OW.  Fig.  4.17  shows  the  measured  /,*  vs.  L  for  the  8  nm  Ino.jsGao.gsAs  single 
QW  lasers  with  100  Jim  stripe  width.  The  open  circle  is  the  averaged  /,;,  for  the  good  laser 
devices  tested  from  the  same  laser  bar  and  the  error  bars  represent  the  span  of  variation  in 

CdZnSe  OW.  Fig.  4. 1 8  shows  the  measured  7,^  vs.  L  for  the  4  nm  Cdo jZiiq  7Se  single 
QW  lasers  with  5  nm  wide  buried-ridge. 

4.4.2  Prediction  '"'       ••^■-'^•' 

In  Section  4.3,  we  have  shown  that  a  certain  value  of  A'  in  QW  is  needed  to  provide 
a  peak  material  gain  g  required  for  lasing  at  a  given  cavity  length  L.  For  diode  lasers,  N  is 
provided  by  injecting  current  through  metal  contacts,  as  depicted  in  Fig.  4.19,  where  elec- 
trons and  holes  are  injected  to  the  active  layer  (or  QW  in  the  laser  structures  studied)  from 
n-  and  p-side,  respectively,  and  then  recombine  with  each  other.  Under  steady  state  condi- 
tion, rate  of  injected  electrons  into  QW  is  equal  to  the  rate  of  recombining  electron  in  QW 
which  can  be  expressed  mathematically  as 

where  T);  is  the  injection  efficiency  which  is  defined  as  the  fraction  of  current  being  injected 
into  QW,  /  is  the  injected  current,  w  is  the  stripe  width,  L  is  the  cavity  length,  d,  is  the  QW 
thickness  and  /f,„  is  the  rate  of  recombining  electron  per  unit  volume  (in  units  of  s'cm'^). 
The  recombination  processes  include  a  spontaneous  recombination  rate  R,.,  a  nonradiative 


91 


photo-detector 


current  control 


VScI 
measuring  box 


detector 
battery  box 


I    V   P. 


n    laser  diode 


SRS  boxcar 


*•/ 


V    '-  J 


computer 


pulser 


P  -- --!.x....  ^ 

■  , .^— © 


oscilloscope 


Figure  4.16  Illustration  of  experimental  setup  for  measuring  the  output  power  /*„  as  a 
function  of  injection  current  /  for  diode  lasers.  By  extrapolating  the  linear 
region  of  P,,  vs.  /  curve  back  to  the  /  axis,  the  cross-over  is  defined  as  the 
threshold  current  /,l. 
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Figure  4. 17  Measured  threshold  currents  I,^  of  various  laser  cavity  lengths  L  for  the  8 
nm  Ino.isGao.ssAs  QW  lasers.  The  error  bars  show  the  variation  ofI,h„  ob- 
tained from  each  laser  bar  and  the  open  circles  are  the  coiresponding  aver- 
age values  of /,j„. 
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Figure  4.18  Measured  threshold  currents  /,^„  of  various  laser  cavity  lengths  L  for  the  4 
nm  Cdo,3Zno.7Se  QW  lasers.  The  error  bars  show  the  variation  of  I,i^  ob- 
tained from  each  laser  bar  and  the  open  circles  are  the  corresponding  aver- 
age values  of /,^„. 
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Figure  4.19  Schematic  sketch  of  a  current-injected  diode  laser  device.  Holes  and  elec- 
trons are  injected  from  the  p-side  and  n-side,  respectively,  and  then  radia- 
tively  recombine  with  each  other  in  the  active  layer. 
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recombination  rate  R„,  a  carrier  leakage  rate  R,  and  a  stimulated  recombination  rate  R„. 
Thus,  we  can  write  [Col95] 

For  a  laser  below  or  at  threshold  where  the  photon  density  is  relatively  low,  R„  is  negligible. 
With  R„  -  0  and  using  Eq.  (4.5)  in  Eq.  (4.4),  we  obtain 

^=   («»,  +  '?„, +  «/)>^^^.-  (4.6) 

Since  the  exact  dependence  oiR„  and  Ri  on  /  is  unknown,  a  conventional  approach  to  deal 
with  this  problem  is  to  define  a  radiative  officiating  T),  as 

n  = -^ .  (4.7) 

Then,  substituting  R,^  +  R„  +  /?,  in  Eq.  (4.6)  in  terms  of  R,p  and  Tl;.  using  Eq.  (4.7),  we  can 
rewrite  Eq.  (4.6)  as 

Ti,.Ti/  =  ed^R    wL,  (4.8) 

which  means  that  a  fi:action  (=  tI/TI,.)  of  /  gives  rise  to  the  spontaneous  recombination  rate 
in  QW.  For  the  convenience,  we  define  a  radiative  current  density  J,ad  2S  the  component 
of/  which  gives  rise  to  the  spontaneous  recombination 

•/r.<i^^^n,n/.  (4.9) 

Combining  Eqs.  (4.8)  and  (4.9),  we  obtain 

Jrad^ed^'^sp-  (4-10) 

Recalling  that  Rsp(.h(a)  is  the  spontaneous  recombination  rate  per  unit  volume  per  unit  en- 
ergy at  hat  (in  units  of  s'cm-^eV'),  R^^  (in  units  of  s'cm'^)  is  then  just  the  integration  of 
R^itiai)  over  aU  to  [Col95] 

R.p  =  ^^R.p(f^^)'i(^'^)-  (4.11) 

As  a  consequence  of  using  Eq.  (4. 11 )  in  Eq.  (4.9),  7^^  can  be  expressed  as 
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-^,..  =  ^^.J^«..  («")''(''")•  .,  (4.12) 

Specifying  N  and  T,  Rspiha)  can  be  calculated  using  Eq.  (3.2)  with  inclusion  of  CS  (without 
CE)  and  Eq.  (3.63)  with  inclusion  of  CE.  Effect  of  CE  on  J^^  can  thus  be  calculated  by 
using  Eq.  (4.12). 

4.4.2. 1  Gain-current  relation 

The  first  step  in  predicting  /,^  for  a  given  laser  is  to  determine  the  value  of  Z,^  re- 
quired to  achieve  g,i,  for  that  laser.  Specifying  N  and  temperature  T,  we  can  calculate  gifiai) 
and  Rsp(fUo)  based  on  the  model  chosen  (see  Chapter  2  and  3).  The  peak  material  gain  g 
can  therefore  be  found  from  g(fi(o)  and  by  using  Eq.  (4.12),  J^^  can  be  calculated  firom 
R^p(h(n).  As  a  consequence,  a  g-Jrad  pai""  is  generated  for  a  given  N.  By  repeating  this  pro- 
cedure for  various  N,  which  is  illustrated  in  Fig.  4.20,  a  g-J,^  curve  (gain-current  relation) 
is  created  for  a  given  laser  structure.  The  cross-over  of  a  g-J,ad  curve  at  the  J,^  axis  is  de- 
fined as  the  transparency  current  density  /„  where  the  material  becomes  transparent  for  g 
=  0. 

InGaAs  OW.  Fig.  4.21  shows  the  g-J,ad  curves  calculated  with  and  without  CE  for 
the  8  nm  lno.15Gao.85As  QW  lasers.  The  model  with  CE  included  predicts  70  A/cm^  for  J„ 
while  the  model  without  CE  included  predicts  65  A/cm^  for  /„.  For  low  J„ui  values,  the 
model  with  CE  predicts  g  values  lower  than  those  without  CE.  For  7,^  values  higher  than 
the  cross  point  of  the  two  curves,  the  model  with  CE  predicts  g  values  a  bit  higher  than 
those  without  CE.  This  difference  is  important  in  predicting  the  threshold  current  /,j  dis- 
cussed in  the  next  subsection. 

CdZnSe  OW.  Fig.  4.22  shows  the  g-Jrad  curves  calculated  with  and  without  CE  for 
the  4  nm  Cdo.3Zno.7Se  QW  lasers.  The  model  with  CE  included  predicts  150  A/cm^  for/,, 
while  the  model  without  CE  included  predicts  120  A/cm^  for  /„..  For  low  J,^  values,  the 


97 

model  with  CE  predicts  g  values  lower  than  those  without  CE.  For  /,^  values  higher  than 
the  cross  point  of  the  two  curves,  the  model  with  CE  predicts  g  values  much  higher  than 
those  without  CE.  This  effect  of  CE  on  the  g-J,^  relation  will  be  shown  to  be  very  impor- 
tant in  predicting  the  threshold  current  /^  accurately  discussed  in  the  next  subsection. 

Comparison.  Comparing  Fig.  4.21  with  Fig.  4.22,  we  observe  that  CE  obviously  has 
a  much  more  profound  effect  on  the  g-J^^  relation  for  the  CdZnSe  QW  lasers  than  for  the 
InGaAs  QW  lasers.  This  again  demonstrates  that,  as  expected,  CE  is  stronger  in  wide  band- 
gap  materials  than  in  narrow  bandgap  materials. 

As  mentioned  in  Section  2.6,  a  CdZnSe  QW  laser  structure  has  a  higher  transparent 
carrier  density  N„  than  an  InGaAs  QW  laser  structure,  resulting  from  the  fact  that  the 
CdZnSe  QW  has  heavier  carrier  effective  masses.  This  larger  A',r  shows  up  as  higher  /„ 
values  in  g-J,ad  curves  for  the  CdZnSe  QW  than  the  InGaAs  QW  (about  twice  higher),  as 
shown  in  Figs.  4.21  and  4.22.  It  looks  like  that  the  CdZnSe  QW  lasers  will  not  be  as  good 
as  the  InGaAs  QW  lasers,  since  they  require  higher  transparency  to  achieve  gain  in  the  QW. 
However,  the  strong  CE  effect  in  the  CdZnSe  QW  lasers  results  a  very  steep  g-J^ad  curve, 
as  shown  in  Fig.  4.22,  which  indicates  that  the  CdZnSe  QW  lasers  can  be  as  good  as  the 
InGaAs  QW  lasers  or  ever  better.  To  understand  why  CE  increase  the  slope  of  a  g-J,^ 
curve,  we  plot  g  and  J,aj  as  a  function  ofN  in  Figs.  4.23  and  4.24  for  the  8  nm  I%,5Ga<,  gsAs 
and  the  4  nm  Cdo.3Zno7Se  QW  lasers,  respectively.  As  shown,  the  enhancements  in  g  due 
to  CE  are  higher  than  those  in  y„j,  especially  for  the  CdZnSe  QW  lasers.  As  a  conse- 
quence, the  slope  of  a  calculated  g-J^ad  curve  is  larger  with  CE  included  than  without  CE. 

4.4.2.2  Predicting  threshold  current 

Since  g,/,  is  known  for  a  given  L  through  Eq.  (4.3),  the  predicted  value  of  J,^  at 
threshold,  defined  as  7^^  ,/,,  can  be  determined  from  the  g-Jrad  curve,  as  depicted  in  Fig. 
4.20.  To  the  extent  that  current  and  carrier  spreading  beyond  the  stripe  area  (lateral  leak- 
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Figure  4.20     Flow  chart  of  the  procedures  for  obtaining  a  theoretical  g  vs.  J,^  curve. 
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Figure  4.21  Effect  of  CE  on  the  relation  between  peak  gain  g  and  radiative  current  den- 
sity J^  for  the  8  nm  Ino.15Gao.85As  QW  lasers.  The  solid  line  is  calculated 
with  inclusion  of  CE  and  the  dotted  Une  is  calculated  without  CE. 
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Figure  4.22  Effect  of  CE  on  the  relation  between  peak  gain  g  and  radiative  current  den- 
sity J,^  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers.  The  solid  line  is  calculated 
with  inclusion  of  CE  and  the  dotted  line  is  calculated  without  CE. 
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Figure  4.23  Effects  of  CE  on  peak  gain  g  and  radiative  current  density  J,^  for  the  8  nm 
Ino.15Gao.85As  QW  lasers,  (a)  g  as  a  function  oiN\  and  (b)  7,^  as  a  function 
of  M  The  solid  lines  are  calculated  with  inclusion  of  CE  and  the  dotted  lines 
are  calculated  without  CE. 
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Figure  4.24  Effects  of  CE  on  peak  gain  g  and  radiative  current  density  J,^  for  the  4  nm 
Cdo.3Zno.7Se  QW  lasers,  (a)  g  as  a  function  of  N;  and  (b)  Jy^d  as  a  function 
of  W.  The  solid  lines  are  calculated  with  inclusion  of  CE  and  the  dotted  lines 
are  calculated  without  CE. 
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age)  can  be  ignored,  J^^d,  ,1,'^^-  had,  :h  is  the  predicted  fraction  of  total  current  giving  rise 
to  spontaneous  emission  in  the  QW  at  threshold.  Using  Eq.  (4.9),  we  obtain  the  predicted 
value  of  /,(i  for  a  given  L 

where  Tl,  is  assumed  to  be  unity  in  our  lasers. 

If  the  high  conductivity  p*  cap  layer  was  not  removed  in  the  area  outside  the  stripe 
region,  like  the  InGaAs  QW  lasers  studied  in  these  experiments  (see  Fig.  4.1),  the  lateral 
leakage  current  /;t  can  be  significant  [Gro94].  As  a  consequence,  7,^  must  be  subtracted 
from  the  measured  value  of /,j,  that  is  /,a„,  in  order  to  compare  with  the  value  of/,,,  predicted 
using  Eq.  (4. 13).  However,  when  /« is  significant,  the  "measured"  value  of  Tli  =  J\^,  deter- 
mined from  the  \lT\d  vs.  L  plot,  must  also  be  corrected  IHu94].  It  can  be  shown  that  when 
/,t  does  not  contribute  to  the  laser  power  (true  for  our  lasers  where  the  lateral  width  of  the 
laser  emission  at  the  facet  is  about  the  same  as  the  stripe  width  w  =  1(X)  |im),  then  the  pre- 
dicted value  of  threshold  current  /^^  is  given  by 

'..-^-  (4.14) 

As  a  consequence,  one  can  compare  7,^^  to  I,km  without  actually  determining  In,. 

InGaAs  QW.  Given  ri,„  =  0.65  and  Tj,  =  1,  the  7,^^  vs.  L  curves  calculated  with  and 
without  CE  for  the  8  nm  Ino.isGa^^^As  QW  lasers  are  plotted  in  Fig.  4.25  along  with  the 
measured  I,^  at  three  cavity  lengths.  It  is  clear  that  the  model  with  CE  does  a  better  job 
than  the  model  without  CE  in  predicting  7,j„. 

CdZnSe  QW.  Given  T|i„  =  0.69  and  t],  =  1,  plots  of  7,,,^  as  a  function  of  L  calculated 
with  and  without  CE  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers  are  shown  in  Fig.  4.26  along  with 
the  measured  I,^  at  three  cavity  lengths.  The  model  without  CE  obviously  overestimate 
the  threshold  currents  by  about  40%  to  50%  of  the  measured  threshold  currents.  With  in- 
clusion of  CE,  the  prediction  shows  improved  agreement  with  the  experiment. 

'  '  ;         -,  ■        1   ;  "l  ?  ^;  ' 
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Figure  4.25  Effect  of  CE  on  the  predicted  threshold  current  I,i^  as  a  function  of  cavity 
length  L  for  the  8  nm  Ino.ijGao  jsAs  QW  lasers.  The  solid  line  is  calculated 
with  inclusion  of  CE  and  the  dotted  line  is  calculated  without  CE.  Measured 
threshold  currents  /,^„  of  three  cavity  lengths  L  =  500,  1000  and  1500  nm 
are  also  shown  as  the  open  circles  with  error  bars. 
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Figure  4.26  Effect  of  CE  on  the  predicted  threshiold  current  7,^^  as  a  function  of  cavity 
length  L  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers.  The  solid  line  is  calculated 
with  inclusion  of  CE  and  the  dotted  line  is  calculated  without  CE.  Measured 
threshold  currents  I,f^  of  three  cavity  lengths  L  =  740,  1220  and  2010  nm 
are  also  shown  as  the  open  circles  with  error  bars. 
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4.5  Temperature  Dependence  of  Lasing  Energy 

4.5.1  Experiment 

A  chip  of  the  4  nm  Cdo.3ZnQ7Se  QW  laser  with  the  structure  shown  in  Fig.  4.2  was 
cleaved  from  the  laser  bar  with  L  =  1219  (im.  Then  using  EPO-TEK  H20E  silver  epoxy, 
the  laser  chip  was  mounted  with  p-side  up  on  the  Cu-block  and  Au-wire  bonded  to  one  of 
the  Au  probing  pads  evaporated  on  the  tops  of  two  ceramic  standoffs,  as  illustrated  in  Fig. 
4.27(a).  To  avoid  thermally  degrading  the  CdZnSe  QW  laser,  this  chip-on-block  (COB) 
package  was  left  in  the  room  temperature  for  36  hours  and  then  baked  in  50°C  for  1  hour 
to  cure  the  silver  epoxy.  As  sketched  in  Fig.  4.27(b),  the  COB  was  loaded  to  a  Cu-fixture 
and  a  Cu-sheet,  which  was  tightened  with  two  plastic  screws  to  the  Cu-fixture  with  a  ce- 
ramic board  in  between,  was  used  to  make  contact  with  the  two  standoffs  of  the  COB.  The 
whole  fixture  was  loaded  into  the  Laser  Photonics  L5736  cryostat  and  the  wire  soldered  on 
the  Cu-sheet  was  connected  to  the  circuit  board  inside  the  cryostat  which  was  externally 
connected  to  the  HP214A  pulse  generator  and  the  Laser  Photonics  L5720  temperature  con- 
troller. After  the  cryostat  was  sealed  and  vacuumed  down  to  20  mtorr,  the  temperature  in- 
side the  cryostat  was  cooled  down  using  liquid  nitrogen  and  controlled  by  the  temperature 
controller.  The  laser  device  was  characterized  using  1  |isec  wide  pulses  at  a  1  kHz  repeti- 
tion rate.  The  laser  beam,  emitting  from  one  of  cleaved  facet  and  passing  through  the  quartz 
window  of  the  cryostat,  was  imaged  to  the  OMA  to  record  the  spectra  of  the  light  output 
and  determine  the  lasing  wavelength,  just  like  the  one  shown  in  Fig.  4.7.  Fig.  4.28  shows 
the  measured  values  of  the  lasing  energy  £  at  various  temperatures  T  in  the  range  from  82 
to  3(X)  K.  It  is  found  that  as  T  reduces,  E  shifts  to  higher  energies  (a  blue  shift).  Using  a 
linear  fit  to  the  data  points  (plotted  as  the  dotted  lines),  the  increasing  rate  of  E  with  de- 
creasing T  is  found  to  be  about  0.48  meV/K  for  temperatures  between  160  K  and  300  K 
while  for  temperatures  below  160  K  and  above  82  K,  the  increasing  rate  is  smaller  and 
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about  0.32  me V/K. 
4.5.2  Prediction 

Using  Eq.  (4.3),  the  required  threshold  material  gain  g,/,  for  the  4  nm  CdojZno^Se 
QW  laser  with  L  =  1219  \im  is  estimated  to  be  around  2159  cm'.  Using  the  model  with 
many-body  effects  included,  the  TE  gain  spectra  at  various  temperatures  with  peak  gain 
value  reaching  the  required  g,i,  are  calculated.  The  TE  gain  spectra  in  Figs.  4.29(a)  are  cal- 
culated with  inclusion  of  CS  at  the  threshold  carrier  densities  N,^  =  7.4,  5.4  and  3.2  x  10" 
cm"^  for  T  =  300,  240  and  160  K,  respectively,  and  under  the  same  conditions,  the  TE  gain 
spectra  in  Figs.  4.29(b)  are  calculated  with  inclusion  of  CS  and  BGR.  In  addition  to  CS  and 
BGR,  the  TE  gain  spectra  in  Figs.  4.29(c)  are  calculated  with  inclusion  of  CE  at  the  N,,,  = 
5.23,  3.65,  2.07  x  10'«  cm-^  for  T  =  300,  240  and  160  K,  respectively.  All  plots  show  that 
the  gain  spectra  blue  shift  to  higher  energies  as  the  T  decreases.  This  is  mainly  due  to  the 
fact  that  the  bandgap  of  semiconductor  material  increases  with  decreasing  temperatiu^e. 

As  discussed  in  section  4.3.2  that  E  is  determined  by  the  spectral  location  of  the  peak 
gain  at  threshold.  Fig.  4.30  shows  plots  of  £  as  a  function  of  T  determined  from  the  calcu- 
lated threshold  gain  spectra  (see  Fig.  4.29)  at  temperatures  in  the  range  from  160  to  300  K. 
With  only  CS  included,  the  predicted  E  values  are  about  65  to  85  meV  higher  than  the  mea- 
sured E  values  and  a  blue  shift  of  49  meV  in  E  is  obtained  for  T  decreasing  from  300  K  to 
160  K.  The  measured  shift  of  E  within  the  same  range  of  T  is  found  to  be  about  67  meV 
from  Fig.  4.28,  a  factor  of  1 .37  larger.  With  BGR  also  included,  the  predicted  E  values  are 
close  to  the  measured  £  values,  but  the  corresponding  shift  of  £  with  T  is  about  76  meV. 
In  addition  to  CS  and  BGR,  the  model  with  CE  included  predicts  a  blue  shift  of  about  71 
meV  which  is  in  the  best  agreement  with  the  experiment  among  all  three  models.  As  a  con- 
sequence, it  is  important  to  include  CE  in  the  theory  in  predicting  the  dependence  of  lasing 
energy  on  temperature  between  160  and  300  K  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers. 
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Figure  4.27      Scliematic  diagrams  of  (a)  mounting  the  laser  chip  on  the  Cu-block  (COB) 
and  (b)  loading  the  COB  to  the  Cu-fixture. 
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Figure  4.28  The  measured  lasing  energy  £  as  a  function  of  temperature  T  for  the  4  nm 
Cdo.3Zno.7Se  QW  laser  with  L  =  1219  nm.  The  open  circles  are  the  mea- 
sured data  and  two  dotted  lines  are  linearly  fitted  to  the  data  over  two  tem- 
perature ranges,  between  82  and  160  K  and  between  160  K  and  300  K. 
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For  T  lower  than  160  K,  the  CE  model  reaches  its  limits  due  to  the  random  phase  ap- 
proximation (RPA)  made  in  considering  the  screening  effects  and  the  Pad6  approximation 
made  in  deriving  the  Coulomb  enhancement  factor,  as  discussed  in  Chapter  3.  It  is  inter- 
esting to  notice  that  160  K  is  also  the  temperature  at  which  the  measured  increasing  rate  in 
E  with  decreasing  T  slows  down.  This  may  imply  that  excitons  become  important  at  tem- 
perature below  160  K.  At  cryogenic  temperatures,  the  role  of  phase  transition  between  ex- 
citons and  electron-hole  plasma  in  the  stimulated  emission  processes  of  CdZnSe  QW  laser 
structures  has  been  discussed  extensively  by  Cingolani  et  al.  in  [Cin96]. 

4.6  Summary 

In  this  chapter,  the  cavity  length  dependences  of  lasing  energy  and  threshold  current 
were  both  studied  using  the  InGaAs  and  CdZnSe  QW  lasers.  It  has  been  demonstrated  that 
predictions  using  the  theory  with  Coulomb  enhancement  included  are  in  good  agreement 
with  the  experiments.  The  importance  of  Coulomb  enhancement  in  the  CdZnSe  QW  lasers 
was  also  justified  by  the  excellent  agreement  found  between  the  measured  temperature  de- 
pendence of  lasing  energy  and  the  prediction  with  CE  included  at  temperatures  between 
160  K  and  300  K. 


8  • .  \ 


•■;  /;»" 


111 


I 


3500 
3000 


I     I     I     I      I     I     I     I     I     I     I     I     I     ]     I     I     I     I     I     I     I     I     I    'I     I      I     I     I 

L  4  nm  CdZnSe  QW 
w/CS 
2500  E-  r=300K     240 K     160 K   -3 


g:h 


3500 
3000 
2500 
2000 


■  •  •  ■  I 

4  nm  CdZnSe  QW 
w/ CS+BGR 

r=300K 


240  K 


160  K 


g,i, 


Figure  4.29  Calculated  TE  gain  spectra  at  r  =  300,  240  and  160  K  for  the  4  nm 
CdfljZno.vSe  QW  laser  with  L  =  1219  |j.m  at  threshold,  where  the  required 
threshold  material  gain  g,^  is  about  2159  cm',  using  models  (a)  with  CS  in- 
cluded; (b)  with  CS  and  BGR  included;  (c)  with  CS,  BGR  and  CE  included. 
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Figure  4.30  Temperature  T  dependence  of  the  lasing  energy  £  for  the  4  nm  Cdo.sZno  7Se 
QW  laser  with  cavity  length  L  =  1219  |im.  The  open  circles  are  the  experi- 
mental measurements,  the  dotted  curve  is  calculated  with  inclusion  of  CS, 
the  dash-dotted  curve  is  calculated  with  inclusion  of  CS  and  BGR  and  the 
solid  curve  is  calculated  with  inclusion  of  CS,  BGR  and  CE. 


CHAPTERS 
ABSORPTION  RESONANCES  DUE  TO  COULOMB  ATTRACTION 


5.1  Introduction 

As  discussed  in  Chapter  3,  exciton  binding  energy  in  II- VI  CdZnSe  QW  is  larger  than 
in  ni-V  InGaAs  QW  for  its  heavier  electron  and  hole  effective  masses  and  smaller  dielec- 
tric constant,  exciton  lasing  mechanism  in  II- VI  CdZnSe  QW  was  ever  suggested  [Din92, 
Din93,  Din94,  Hag92].  This  excitonic  origin  of  lasing  is  very  different  from  the  conven- 
tional EHP  theory  saying  that  the  optical  gain  in  semiconductor  lasers  is  provided  by  the 
recombination  of  electrons  and  holes  in  the  active  layers  or  quantum  wells  of  our  interests. 

In  this  chapter,  first,  we  will  review  some  of  the  experiments  that  were  claimed  as  the 
evidences  for  the  excitonic  gain  in  CdZnSe  QW  lasers  and  excitonic  emission  in  CdZnSe 
light  emitting  diodes  (LEDs).  Then,  we  will  use  the  models  developed  in  Chapters  2  and 
3  to  show  many-body  effects  (CS,  BGR  and  CE)  also  modify  the  absorption  spectra  at  low 
carrier  densities  as  well  as  gain  spectra  at  high  carrier  densities.  Finally,  we  will  compare 
the  calculation  results  to  the  experimental  observations  and,  without  appealing  to  excitons, 
we  attempt  to  explain  several  spectral  features  claimed  as  proofs  for  excitonic  gain  and 
emission  in  CdZnSe  QW. 

5.2  Claims  of  Evidence  for  Excitonic  Gain 

As  shown  in  Fig.  5.1,  M.  Hagerott  etal.  [Hag92]  observed  close  spectral  correlation 
between  the  so-called  exciton  absorption  and  the  spontaneous  emission  spectra  of  the 
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CdZnSe  QW  light  emitting  diodes  (LEDs)  mT  =  11  and  300  K.  It  was  suggested  that  the 
room  temperature  LED  emission  of  the  CdZnSe  QW  is  of  excitonic  origin.  J.  Ding  et  al. 
[Din92]  reported  the  first  observation  of  so-called  excitonic  gain  and  laser  emission  in 
ZnSe-based  QWs  which  drew  a  great  deal  of  attention  in  the  CdZnSe  research  society  at 
that  time.  Fig.  5.2  shows  their  experimental  data  measured  at  7  =  77  and  295  K  where  laser 
emission  evolved  from  the  low-energy  tail  of  a  resonance-like  absorption  peak.  This  ab- 
sorption peak  was  attributed  to  the  n  =  1  HH  exciton  absorption  resonance  and  the  close 
spectral  following  of  laser  emission  on  its  low-energy  side  was  interpreted  as  excitonic  las- 
ing  in  nature.  A  phenomenological  model  of  optical  gain  was  proposed  based  on  partial 
phase-space  filling  (PSF)  of  an  inhomogeneously  broadened  exciton  resonance  and  includ- 
ing the  effects  of  exciton-phonon  interaction.  This  phenomenological  gain  model  was  used 
to  explain  their  pump-probe  spectroscopic  experiments  at  cryogenic  temperature  up  to 
220K  [Din93].  The  gain/absorption  spectra  of  Fig.  5.3  reconstructed  from  transient  pump- 
probe  techniques  show  the  coexistence  of  optical  gain  with  the  so-called  exciton  absorption 
resonance  in  time  frame  [Din93].  This  was  said  to,  quoted  from  the  paper, "...  provide  the 
ultimate  proof  that  the  well-defined  n  =  \  HH  exciton  resonance  (absorption  peak)  provides 
optical  gain  at  its  low-energy  side."  The  room  temperature  gain/absorption  spectra  of 
CdZnSe  QW  diode  lasers,  which  were  determined  via  a  transformation  from  the  unampli- 
fied  spontaneous  emission  [HenSO,  Kes94]  through  an  opening  in  the  top  contact  of  lasers, 
were  shown  in  Fig.  5.4  and  two  resonance-like  absorption  peaks  appearing  at  low  excitation 
level  were  attributed  to  absoiption  from  n  =  1  HH  and  LH  excitons  [Din94]. 

5.3  Calculation  of  Gain/Absorption  Spectra 

In  Chapter  3,  we  have  shown  how  the  gain  spectra  for  carrier  density  N  greater  than 
the  transparent  carrier  density  N„  change  when  various  many-body  effects  (CS,  BGR  and 
CE)  are  included  one  by  one  in  the  calculations.  As  lowering  N  smaller  than  Af,„  gain  spec- 
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Figure  5.1  Comparison  of  absorption  and  spontaneous  emission  spectra  near  the  so- 
called  n  =  1  HH  exciton  resonance  of  a  CdZnSe  QW  LED  at  T  =  77  K  and 
300  K  [Hag92]. 
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Figure  5.2  Evolution  of  the  edge  emission  of  a  CdZnSe  QW  laser  at  T  =  77  and  295  K. 
The  excitation  level  is  indicated  by  the  percentage  of  what  is  just  above  the 
threshold.  The  so-called  n  =  1  HH  exciton  absorption  resonance  is  shown 
for  reference  and  the  close  spectral  following  of  laser  emission  on  its  low- 
energy  side  was  interpreted  as  excitonic  lasing  in  nature  [Din92]. 
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Figure  5.3  The  so-called  n  =  1  HH  exciton  absorption  resonance  of  a  CdZnSe  QW  laser 
at  7"  =  77  K  without  photoexcitation  (circles),  300  fs  (squares)  and  10  ps  (tri- 
angles) after  the  photoexcitation.  The  laser  emission  spectrum  is  also 
shown.  Coexistence  of  gain  and  resonance-like  absorption  was  elucidated 
as  that  optical  gain  in  CdZnSe  QW  lasers  is  provided  by  excitons  [Din93]. 
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Figure  5.4  Gain/absorption  spectrum  of  a  CdZnSe  QW  diode  laser  at  7  =  300  K,  show- 
ing the  effects  of  increasing  current  injection  on  the  so-called  n  =  1  HH  and 
LH  QW  exciton  states.  The  inset  shows  the  experimental  geometry  for  re- 
cording the  spontaneous  and  stimulated  emission  spectra  [Din94]. 
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tra  become  absorption  spectra.  In  this  section,  we  will  show  how  many-body  effects  mod- 
ify the  absorption  spectra. 

Single-Particle  Model.  The  TE  gain/absorption  spectra  forW  =  2.2  x  10"  and  5.4  x 
10"  cm-3  at  7"  =  300  K,  calculated  using  Eq.  (2.25),  are  shown  in  Fig.  5.5(a).  Two  discon- 
tinuities, resulting  from  the  stair-step  density  of  states  in  QW  structure,  occur  at  photon  en- 
ergies Eg,  ^  and  Eg,^ ,,,  which  coirespond  to  the  energy  gaps  between  CI  and  HH 1  subbands 
and  between  CI  and  LHl  subbands,  respectively.  VorN  =  2.2  x  10'*  cm',  the  whole  spec- 
trum is  under  absorption.  For  N=5.4x  10'*  cm'',  the  spectrum  experiences  gain  at  photon 
energies  between  E^,  ^  and  the  quasi-Fermi  level  separation.  Since  the  measured  gain/ab- 
sorption spectra  shown  in  Fig.  5.4  do  not  have  discontinuous  features,  it  is  clear  that  this 
model  is  incomplete.  Next  we  will  show  when  many-body  effects  are  incorporated  into  the 
single-particle  model,  realistic  gain/absorption  spectra  are  obtained. 

Manv-bodv  effects  -  CS.  The  CS-modified  single-particle  TE  gain/absorption  spec- 
tra are  plotted  in  Fig.  5.5(b)  for  W  =  2.2  x  10'*  and  5.4  x  10"  cm'.  It  is  evident  that  these 
spectra  are  broadened  and  smoothed  relative  to  their  single-particle  counterparts  plotted  in 
Fig.  5.5(a).  The  peak  gain  for  W  =  5.4  x  10'*  cm'  is  reduced  significantly  and  its  spectral 
location  is  blue-shifted  to  higher  photon  energy  (in  this  case  by  about  11 .5  meV).  For  A'  = 
2.2  x  10'*  cm',  the  two  absorption  steps  due  to  Cl-HHl  and  Cl-LHl  transitions  remain 
visible,  but  there  is  no  sign  of  resonance-like  absorption  peaks  as  the  ones  observed  in  ex- 
periment (see  Fig.  5.4). 

Many-body  effects  -  BGR.  The  corresponding  BGR  shifted  spectra  of  Fig.  5.5(b)  are 
shown  in  Fig.  5.5(c)  where  low-energy  tails  of  gain  and  absorption  specti-a  are  well  below 
(~  80  meV)  the  unexcited  Cl-HHl  energy  gap  E^,^  u,.  This  situation  is  often  mistakenly 
elucidated  as  the  evidence  for  existence  of  excitons  in  CdZnSe  QW  lasers  even  at  room 
temperatiu-e  [Din94]. 

Manv-bodv  effects  -  CR.  Plots  of  calculated  TE  gain  spectra  with  CE  included  for 
N  =  2.2x  10'*  and  5.4  x  10'*  cm'  are  shown  in  Fig.  5.5(d).  Comparing  Fig.  5.5(d)  to  (c). 


not  only  does  CE  increase  g  for  high  N,  but  also  enhances  the  absorption  steps  to  resonance- 
like absorption  peaks  for  low  N.  Much  resemblance  between  the  theoretical  results  and  ex- 
perimental observation  will  be  discussed  in  the  following  section  in  details. 

5.4  Comparison  between  Calculations  and  Experiments 

5.4.1  Spectral  Correlation  between  Peaks  of  Absorption  and  Spontaneous  Emission 

As  shown  in  Fig.  5.1,  the  peak  positions  of  the  absorption  peak  and  spontaneous 
emission  peak  were  found  to  be  very  close  to  each  other  at  both  T  =  ll  and  300  K  and  this 
close  spectral  correlation  was  interpreted  as  the  strong  evidence  for  that  the  room  temper- 
ature LED  emission  is  of  excitonic  origin.  Here  we  use  Eq.  (3.63)  (with  CE  included)  to 
calculate  the  spontaneous  einission  spectrum  for /V  =  4.2  x  10"  cm',  as  shown  in  Fig.  5.6. 
The  absorption  spectrum  for  /V  =  2.2  x  10'*  cm"^  in  Fig.  5.4(d)  is  replotted  in  Fig.  5.6  with 
an  inverse  sign  in  arbitrary  units  for  reference.  We  observe  that  the  peak  position  of  the 
spontaneous  emission  spectrum  is  very  close  to  the  position  of  the  resonance-like  absorp- 
tion peak.  It  is  clear  by  comparing  to  Fig.  5. 1  that  we  are  able  to  predict  this  correlation  at 
T=  3(X)  K  without  appealing  to  excitons. 

5.4.2  Energy  Arrangement  of  Some  Key  Spectral  Features 

In  experiment,  the  lasing  emission  spectra  of  CdZnSe  QW  lasers  was  found  to  overlap 
with  the  low-energy  tails  of  their  resonance-like  absorption  spectra,  as  shown  in  Fig.  5.2, 
which  were  often  identified  as  exciton  absorption  spectra  [Din92,  Din93].  The  lasing  en- 
ergy E,  defined  as  the  photon  energy  at  the  peak  of  the  stimulated  emission  spectra,  is  red 
shifted  with  respect  to  E,  and  £,;,  defined  as  the  spectral  position  of  the  «  =  1  HH  exciton 
absorption  peak  and  the  energy  gap  between  the  lowest  subbands  in  CB  and  VB,  on  the  or- 
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Figure  5.5  Calculated  TE  gain/absorption  spectra  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers 
at  carrier  densities  indicated  by  the  numbers  (x  10'*  cm'')  using  (a)  single- 
particle  model;  (b)  with  CS;  (c)  with  CS  and  BGR;  (d)  with  CS,  BGR  and 
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Figure  5.6  Calculated  spontaneous  emission  spectrum  for  the  4  nm  CdojZno.vSe  QW 
LEDs  at  A'  =  4.2  X  10"  cm''  with  CE  included.  The  absorption  spectrum  for 
N  =  2.2x  10"  cm''  in  Fig.  5.4(d)  is  replotted  as  the  dotted  curve  with  an  in- 
verse sign  and  in  arbitrary  units  for  reference. 


125 

der  of  20  and  60  meV  respectively  [Din93]. 

In  theory,  for  typical  cleaved-facet  lasers,  £  can  be  predicted  by  the  spectral  location  of 
the  peak  of  gain  spectrum  as  discussed  in  Section  4.3.  For  example,  E  =  2.427  eV,  deter- 
mined from  the  calculated  gain/absorption  spectrum  forN  =  5.4  x  10"  cm'  in  Fig.  5.5(d), 
corresponds  to  the  lasing  energy  for  a  CdZnSe  QW  laser  with  L  =  1082  nm.  Comparing  to 
the  calculated  spectrum  for  N  =  2.2  x  10"  cm"'  in  Fig.  5.5(d),  E  is  related  to  £,  ci-hhi  = 
2.455  eV  and  Eg,  ^  =  2.482  eV,  defined  as  the  spectral  position  of  the  Cl-HHl  absorption 
peak  and  the  Cl-HHl  subband  energy  gap,  by  a  red  shift  of  about  28  and  55  meV  respec- 
tively. These  numbers  are  comparable  to  the  corresponding  experimentally  measured  val- 
ues mentioned  earlier  (20  and  60  meV  respectively).  In  our  gain  model  based  on  EHP,  the 
red  shifts  are  mainly  due  to  bandgap  renormalization  from  the  screening  of  Coulomb  repul- 
sion among  charge  carriers  and  the  resonance-like  absorption  peaks  are  due  to  the  enhanced 
transition  matrix  element  from  Coulomb  attraction  between  electron-hole  pairs. 

5.4.3  Coexistence  of  Gain  with  Absorption 

The  coexistence  of  optical  gain  with  the  so-called  exciton  absorption  observed  at  low 
temperature  (7"  =  77  K),  as  shown  in  Fig.  5.3,  was  explained  as  the  "ultimate  proof  that  the 
optical  gain  in  CdZnSe  QW  is  provided  by  the  well-defined  n  =  1  HH  exciton  resonance 
[Din93].  For  this  reason,  it  is  interesting  for  us  to  perform  calculation  at  temperatures  lower 
than  300  K  with  CE  included.  As  shown  in  Fig.  5.7(a),  the  TE  gain/absorption  spectrum, 
calculated  at  7  =  200  K  for  W  =  3  x  10'*  cm',  shows  the  existence  of  gain  while  the  Cl- 
HHl  absorption  peak  remains  well  defined  at  higher  energy.  As  decreasing  the  tempera- 
ture further  down  to  160  K,  observed  from  the  calculation  result  plotted  in  Fig.  5.7(b)  for 
A/  =  2.1  X  10"  cm',  the  Cl-HHl  absorption  peak  becomes  sharper  and  narrower  with  co- 
existence of  gain  on  its  low-energy  side.  As  mentioned  in  Section  4.5.2,  the  CE  model  used 
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in  this  work  reaches  it  limits  at  temperature  lower  than  160  K  where  excitons  may  become 
important  Although  we  are  not  able  to  make  predictions  at  77  K  using  this  CE  model,  co- 
existence of  gain  with  resonance-like  absorption  in  CdZnSe  QW  lasers  has  been  predicted 
for  temperature  down  to  160  K,  without  appealing  to  excitons. 

5.4.4  Evolution  of  Gain/Absorption  Spectra  with  Injection  Level 

Fig.  5.8  shows  the  TE  gain/absorption  spectra  at  7  =  300  K  for  various  values  of  N 
which  correspond  to  various  injection  levels.  At  lower  injection  levels  {N  =  2.2  and  3.0  x 
10"  cm"^),  both  Cl-HHl  and  Cl-LHl  absorption  peaks  are  well  defined.  As  increasing  in- 
jection level  (N  =  3.8  and  4.6  x  10"  cm"'),  the  Cl-HHl  absorption  peak  begins  to  saturate. 
For  iV  >  3.8  X  10"  cm"'  gain  appears  on  the  low-energy  side  of  the  spectra.  At  even  higher 
injection  levels  corresponding  to  lasing  threshold  (N  =  5A  and  6.2  x  10"  cm"'),  the  Cl-HHl 
absorption  peak  becomes  hardly  identifiable  while  the  Cl-LHl  absorption  peak  remains  re- 
solvable. These  characteristics  of  the  gradual  evolution  of  gain/absorption  spectra  were 
also  experimentally  observed  at  room  temperature  for  CdZnSe  QW  lasers  [Din94]  (see. 
Fig.  5.4),  but  the  optical  gain  was  said  to  be  excitonic  in  nature. 

5.5  Summary 

As  just  being  briefly  reviewed  in  Section  5.2,  several  key  spectral  featiues  observed 
in  experiments  were  identified  as  the  evidence  for  existence  of  excitons  and  the  optical  gain 
mechanism  in  wide  bandgap  CdZnSe  QW  lasers  was  further  said  to  be  excitonic  in  nature, 
unlike  the  narrow  bandgap  III-V  QW  lasers  for  which  optical  gain  comes  from  EHP 
[Col95,  Coi93].  To  our  best  knowledge,  no  complete  theory  has  ever  been  reported  for  the 
physical  mechanism  responsible  for  absorption  as  well  as  gain  in  II- VI  CdZnSe  QW  diode 
lasers  at  room  temperature.  In  this  chapter,  without  appealing  to  excitons,  we  have  success- 


127 


10000 

0 

-10000 

B      -20000 


3      -30000  E- 


-40000    - 


-50000    - 


-60000 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  [  I  I  I  I 


r=200K 

4  nm  CdZnSe  QW 

w/CS+BGR+CE 


■■■■'■■ 


Cl-HHl 


Cl-LHl 


■  I  ■ 


J_ 


_L 


2.35        2.40        2.45        2.50        2.55        2.60        2.65 
h(0  (eV) 

(a) 


10000 
0 

-10000  F- 

-20000    - 


I      -30000 


-40000 
-50000 


-60000 


1  1  1  1  1  1  ,  1  ,  .  .  .  , 

: 

- 

Cl-HHl            \ 

j 

r  r  =  200K 

\ 

/    r 

:    4nmCdZnSeQW 

\     / 

■ 

:    w/ CS+BGR-i-CE 

Cl-LHl  I  / 

" 

1  1  1  1  1  1  .1  1  1  1  . 

,  1  .  .  ,  .  1  .  .  ,  ,  1 

2.35        2.40        2.45        2.50        2.55 
hca  (eV) 
(b) 


2.60 


2.65 


Figure  5.7  Calculated  TE  gain/absorption  spectra  for  the  4  nm  Cdo  aZn,,  ySe  QW  lasers 
with  CE  at  (a)  r  =  200  K,  W  =  3  X  10'*  cm'^  and  (b)  T=  160  K,  A^  =  2.1  x 
10'*  cm',  showing  the  coexistence  of  optical  gain  and  the  resonance-like 
Cl-HHl  absorption  peak. 
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Figure  5.8  Calculated  TE  gain/absorption  spectra  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers 
with  CE  at  carrier  densities  indicated  by  the  numbers  (x  10'*  cm''),  showing 
the  effects  of  increasing  carrier  density  (or  increasing  injection  level)  on  the 
resonance-like  Cl-HHl  and  Cl-LHl  absorption  peaks. 
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fully  demonstrated  that  these  key  spectral  features  at  room  temperature  can  be  reproduced 
by  using  the  conventional  electron-hole  plasma  theory  modified  by  many-body  effects.  Es- 
pecially, Coulomb  attraction  effect  between  elections  and  holes  in  CdZnSe  QW  lasers  has 
been  shown  being  very  important  at  room  temperature. 


CHAPTER  6 
CARRIER-INDUCED  EFFECTS  IN  II- VI QW  LASERS 


6.1  Introduction 

The  optical  gain  in  semiconductor  lasers  is  provided  by  carrier  injection  into  the  ac- 
tive layer  and  the  sequential  recombinations  of  electron-hole  pairs.  The  relationship  be- 
tween the  optical  gain  and  the  carrier  density  in  the  QW  structure  has  been  described  in 
Chapters  2  and  3  using  the  EHP  theory  with  various  many-body  effects  included.  Since  the 
imaginary  part  of  the  refractive  index  is  linearly  related  to  the  optical  gain  [Col95],  any 
change  in  the  optical  gain  due  to  the  change  in  the  carrier  density  would  resuh  in  a  change 
in  the  imaginary  refractive  index.  In  a  conventional  approach,  the  change  in  the  real  part 
of  the  refractive  index  can  then  be  determined  from  the  change  in  its  imaginary  part  using 
a  Kramers-Kronig  (KK)  transformation  [HenSl].  It  has  been  shown  both  theoretically  and 
experimentally  that  in  ni-V  GaAs-based  bulk  and  QW  lasers,  the  real  refractive  indices  at 
gain  peaks  (or  at  lasing  energies)  reduces  with  increasing  carrier  density  and  thus  with  in- 
creasing optical  gain  [HenSO,  HenSl,  Kes90,  Kes91,  Tie92].  For  the  gain-guided  GaAs- 
based  QW  laser,  the  carrier-induced  reduction  in  the  real  refractive  index  in  the  QW  has 
been  shown  to  affect  the  quality  of  the  laser  beam  [Col91,  Shi89].  However,  on  the  con- 
trary, a  theoretical  calculation  reported  that  the  significant  Coulomb  enhancement  in  room- 
temperature  II- VI  CdZnSe  QW  lasers  results  in  increases  in  the  real  refractive  indices  at 
gain  peaks  with  increasing  carrier  density  [Cho94].  If  this  prediction  is  coirect,  the  beam 
quality  of  the  CdZnSe  QW  lasers  should  be  different  from  that  of  the  GaAS-based  QW  la- 
sers. 
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This  chapter  is  devoted  to  investigating  the  effects  of  injected  carriers  on  the  refrac- 
tive index  and,  therefore,  on  the  beam  quality  of  the  II- VI  CdZnSe  QW  lasers  at  room  tem- 
perature. First,  the  relations  between  the  optical  gain,  the  complex  dielectric  function  and 
the  complex  refractive  index  will  be  derived.  Without  using  the  KK  transformation  method 
[HenSl],  the  real  refractive  index,  as  well  as  the  imaginary  refractive  index,  as  a  function 
of  photon  energy  will  be  calculated  using  the  EHP  theory  with  many-body  effects  included. 
Effects  of  carriers  in  the  QW  on  the  refractive  index  will  be  presented.  A  qualitative  de- 
scription of  the  effect  of  canier-induced  refractive  index  change  on  the  lateral  profile  of  a 
laser  beam  at  the  laser  facet  (near-field  pattern)  will  be  given.  Finally,  experimental  results 
from  the  measurement  of  the  near-field  patterns  using  CdZnSe  QW  lasers  will  be  present- 
ed. 

6.2  Effect  of  Injected  Carriers  on  the  Refractive  Index 

6.2.1  Complex  Optical  Dielectric  Function  ^' 

The  complex  optical  dielectric  function  e(fico)  is  related  to  the  complex  optical  sus- 
ceptibility x(fto))  by  [Cho94,  Hau90,  Kli97] 

e(fia))  =  e^[l+z(Sa))]  ,  (6.1) 

where  e„  is  the  background  dielectric  constant.  In  Eq.  (3.57),  the  relation  between  the  op- 
tical gain  function  g  (hoi)  and  the  imaginary  part  of  %  {h(o)  is  given  and  in  an  analogy,  one 
can  define  a  complex  function  g  (hw)  as 

giha)  =  ^ixinm),  (6.2) 

where  the  real  part  of  g  (h(o)  is  die  optical  power  gain  as  the  one  given  in  Eq.  (3.57)  and  the 
imaginary  part  is  the  carrier-induced  phase  shift  [Cho94].  Using  Eq.  (6.2)  in  Eq.  (6. 1),  we 
obtain        ... 


132 

E(fico)  =  e„ri-i-^g(ftto)l.  (6.3) 

L        con  J 

Equating  the  real  and  imaginary  parts  of  the  left-hand-side  (LHS)  to  those  of  the  right-hand- 
side  (RHS)  with  the  use  of  the  following  definitions  for  e  (to)  and  g  (hm), 

e(Sco)  =e^(fico) +ie;(Sco)  ,  (6.4) 

g  ihw)  =  g  ihw)  +  ig-  ihuy)  ,  (6.5) 

we  have 

e^(ftco)  =e„  +  ^«,.(fia)) 

(6.6) 
e;(fia))  =  —^(^ico) 

where  the  relation  e^  =  «^  has  been  used.  Eq.  (6.6)  shows  that  the  imaginary  part  of 
g  (fico)  gives  rise  to  the  real  part  of  e  {tun)  and  the  real  part  of  g  (tuo)  gives  rise  to  the  imag- 
inary part  of  e  (fia). 

6.2.2  Complex  Refractive  Index 

Since  the  square  of  the  complex  refractive  index  is  equal  to  the  complex  dielectric 
function  as  shown 

P{h(a)  =  I  (SCO)  ,  (6.7) 

where 

n  (ftco)  =  n^  (fico)  +  in-  (ha)  .  (6.8) 

Using  Eqs.  (6.8)  and  (6.4)  in  the  LHS  and  RHS  of  Eq.  (6.7),  respectively,  and  equating  the 
real  and  imaginary  parts  on  both  sides,  respectively,  we  have 

n}(h(ii)  -nf(tiu>)  =  e,  (Sco) 

(6.9) 
2n^  (Sco) /),.  (Sco)  =  e;(ftco) 

Solving  the  dual  equations  in  Eq.  (6.9)  for  the  real  refractive  index  n  (fuo)  in  terms  of 
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e^  (tuo)  and  e^  (ftco),  we  obtain 


n^iha)  =  Je,(ft03)  +\jt'i{hm)  +e?(^M)  .  (6.10) 

For  the  imaginary  refractive  index  n,.  (ftco),  we  can  also  solve  it  from  Eq.  (6.9)  and  obtain 


«;(^(o)  =  +J-e,(^a))  +ye}(fi(a)  +e?(ficD) ,  (6.11) 

or  we  can  use  the  basic  definition  that  relates  the  power  gain  to  the  imaginary  refractive 
index  as  [Col95] 

^.-V.Cio,)   =-^«,.(ftco).  (6.12) 

where  k„  is  the  propagation  constant  in  vacuum.  As  a  consequence,  the  choice  of  positive 
or  negative  root  in  Eq.  (6.1 1)  for  n-  (fiw)  depends  on  the  sign  of  g  {fua)  due  to  the  fact  that 
they  are  opposite  in  sign  as  shown  in  Eq.  (6. 12). 


6.2.3  Simulation  results 

The  procedure  for  calculations  of  the  real  and  imaginary  refractive  indices  as  a  func- 
tion of  photon  energy,  which  are  n,  (ha)  and  n,-  (tua),  respectively,  is  summarized  as  fol- 
lows: (1)  Calculate  the  complex  function  g(hoi)  by  dropping  the  notations  Real  in  Eqs. 
(3.1)  and  (3.63)  for  models  with  and  without  CE,  respectively;  (2)  From  g  (ha),  calculate 
die  complex  dielectric  function  using  Eq.  (6.3);  (3)  Calculate  n,  (ha)  using  Eq.  (6.10)  and 
n,(»(0)  using  Eq.  (6.1 1)  or  Eq.  (6.12). 

Manv-bodv  effects  -  CS.  As  shown  in  Fig.  4.1 3(a),  gain  spectra  were  calculated  with 
inclusion  of  CS  for  the  4  nm  Cdo.3Zno.7Se  QW  lasers  at  N  =  6.7,  7.4,  8.3,  9.2  x  10'*  cm"' 
which  are  the  carrier  densities  required  to  reach  the  laser  action.  At  these  values  of  A',  the 
spectra  of  the  real  and  imaginary  refractive  indices  are  calculated  with  CS  and  shown  in 
Fig.  6. 1.  The  real  refractive  indices  around  where  the  gain  peaks  locate,  as  indicated  by  the 
open  circles  in  Fig.  6.1(a),  are  found  to  decrease  when  N  increases.   The  spectra  of  the 
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imaginary  refractive  index  in  Fig.  6.1(b)  and  the  gain  spectra  in  Fig.  4.13(a)  are  very  alike 
with  opposite  signs  to  each  other,  as  expected  from  Eq.  (6.12). 

Manv-bodv  effects  -  BGR.  With  the  shift  of  BGR  shown  in  Fig.  3.6,  the  correspond- 
ing rigidly-shifted  spectra  in  Figs.  6.1(a)  and  (b)  are  shown  in  Fig.  6.2(a)  and  (b),  respec- 
tively. Although  the  spectral  location  of  peak  gain  shows  a  red  shift  with  increasing  N,  the 
real  refractive  indices  around  gain  peaks  remain  decreasing  as  N  increases. 

Many-bodv  effects  -  CE.  With  inclusion  of  CE,  in  addition  to  CS  and  BGR,  Figs. 
6.3(a)  and  (b)  show  the  spectra  of  n,  (ftco)  and  n,  (fico),  respectively,  calculated  at  W  =  5.0, 
5.4, 5.7, 6.0  X  10"  cm''  which  are  the  required  N  values  to  reach  laser  action  in  this  model, 
as  discussed  in  section  4.3.2  and  shown  in  Fig.  4.15.  As  comparing  to  Figs.  6.1(a)  and 
6.2(a),  the  real  refractive  indices  in  Fig.  6.3(a)  are  a  bit  higher  at  gain  peaks  and  sweep  over 
a  wider  range  around  the  gain  peaks.  A  reduction  in  the  real  refractive  index  at  gain  peak 
is  also  found  with  increasing  N. 

6.3  Beam  Quality 

6.3.1  Antiguiding  Factor 

In  the  former  section,  the  changes  in  the  real  and  imaginary  refractive  indices  due  to 
the  injected  carriers  are  shown.  Since  the  refractive  index  is  important  in  determining  the 
mode  profile  in  a  laser  structure,  the  change  in  the  refractive  index  will  affect  the  quality  of 
a  laser  beam.  Starting  with  a  relationship  between  how  the  real  and  imaginary  indices  are 
affected  by  the  carrier  density,  it  is  described  using  what  is  referred  to  as  the  antiguiding 
factor  ;  ,      C'  I '. 

■,.    .''•    .  ('  ■■  ■  -.    ;•-'       dn^(ha3) 
dN 
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Figure  6. 1  Calculated  spectra  of  the  (a)  real  and  (b)  imaginary  refractive  indices  for  the 
4  nm  Cdo.3Zno7Se  QW  lasers  with  CS  included.  The  open  circles  indicate 
the  values  of  the  real  and/or  imaginary  indices  at  gain  peaks  for  the  carrier 
densities  indicated  by  the  numbers  (x  10'*  cm"^). 
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Figure  6.2        Calculated  spectra  of  the  (a)  real  and  (b)  imaginary  refractive  indices  for  the 
4  nm  Cdo.3Zno,7Se  QW  lasers  with  CS  and  BGR  included.  The  open  circles 
indicate  the  values  of  the  real  and/or  imaginary  indices  at  gain  peaks  for  the 
*      carrier  densities  indicated  by  the  numbers  (x  10'*  cm"'). 
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Figure  6.3  Calculated  spectra  of  the  (a)  real  and  (b)  imaginary  refractive  indices  for  the 
4  nm  Cdo.3Zno.7Se  QW  lasers  with  CS,  BGR  and  CE  included.  The  open  cir- 
cles indicate  the  values  of  the  real  and/or  imaginary  indices  at  gain  peaks  for 
the  carrier  densities  indicated  by  the  numbers  (x  10'*  cm"^). 
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Substituting  n,  (fica)  with  Eq.  (6. 1 2),  Eq.  (6. 1 3)  can  be  rewritten  as 

dn^ihw)  dn^ihw) 

a(fito)   = ,    ,.    ,      =  -^ — ,   ,.    ,     ■  (0.14) 

,  c     dgihoi)  X     dg  (tia) 

.   -      -  '  dN  (/N 

where  the  changes  in  optical  gain  and  real  refractive  index  are  now  combined  to  determine 
the  antiguiding  factor  [Cho94].  The  paraineter  a  (Sco)  is  often  also  referred  to  as  the  Une- 
width  enhancement  factor  [Kes90,  Kes91,  Tie92],  since  it  is  a  key  quantity  in  determining 
the  linewidth  and  chirp  of  single  mode  lasers  to  be  used  for  optical  communications.  This 
topic  is  beyond  the  scope  of  present  work  and  left  to  the  discussion  in  [Osi87].  The  focus 
of  this  chapter  is  placed  on  how  the  antiguiding  factor  is  related  to  the  beam  quality. 

Usually,  the  term  dg  (ti(a)/dN  in  Eq.  (6.14)  is  positive  since  g  (tiO))  increases  with  N 
as  shown  in  Figs.  4.10(a),  4.1 1(a)  and  4.12(a)  for  the  InGaAs  QW  lasers  and  Figs.  4.13(a), 
4.14(a)  and  4.15(a)  for  the  CdZnSe  QW  lasers.  If  the  real  refractive  index  decreases  with 
increasing  N,  then  dn^  (h(i>)/dN  is  negative  which  results  in  a  positive  value  for  a  (ftco). 
From  the  spectra  of  n,  (ha)  shown  in  Figs.  6. 1  (a),  6.2(a)  and  6.3(a),  values  of  n,  (Sw)  around 
gain  peaks  decrease  with  increasing  N  in  the  QW.  Asa  consequence,  the  spectra  of  a  (Sco) 
shown  in  Figs.  6.4(a),  (b)  and  (c)  for  all  three  models  have  positive  antiguiding  factors  at 
where  gain  exists.  The  calculation  with  CE  included  still  shows  positive  antiguiding  fac- 
tors for  the  CdZnSe  QW  lasers  like  what  have  been  observed  and  predicted  in  the  GaAs- 
based  lasers  [HenSO,  HenSl,  Kes90,  Kes91].  However,  this  calculation  result  is  opposite 
to  the  previous  prediction  where  the  strong  Coulomb  enhancement  in  the  CdZnSe  QW  la- 
sers results  in  negative  values  for  the  antiguiding  factor  around  gain  peaks  [Cho95b].  A 
negative  antiguiding  factor  means  that  the  real  refractive  index  increases  with  the  carrier 
density  according  Eq.  (6.14). 


139 


I 


3  i-T 


0 


I    '    I    '  \V\ '    I    '    I    '    I    '    I    '    I    '    I 


r=300K 
L  4  nm  CdZnSe  QW 
:    w/CS 

I    ■    I    I    I    I    I 


.    I    .    r    ■    I 


2.45  2.47  2.49  2.51 

ftco  (eV) 
(a) 


2.53 


2.55 


r=300K 
L  4nmCdZnSeQW 
w/CS+BGR 
I    ■    I    ■    I 


J_ 


O- 


J_ 


2.35 


2.37 


2.39 


2.41 


2.43 


SCO  (eV) 
(b) 


2.45 


Figure  6.4  Calculated  spectra  of  the  antiguiding  factor  for  the  4  nm  Cdo.sZnojSe  QW 
lasers  (a)  with  CS,  (b)  with  CS  and  BGR,  and  (c)  with  CS,  BGR  and  CE  in- 
cluded. The  open  circles  indicate  the  values  of  antiguiding  factors  at  gain 
peaks  for  the  carrier  densities  indicated  by  the  nutnbers  (x  10"  cm'^). 
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6.3.2  A  Qualitative  Description 

In  this  section,  a  qualitative  description  will  be  given  on  how  the  antiguiding  factor, 
which  is  associated  with  the  carrier-induced  refractive  index  change  discussed  earlier,  has 
important  implications  for  the  laser  beam  quality.  For  mathematical  models  on  this  subject, 
we  leave  them  to  the  literatures  [Coo75,  Kir77,  Nas73]. 

Fig.  6.5(a)  is  a  simplified  cross-section  sketch  of  a  stripe-geometry  gain-guided  QW 
laser  with  a  stripe  width  of  w,  where  j:-direction  is  the  lateral  direction,  y-direction  is  along 
the  axis  of  the  stripe  and  z-direction  is  the  growth  direction  (normal  to  the  plane  of  QW). 
In  the  QW  layer,  the  region  beneath  the  stripe  is  defined  here  as  the  "pumped  section"  and 
the  "unpumped  section"  is  referred  to  as  the  region  on  the  either  side  of  the  pumped  section. 
Assuming  that  current  spreading  can  be  neglected,  injection  current  is  then  confined  to  the 
pumped  section  and  the  lateral  current  profile  is  sketched  in  Fig.  6.5(b).  Below  threshold, 
the  injection  current  creates  a  lateral  carrier  density  profile,  as  shown  in  Fig.  6.5(c),  where 
the  injected  carrier  density  N  <  N,^  only  exists  in  the  pumped  section  when  lateral  carrier 
diffusion  is  ignored.  The  injected  carriers  go  into  spontaneous  emission  since  the  stimulat- 
ed emission  is  negligible  below  threshold.  If  the  antiguiding  factor  is  positive,  the  real  re- 
fractive index  will  have  an  inverse  profile  of  the  carrier  density  profile.  As  a  consequence, 
the  real  refractive  index  is  lower  beneath  the  stripe  than  the  edges  producing  a  so-called  car- 
rier-induced antiguiding  effect,  as  shown  in  Fig.  6.5(d).  This  antiguiding  effect  has  been 
observed  leading  to  worse  lateral  beam  confinements  and  higher  threshold  currents  in  nar- 
row-stripe gain-guided  InGaAs  QW  lasers  [Bee91,  Shi89]  than  in  their  ridge-guided  coun- 
terparts [Off89].  Fig.  6.5(e)  shows  that  if  the  antiguiding  factor  is  negative,  which  is  very 
unusual,  the  real  refractive  index  in  the  pumped  section  is  higher  than  those  in  the  un- 
pumped sections  on  both  sides,  resulting  in  a  better  lateral  beam  confinement  for  narrow- 
stripe  lasers. 

Above  threshold,  that  is  /  >  /,;,,  a  narrow-stripe  laser  (w  <  10  nm)  typically  lases  in  a 
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fundamental  lateral  mode  where  the  optical  intensity  distribution  (lateral  near-field  pattern) 
of  the  stimulated  emission  fills  the  entire  width  of  the  pumped  section,  as  illustrated  in  Fig. 
6.6(a).  For  the  piupose  of  this  qualitative  explanation,  the  cosine-form  of  intensity  profile 
is  replaced  by  a  rectangular  profile.  Since  all  extra  carriers  above  N,;,  go  into  stimulated 
emission,  carrier  density  in  the  pumped  section  is  clamped  at  N^,  as  shown  in  Fig.  6.6(b). 
Again,  if  the  antiguiding  factor  is  positive,  the  real  refractive  index  profile  in  Fig.  6.6(c) 
will  have  a  lower  value  in  the  pumped  section  than  in  the  either  side  and  the  optical  field 
will  experience  a  carrier-induced  antiguiding  effect,  in  addition  to  the  gain-guiding  mech- 
anism provided  by  the  gain  generated  along  the  axis  of  the  stripe.  Much  of  the  light  will  be 
absorbed  in  the  lossy  unpumped  sections.  If  the  antiguiding  factor  is  negative,  the  beam 
will  have  a  better  lateral  confinement  due  to  the  convex  real  index  profile  in  Fig.  6.6(d). 

For  wide-stripe  gain-guided  lasers  {w  >  20  \xm),  as  the  current  increases  above  thresh- 
old, they  usually  will  lase  in  higher  order  lateral  mode,  which  has  more  than  one  peak  in 
the  lateral  optical  intensity  profile  as  illustrated  in  Fig.  6.7(a).  Each  peak  corresponding  a 
lasing  filament  and  the  total  width  of  lasing  filaments  are  less  than  the  stipe  width.  The 
carrier  density  is  clamped  at  N,/,  inside  a  lasing  filament  and  continues  to  rise  above  N,,,  on 
the  either  side  of  the  filament.  This  phenomenon  is  referred  to  as  spatial  hole  burning  and 
illustrated  in  Fig.  6.7(b)  based  on  the  assumptions  of  the  rectangular  intensity  profile  and 
no  lateral  carrier  diffusion.  In  the  case  of  the  positive  antiguiding  factor,  as  shown  in  Fig. 
6.7(c),  the  real  refractive  index  variation  is  inverse  to  the  carrier  density  variation.  There- 
fore, the  real  index  inside  the  filament  is  higher  than  that  in  the  either  side  of  the  filament, 
which  forms  a  self-focusing  carrier-induced  index  waveguide.  As  a  consequence,  the  size 
of  the  self-focused  filament  will  be  very  small  (<  5  nm)  [Kir77].  If  the  antiguiding  factor 
is  negative,  the  real  refractive  index  variation  follows  the  variation  of  the  carrier  density,  as 
shown  in  Fig.  6.7(d),  which  produces  an  antiguiding  effect  for  the  filaments  and  suppresses 
filamentary  lasing  to  achieve  a  better  beam  quality. 
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Figure  6.5  Below  threshold,  effect  of  carrier-induced  real  refractive  index  change  illus- 
trated by  (a)  cross-section  sketch  of  a  gain-guided  stripe-geometry  QW  laser 
and  sketched  lateral  profiles  of  (b)  current;  (c)  carrier-density;  (d)  real  re- 
fractive index  when  the  antiguiding  factor  a  >  0;  (e)  real  refractive  index 
when  a  <  0. 
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Figure  6.6  Above  threshold,  effect  of  carrier-induced  real  refractive  index  change  in 
narrow-stripe  gain-guided  QW  lasers  illustrated  by  the  sketched  lateral  pro- 
files of  (a)  optical  intensity;  (b)  carrier  density;  (c)  real  refractive  index 
when  the  antiguiding  factor  a  >  0;  (d)  real  refractive  index  when  a  <  0. 
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Figure  6.7  Above  tlireshold,  effect  of  carrier-induced  real  refractive  index  change  in 
wide-stripe  gain-guided  QW  lasers  illustrated  by  the  sketched  lateral  pro- 
files of  (a)  optical  intensity;  (b)  carrier  density;  (c)  real  refractive  index 
when  the  antiguiding  factor  a  >  0;  (d)  real  refractive  index  when  a  <  0. 
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6.4  Near-Field  Measurement 

One  important  measurement  indicating  the  quality  of  a  laser  beam  is  measuring  its 
optical  intensity  distribution  at  the  laser  facet  (sketched  as  an  oval  shape  in  Fig.  4.19), 
which  is  usually  referred  to  as  the  near-field  measurement.  In  this  work,  near-field  patterns 
at  various  injection  levels  are  measured  using  two  configuration  types  of  CdZnSe  QW  la- 
sers fabricated  from  the  same  wafer  with  an  epitaxial  structure  as  sketched  in  Fig.  2. 2.  One 
type  is  the  narrow-stripe  buried-ridge  index-guided  lasers  with  a  stripe  width  w  =  5  |J,m  and 
a  cavity  length  L  =  1219  nm,  as  depicted  in  Fig.  4.2,  and  the  other  type  is  the  wide-stripe 
gain-guided  lasers  with  w  =  20  |J.m  and  L  =  1300  nm,  as  shown  in  Fig.  6.8  with  a  sketch  of 
a  near-field  pattern.  Both  types  of  lasers  were  characterized  in  bar  forms  using  current  puls- 
es at  room  temperature.  Fig.  6.9  illustrates  the  experimental  setup  of  the  near-field  mea- 
surement. The  near-field  pattern  at  the  laser  facet  was  imaged  onto  the  paper  screen  using 
a  SOX  microscope  objective.  The  image  on  the  paper  screen  was  captured  by  the  Hamamat- 
su  C2400  CCD  camera  and  monitored  on  a  TV  screen.  The  image  of  the  measured  near- 
field  pattern  was  also  sent  to  a  Pentium  120  MHz  PC  equipped  with  a  viewgraber  interface 
for  image  processing  and  data  analysis  using  the  software  HLImage**. 

Fig.  6.10  shows  plots  of  measured  lateral  near-field  intensity  profiles  for  the  20  |im 
wide-stripe  gain-guided  CdZnSe  QW  laser  with  f,^  =  150  mA.  At  /  =  250  mA,  several 
bright  spots  are  found  along  the  lateral  direction  in  the  near-field  pattern  which  clearly  show 
up  as  several  sharp  peaks  in  the  analyzed  lateral  near-field  intensity  profile.  This  filamen- 
tary lasing  phenomenon  in  the  wide-stripe  gain-guided  lasers  indicates  that,  as  discussed  in 
Section  6.3.2  and  illustrated  in  Fig.  6.7,  CdZnSe  QW  lasers  have  positive  antiguiding  fac- 
tors at  lasing  energies.  This  agrees  with  the  calculation  results  shown  in  Fig.  6.4  and  thus 
invalidates  a  previous  prediction  that  CdZnSe  QW  lasers  should  have  negative  antiguiding 
factors  at  where  optical  gain  exists  due  the  strong  Coulomb  attraction  effect  [Cho94].  In  a 
comparison  with  the  wide-stripe  gain-guided  lasers,  near-field  patterns  for  the  5  nm  nar- 
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Figure  6.8        Cross-section  sketch  of  the  20  |im  wide  gain-guided  CaZnSe/ZnSSe  SCH 
single  QW  laser  devices  provided  by  the  3M  Company. 
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Figure  6.9  Illustration  of  experimental  setup  for  measuring  the  near-field  pattern  of 
CdZnSe  QW  diode  lasers.  The  near-field  pattern  at  the  laser  facet  was  im- 
aged to  a  paper  screen  using  a  microscope  objective.  The  image  was  cap- 
tured by  a  CCD  camera  and  displayed  on  a  TV  monitor.  The  image  was  then 
analyzed  using  the  software  HLImage**  in  a  computer. 
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row-stripe  buried-ridge  index-guided  CdZnSe  QW  laser  with  /^  =  25  mA  were  also  mea- 
sured at  various  injection  currents.  The  analyzed  lateral  near-field  intensity  profiles,  as 
shown  in  Fig.  6.11(b),  show  smooth  single-lobe  fundamental  mode  profiles  with  no  sign  of 
filamentary  lasing  due  to  the  built-in  lateral  index  waveguide  of  the  buried-ridge  structure. 

6.5  Summary 

Complex  refractive  index  as  a  function  of  photon  energy  was  derived  from  the  com- 
plex optical  susceptibility  using  the  EHP  theory  with  various  many-body  effects  included. 
Spectra  of  real  and  imaginary  refractive  indices  were  calculated  for  the  CdZnSe  QW  struc- 
ture at  various  carrier  densities.  It  was  found  that  the  real  refracuve  indices  around  gain 
peaks  decrease  with  increasing  carrier  density.  As  a  consequence,  the  associated  antiguid- 
ing  factors  were  found  having  positive  values  at  gain  peaks  for  the  CdZnSe  QW  structure, 
even  when  the  significant  Coulomb  enhancement  was  taken  into  account  in  the  calculadon. 
This  result  was  opposite  to  a  previous  work  where  negative  antiguiding  factors  at  gain 
peaks  were  predicted  for  the  CdZnSe  QW  structure  [Cho94].  The  filamentary  lasing  phe- 
nomenon observed  in  the  measured  near-field  pattern  for  the  wide-stripe  gain-guided 
CdZnSe  QW  laser  indicated  agreement  with  our  prediction  that  CdZnSe  QW  lasers  should 
have  positive  antiguiding  factors  at  lasing  energies. 
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Figure  6.10  Analyzed  lateral  intensity  profiles  of  the  measured  near-field  patterns  for 
the  20  |xm  wide  gain-guided  CdZnSe  QW  lasers  with  I,,,  =  150  mA  at/  =  50 
and  250  mA. 
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Figure  6.11  Analyzed  lateral  intensity  profiles  of  the  measured  near-field  patterns  for 
the  5  Jim  wide  buried-ridge  index-guided  CdZnSe  QW  lasers  with  I,^  =  25 
mA  at  /  =  10, 20,  30  and  40  mA. 


CHAPTER? 
CONCLUSIONS  AND  FUTURE  WORKS 


7.1  Conclusions 


Effects  of  Coulomb  attraction  between  injected  electrons  and  holes  in  semiconductor 
quantum  well  lasers  were  studied.  The  investigation  was  based  on  a  realistic  model  pre- 
sented in  Chapter  2  and  3,  which  was  developed  from  the  conventional  electron-hole  plas- 
ma theory  [Col95,  Cor93],  modified  by  many-body  effects  [Asa93,  Cho94,  Hau90], 
generalized  for  the  muW-subband  system  and  incorporated  with  the  effective-mass  anisot- 
ropy  due  to  the  biaxial  strain  in  the  quantum  well.  Numerical  calculations  showed  that 
Coulomb  attraction  enhances  the  optical  gain  and  spontaneous  emission  in  quantum  well 
laser  structures.  This  so-called  Coulomb  enhancement  was  found  to  be  stronger  in  blue- 
green  CdZnSe  QW  lasers  than  in  infrared  InGaAs  QW  lasers. 

With  inclusion  of  Coulomb  enhancement,  the  spectral  location  of  peak  gain  was 
found  to  shift  to  higher  energy  with  increasing  injected  carrier  density  in  either  InGaAs  or 
CdZnSe  QW  laser  structures.  The  importance  of  Coulomb  enhancement  on  predicting  the 
carrier-induced  energy  shift  in  peak  gain  location  was  justified  by  good  agreement  found 
between  the  measured  (see  Figs.  4.5  and  4.8)  and  predicted  (see  Figs.  4.12  and  4.15)  in- 
creases in  lasing  energy  with  decreasing  cavity  length  (increasing  carrier  density).  For  In- 
GaAs QW  lasers,  the  inclusion  of  Coulomb  enhancement  in  the  prediction  [Hsu97a]  has 
solved  a  problem  raised  in  previous  studies  on  GaAs-based  QW  lasers  [Che93,  Chi88]  that 
the  ignorance  of  Coulomb  enhancement  had  yielded  essentially  a  constant  spectral  location 
of  peak  gain  with  increasing  carrier  density.  This  is  due  to  the  fact  that  the  bandgap  renor- 
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malization  red  shift  and  band  filling  blue  shift  offset  each  other  in  GaAs-based  QW  lasers 
(see  Fig.  4.6).  For  CaZnSe  QW  lasers,  our  study  [Hsu97b]  has  proved  a  previous  predic- 
tion to  be  incorrect  [Cho95b]  which  showed  a  red  shift  in  the  peak  gain  energy  with  increas- 
ing carrier  density. 

In  studying  the  temperature  dependence  of  lasing  energy  using  CdZnSe  QW  lasers, 
good  agreement  was  found  between  the  prediction  with  Coulomb  enhancement  and  exper- 
imental data  at  temperature  ranging  from  160  K  to  300  K  (see  Fig.  4.30).  This  has  demon- 
strated that  electron-hole  plasma  is  responsible  for  lasing  action  in  CdZnSe  QW  lasers  at 
room  temperature  down  to  160  K.  Below  160  K,  excitons  may  play  a  more  important  role 
in  lasing  mechanism  for  CdZnSe  QW  lasers. 

Gain  versus  radiative  current  density  {g-J,J)  relationship  was  calculated  for  both  In- 
GaAs  and  CdZnSe  QW  lasers  (see  Fig.  4.21  and  4.22).  With  inclusion  of  Coulomb  en- 
hancement, the  calculated  g-J,ai  curves  were  found  to  have  higher  transparencies  and  high- 
er slopes  than  the  ones  predicted  without  Coulomb  enhancement  included.  Using  these  g- 
J,^  curves,  the  predicted  threshold  currents  for  both  InGaAs  and  CdZnSe  QW  lasers  with 
Coulomb  enhancement  included  (see  Figs.  4.25  and  4.26)  were  found  to  give  good  agree- 
ment with  the  measured  threshold  currents  at  various  cavity  lengths  (see  Figs.  4.17  and 
4.18).  The  predictions  without  Coulomb  enhancement  were  found  to  overestimate  the 
threshold  currents.  This  study  has  corrected  a  previous  work  on  the  calculation  of  gain-cur- 
rent characteristics  in  CdZnSe  QW  structures  including  many  body  effects  [Ree95a]  where 
the  predicted  threshold  currents  were  much  too  high  comparing  to  our  measured  data. 

For  CdZnSe  QW  laser  structure.  Coulomb  attraction  effects  were  shown  to  not  only 
increase  the  peak  gain  but  also  enhance  the  absorption  steps  to  resonance-like  absorption 
peaks  (see  Fig.  5.5).  The  experimentally  observed  absorption  peaks  in  CdZnSe  QW  lasers 
at  room  temperature  (see  Fig.  5.4),  which  were  mistakenly  attributed  to  exciton  absorption 
resonances  [Din94],  were  thus  successfully  predicted  for  the  first  time  by  using  the  elec- 
tron-hole plasma  theory  with  Coulomb  attraction  and  other  many-body  effects  taken  into 
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account 

Spectra  of  complex  refractive  index  were  calculated  including  many-body  effects  for 
the  CdZnSe  QW  laser  structure  at  various  carrier  densities  (see  Figs.  6. 1, 6.2  and  6.3).  The 
real  refractive  indices  around  gain  peaks  were  found  to  reduce  with  increasing  carrier  den- 
sity in  the  QW,  even  when  Coulomb  enhancement  was  included.  The  resulting  positive  an- 
tiguiding  factors  (see  Fig.  6.4)  disagreed  with  the  negative  antiguiding  factors  predicted  for 
CdZnSe  QW  lasers  in  a  previous  work  [Cho95b].  However,  the  filamentary  lasing  phe- 
nomenon observed  in  the  20  ixm  wide  gain-guided  CdZnSe  QW  laser  (see  Fig.  6.10)  im- 
plied tiiat  CdZnSe  QW  lasers  should  have  positive  antiguiding  factors  at  gain  peaks. 

In  conclusion,  Coulomb  attraction  effects  are  very  important  in  semiconductor  quan- 
tum well  lasers,  especially  in  wide-band  gap  laser  materials  with  small  dielecQic  constants 
such  as  CdZnSe  QW  lasers.  At  room  temperature,  recombinations  of  electron-hole  plasma 
are  responsible  for  lasing  action  in  CdZnSe  QW  lasers,  not  excitons.  The  theory  with  Cou- 
lomb enhancement  works  well  in  CdZnSe  QW  lasers  down  to  160  K.  As  reducing  the  tem- 
perature below  160  K,  excitons  may  become  more  important 

7.2  Future  Works 

Several  related  future  works  are  suggested  as  follows: 

1.  Theoretical  models  developed  in  this  work  can  be  applied  to  other  material  systems 
and  utilized  for  optimization  of  laser  devices.  Recently  there  has  been  a  considerable  ex- 
pansion of  research  activities  on  Hl-V  GaN-based  laser  materials.  This  expansion  is  attrib- 
uted to  the  commercialization  of  high-brightness  ni-V  GaN-based  blue-green  LEDs  by 
Nichia  Chemical  Industries  and  the  first  demonstration  of  blue  lasing  emission  from  laser 
diodes  fabricated  from  this  material  system  [Nak96].  Several  theoretical  studies  based  on 
the  theory  including  many-body  effects  have  shown  that  Coulomb  attraction  effects  are 
strong  in  these  GaN-based  laser  structures  [Cho95c,  Ree95b].  However,  many  unrealistic 
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assumptions  were  made  in  these  works  due  to  lack  of  many  parameters  for  the  hexagonal 
GaN-based  materials.  In  addition,  similar  works  done  for  the  zinc-blend  CdZnSe  material 
system  by  the  same  research  groups  [Cho95b,  Ree95a]  were  shown  in  this  dissertation  to 
conflict  with  experimental  measurements.  It  will  be  interesting  to  simulate  GaN-based  QW 
lasers  using  our  Coulomb  enhancement  model  developed  in  this  work,  once  the  parameters 
for  GaN-based  materials  become  available  in  the  future. 

2.  The  Coulomb  enhancement  model  developed  in  this  dissertation  reaches  its  limit 
for  carrier  density  much  lower  than  the  Mott  density.  This  short  coming  is  mainly  due  to 
the  fact  that  the  Padfi  approximation,  made  in  deriving  the  Coulomb  enhancement  factor, 
is  an  approximation  for  carrier  densities  higher  than  the  Mott  density.  If  this  limitation  can 
be  removed,  then  the  absorption  and  refractive  index  spectra  in  the  unpumped  section  of  a 
laser  device  can  be  calculated.  This  may  be  beneficial  to  the  understanding  and  designing 
of  laser  devices  which  utilize  index  modulation  between  the  pumped  and  unpumped  sec- 
tions such  as  antiguiding  array  laser  devices. 

3.  In  the  recent  years,  a  new  type  of  semiconductor  lasers,  which  utilizes  transitions 
between  subbands  in  the  conduction  band  quantum  wells,  is  under  fast  development  [Fai94, 
Fai95,  Hsu97c].  In  such  so-called  intersubband  laser  devices,  photons  are  generated  from 
electrons  making  transitions  from  a  higher  subband  to  a  lower  subband  in  the  conduction 
band,  as  shown  in  Fig.  7.1,  which  is  unlike  the  conventional  interband  QW  lasers  discussed 
in  this  work  that  photons  are  generated  from  recombinations  of  conduction  band  electrons 
and  valence  band  holes.  Since  electrons  are  the  only  type  of  carriers  participating  in  the 
photon  emission  process  in  intersubband  lasers.  Coulomb  attraction  effects  are  not  expect- 
ed to  be  important.  But  other  many-body  effects,  such  as  carrier  scattering  with  phonons 
and  carrier  screening,  still  need  to  be  investigated  to  improve  the  design  of  intersubband 
lasers. 
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Figure  7. 1  Conduction  band  energy  diagram  of  a  portion  of  the  waveguide  core  region 
of  a  vertical  transition  quantum  cascade  laser  [Fai95]  under  a  threshold  for- 
ward bias  electric  field. 


'•-'•■.  APPENDIX  A 

MATERIAL  PARAMETERS  II- VI  TERNARIES 

A.l  Interpolation  Method 

The  interpolation  scheme  is  known  to  be  a  useful  tool  for  estimating  some  physical 
parameters  of  semiconductor  compounds.  If  one  uses  the  linear  interpolation  scheme,  the 
material  parameter  P  of  the  ternary  alloy  Afiifi  can  be  derived  from  those  of  its  relevant 
binaries  AC  and  BC  using 

where  x  is  the  mole  fraction  of  the  composition  A  and  the  subscripts  indicate  the  materials 
of  interest.  Some  material  parameters,  however,  are  shown  in  the  experiments  to  deviate 
from  the  linear  relation  of  Eq.  (A.l)  and  have  an  approximately  quadratic  dependence  on 
the  mole  fraction  x  as 

P^fi,  _,c  =  ^Pac  +(\-x)PBc-b,x(\-x),  (A.2) 

where  b„  is  referred  to  as  the  bowing  parameter.  For  the  Cd^^Uy.^Se  and  ZnSjSc;.,  ternaries, 
the  material  parameters  such  as  the  lattice  constant  a^,  electron  effective  mass  m^,  elastic 
moduli  Cn  and  C12,  spin-orbit  splitting  energy  A,  static  dielectric  constant  e,,  background 
refractive  index  n  are  assumed  to  vary  linearly  with  the  mole  fraction  x  and,  therefore,  Eq. 
(A.l)  is  used  to  estimate  these  parameters.  On  the  other  hand,  the  parameters  such  the  bulk 
bandgap  Eg  are  obtained  using  Eq.  (A.2).  Some  other  ternary  parameters  are  available  in 
the  literatures  such  as  the  Luttinger  parameters  Yi  and  Y2  in  [Wu94b],  hydrostatic  deforma- 
tion potential  a  and  shear  deformation  potential  b  in  [Cai  95]. 
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A.2  Binary  Parameters 

The  material  parameters  for  the  relevant  11- VI  binary  alloys  used  to  estimate  the  11- 
VI  ternary  CdZnSe  and  ZnSSe  alloys  are  listed  in  Table  A.l  [Mad87,  Mad92].  The  tem- 
perature dependence  of  some  parameters  which  are  available  in  the  literatures  such  as  Eg 
are  also  listed. 

A.3  Ternary  Parameters 

In  this  section,  the  material  parameters  of  CdZnSe  and  ZnSSe  ternary  alloys,  either 
estimated  using  the  interpolation  scheme  of  the  binary  parameters  listed  in  Table  A.l  or 
quoted  from  the  literatures,  are  summarized. 

1.  Lattice  constant.  Using  the  lattice  constants  a„  of  CdSe  and  ZnSe  listed  in  Table 
A.l  inEq.  (A.l),  a„ of  Cd;^nj.^e  is  obtained 

a„(Cd^Zn;_^Se)  =  5.6676 -h  0.3844x  (A).  (A.3) 

2.  Bulk  band^ap.  The  bowing  parameters  b„  for  bulk  bandgaps  E^  of  CdZnSe  and 
ZnSSe  ternaries  are  0.22  and  0.49  eV,  respectively  [Cai95].  Plugging  these  values  for  fe„ 
and  the  values  for  Eg  and  dEJdT  of  CdSe,  ZnSe  and  ZnSe  binaries  into  Eq.  (A.2),  we  obtain 
the  temperature  dependence  of  £j  for  Cd^^n^  ..Se  and  ZnS^e^ .,  ternaries  as 

E  (Cd  Zn,      Se)   =  2.698- 1.149;c-^0.22x2 
«        ^     '-^  (eV),  (A.4) 

-  (4.5  -  0.97x)  X  10-4  (r  -  300) 

£f  ZnSSe,    J  =  2.698 -^  0.49  j:  +  0.49^^ 
*  ^     '-'^  (eV).  (A.5) 

.     ,,  .,  -(4.5  +  0.2X)  xlO^(r-300) 

where  T  is  temperature  in  K. 

3.  Spin-orbit  sphtting  energy.  The  spin-orbit  splitting  energy  A  for  Cd;^n,.,Se  is  de- 
rived using  Eq.  (A.l)  with  values  of  A  for  CdSe  and  ZnSe  in  Table  A.l  and  given  by 
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Parameters  CdSe  ZnSe  ZnS 


lattice  constant  a„  (A) 
bulk  bandgap  Eg  (eV) 

temperature  dependence  of 
bulk  bandgap  dEJdT (x  10"  eV/K) 
spin-orbit  splitting  energy  A  (eV) 
electron  effective  mass  m^  (x  m^ 
elastic  moduli  Cu  (x  10'°  N/m^) 
elastic  moduli  C,2  (x  10'°  N/m^) 
static  dielectric  constant  e^ 
background  refractive  index  n 


6.052 

5.6676 

1.771 

2.7 

3.68 

at  293  K 

at  295  K 

at  295  K 

-3.53 

-4.5 

-4.7 

0.416 

0.4 

0.112 

0.150 

0.184 

6.13 

.  9.00 

4.96 

5.43 

10.16 

8.8+1.5X  lO'T 

341« 

2.746 

Table  A.  1        List  of  material  parameters  for  II- VI  binary  alloys  CdSe,  ZnSe  and  ZnS. 


m 

A(Cd^Zn^_^Se)   =0.4  + 0.016a:  (eV).  (A.6) 

4.  Electron  effective  mass.  The  electron  effective  mass  m^  for  ternaries  Cd;,Znj.,Se 

and  ZnSjScj.,  are  estimated  using  Eq.  (A.l)  with  values  of  m^  for  CdSe,  ZnSe  and  ZnS 

mfCdZn,      Se)  =  0.15 -0.038a: 

"^  ,  (A.7) 

m^(ZnS^Se^_P   =  0.15  +  0.034a: 

where  m^  is  in  the  unit  of  electron  rest  mass  m„. 

5.  Luttinger  parameters.  The  Luttinger  parameters  Yi  and  Y2  for  Cd^^n^.^Se  and 
ZnS;^e/.,  are  given  by  [Wu94b] 

Yj  (Cd^Zn^  _^Se)   =  4.3  +  0.65a: 

Y,(Cd  Zn,      Se)   =  1.14.  +  0.22a: 

,  (A.8) 

Yi(ZnS^Se^_P   =4.3 

Y2(ZnS^Se^_P   =  1.14 

where  Yi  and  Y2  are  in  the  units  of  h^  I  (2m„).  Since  values  of  Yi  and  Y2  for  ZnS  are  not  avail- 
able in  the  literatures,  values  of  Yi  and  Y2  for  ZnSe  are  used  for  ZnS^Sc;.,  which  is  reason- 
able in  the  case  of  ZnSo.06Seo.94  where  x  =  0.06. 

6.  Elastic  moduli.  The  elastic  moduli  Cn  and  Cn  for  Cd^^n^.^Se  are  derived  using 
Eq.  (A.l)  with  values  of  Cn  and  C^  for  CdSe  and  ZnSe  in  Table  A.l  and  given  by 

C„(CdZn,    ,Se)  =   (9-2.87;c)  x  lO'" 

,„  (N/m^).  (A.9) 

Cj2(Cd^Znj_^Se)  =  (5.34-0.38;c)  x  lO^o 

7.  Static  dielectric  constant  The  temperature  dependence  of  static  dielectric  constant 
Ej  for  CdjZny.jSe  is  derived  using  Eq.  (A.l)  with  values  of  e,  for  CdSe  and  ZnSe  in  Table 
A.l  and  given  by  -  - 

e^(Cd^Zn^_^Se)  =  8.8  +  1.36x+ 1.5  x  10-3(1 -x)T.  (A.IO) 

8.  Background  refractive  index.  Plugging  the  background  refractive  indices  n  of 
CdSe  and  ZnSe  listed  in  Table  A.  1  into  Eq.  (A.  1),  n  of  Cd^n;.,Se  is  obtained 
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S(Cd^Zn;_^Se)  =  2.146 +  0.61  x.  (A.ll) 

9.  Deformation  potential.  The  hydrostatic  and  shear  components  of  the  deformation 
potential  for  Cd^Zn^.^^e  are  given  by  a  =  -5.04  eV  and  i)  =  -1.12  eV,  respectively  [Cai95]. 


t  -■ . 


APPENDIX  B 
TRANSITION  MATRIX  ELEMENT 


B.l  Polarization-Dependent  Effects 


FromEqs.  (57a-c)  in  [Coi93],  the  polarization-dependent  of  the  transition  matrix  el- 
ements for  C-HH  and  C-LH  transitions  are  given  by 

if^  n        f  1  -  |k  •  ep      for  C-HH  transition 

I     T\     _    J  CB  1^ 

|A/p   ~   U-nlkeP       for  C-LH  transition  ' 

where  k  is  a  unit  vector  directed  along  the  electron  wavevector  k  and  e  is  the  unit  polar- 
ization vector.  Fig.  B.l(a)  and  (b)  show  plots  of  k  and  e  in  the  three-dimensional  (3D) 
rectangular  coordinates  and,  therefore,  k  and  e  can  be  expressed  as 

k  =  sin9cos<px-H  sin9sin(py -H  cos6z,  (B.2) 

e=  sinO'coscp'x -H  sinO'sintp'y -H  cosS'z.  (B.3) 

where  it  and  y  are  chosen  as  the  unit  vectors  in  the  in-plane  directions  (in  the  plane  of  the 

QW),  z  is  chosen  as  the  unit  vector  in  the  growth  direction  (normal  to  the  plane  of  the  QW), 

9  and  6'  are  angles  measured  from  the  z-axis  while  (p  and  (p'  are  angles  measured  from 

the  X-axis.  Using  Eqs.  (B.2)  and  (B.3),  we  obtain 

|k  •  er  =  (sin0cos9sin9'cos(p' +  sin9sin(psine'sin(p'-i- cosScose')^ 
=  sin^9cos^(psin^9'cos^(p' 
+  sin^9sin^c()sin^9'sin2<p' 

-I-  cos2ecos29'  •  (B-4) 

+  2  sin^  9  sin  (p  cos  cp  sin^  9'  sin  cp'  cos  (p' 
+  2sin9cos9sin(psin9'cos9'sin(p' 
-t- 2  sin  9  cos  9  cos  (p  sin  9' cos  9' cos  (p' 


tea 
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Since  photons  of  a  given  polarization  interact  with  electrons  of  k  vectors  pointing  in  differ- 
ent directions,  |k  •  ep  in  Eq.  (B.4)  needs  be  averaged  over  these  interactions  by  a  proper 
average  over  k  vector.  Fig.  B.l(c)  shows  a  plot  of  k  in  the  3D  rectangular  coordinates 
where  kj^  and  kf,  are  the  components  of  k  vector  in  the  growth  and  in-plane  directions.  The 
j  subscript  indicates  theyth  quantized  subband  since  electrons  are  quantized  in  the  growth 
direction  (z  direction).  For  electrons  in  a  given  subband,  kj±  is  calculated  from  the  associ- 
ated quantized  energy  level  using  Eq.  (3.9).  At  the  subband  edge,  k,,  =  0  and  9^  =  0.  As 
moving  away  from  the  subband  edge,  k  vector  with  a  given  ki,  sweeps  over  a  cone  angle  of 
2Qj  with  respect  to  the  z-axis.  As  a  consequence,  we  average  Eq.  (B.4)  over  9  from  0  to 
2jc  to  obtain 


|k' 

■e 

=  -sin^eysin^e'cos^tp' 

+  -sin^9  sin^e'sin^(p' 

-1-  cos^ 6  COSTS' 
-hO 
+  0 
+  0 

(B.5) 


where  8  is  replaced  with  9^  indicating  theyth  subband  transition  pair.  For  a  TE  polarization, 
e  =  X  or  y  and  therefore  plugging  9'  =  -  and  ip'  =  0  or  -  in  Eq.  (B.5),  we  have 

lit  •  ep  =  Tsin^9;  =  ---cos^9  for  TE  polarization .  (B.6) 

For  a  TM  polarization,  e  =  z  and  therefore  plugging  9'  =  0  in  Eq.  (B.5),  we  have 

|kep  =  cos29.  for  TM  polarization .  (B.7) 

For  spontaneous  emission,  e  is  arbitrarily  polarized.  In  that  case,  we  need  to  average  Eq. 
(B.5)  over  9'  from  0  to  Jc/2  and  p'  from  0  to  2n  in  spherical  coordinates  and  the  result  is 
given  by 

I'        19  1 

|k  ■  er  =  T  for  averaged  polarization  (B.8) 
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(«) 


(b) 


(c) 


Figure  B.  1  The  three  dimensional  plots  of  various  vectors  in  the  rectangular  coordi- 
nates, (a)  the  unit  vector  k  directed  along  the  electron  wavevector  k;  (b)  the 
unit  polarization  vector  e ;  and  (c)  the  electron  wavevector  k. 
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Using  Eqs.  (B.6)  and  (B.7)  in  Eq.  (B.l),  the  transition  matrix  elements  \Mj.\  ^  for  the  optical 
gain  with  TE  and/or  TM  polarizations  can  be  expressed  as  a  function  of  6^ 

'1  +  1  cos^e^       for  C-HH  TE  gain 

1  -cos^  e^  for  C-HH  TM  gain 

5   1,  •  (»•') 

2-i  cos^  e^.        for  C-LH  TE  gain 

1  +  cos2  e^         for  C-LH  TM  gain 

Using  Eq.  (B.8)  in  Eq.  (B.l),  the  average  transition  matrix  \Mj\^  for  the  spontaneous  emis- 
sion can  be  also  be  expressed  as  a  function  of  Qj 

|mJ2       2 

'      '     =  -     for  both  C-HH  and  C-LH  spontaneous  emission .  (B.IO) 

Eqs.  (2.27)  and  (2.32)  are  therefore  given  as  Eqs.  (B.9)  and  (B.IO),  respectively. 


B.2  Angle-Dependent  Effects 

In  the  former  section,  we  have  derived  the  expression  for  IWj-P  as  a  function  of  By 
without  knowing  how  to  calculate  9y.  From  Fig.  B.  1  (c),  we  use  the  definition  of  the  cosine 
function  of  angle  6^  for  a  right  triangle  to  obtain 

Multiplying  both  the  denominator  and  numeration  at  the  right-hand  side  (RHS)  of  Eq. 
(B.l  1)  by  SV(2m,yii),  it  becomes 

where  m,f^  is  the  in-plane  reduced  effective  mass  for  theyth  subband  transition  pair  defined 
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in  Eq.  (2.31).  Recalling  Eqs.  (2.21)  and  (2.22),  the  transition  energy  £„  can  be  rewritten  as 

'  ^  ,  (B.13) 

■      =  E  +E  ■  +  £  .  + 

i       ^J       ^'     2m,^„ 

where  definition  for  mrf\  in  Eq.  (2.31)  is  used  in  the  last  expression  for  Eq.  (B.13).  Rear- 
ranging Eq.  (B.13),  we  obtain 

=  E    -E  -E    -E  ..  (B.14) 

2m,y||  cv  g  CJ  v; 

Recalling  the  relationship  between  the  ;'th  quantized  energy  level  £„;  and  the  associated 
equivalent  growth-direction  wavevector  ;fc„^x  from  Eq.  (3.11),  the  summation  of  the  y'th 
quantized  energies  in  the  CB  and  VB  are  given  by 

»^*?/i     »^*v;i     »^*?± 

£.  +  £.  =  EZi  +  ^its-— i±,  (B.15) 

^1       ^1       Im^.^     Im^^^     Im^.^ 

where  the  last  expression  is  approximated  under  a  typical  case  that  k^^  and  k^j^  for  theyth 
CB  and  VB  subbands  are  approximately  equal.  Plugging  Eq.  (B.15)  into  the  following 
equation,  we  obtain 

ir^  =  :^--^  =  (E.  +  E.)-^.  (B.16) 

We  then  use  Eqs.  (B.14)  and  (B.16)  in  Eq.  (B.12)  to  obtain 

(£.  +  £.)  -^ 

cos^e,.  =  '-^ .  (B.17) 

(E    -E  -E    -E  ■)  +  (E  ■  +  £  ^-^ 

Eqs.  (2.9)  or  (B.17)  is  therefore  derived  for  the  anisotropic  effective-mass  model  described 
in  Chapter  2.  For  an  isotropic  effective-mass  model  where  the  in-plane  effective  masses 
are  equal  to  their  growth-direction  counterparts  and  m  n  =  m-^ ,  Eq.  (B.  17)  reduces  to 
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cos^e,.  =  ^r^'.  (B.18) 

Eq.  (B.18)  is  usually  used  in  the  literatures  [Chi88]  where  the  anisotropy  in  effective  mass- 
es is  ignored. 
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